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1. 


INTRODUCTION 


A prevalent  problem  encountered  by  Army  commanders 
is  the  lack  of  precise  information  on  the  location  of  their 
troops  and  equipment  as  well  as  that  of  adversary  in  a 
battlefield  environment.  The  projected  NAVSTAR  Global  Posi- 
tioning System  (GPS)  is  capable  of  providing  this  location 
information  accurately  through  use  of  passive  receiving 
equipment  by  the  Army  troops.  However,  as  depicted  in  Fig. 
1-1,  there  are  several  potential  problems  associated  with 
the  use  of  GPS  in  a battlefield  environment.  These  problems 
include  attenuation  of  the  navigation  signal  by  heavy  foliage, 
intentional  jamming  of  the  navigation  signal  by  the  adver- 
sary and  adverse  environmental  conditions  affecting  user 
equipment  operation  and  navigation  accuracy. 


Figure  1-1  Army  GPS  Users  in  a Jamming  Environment 


This  report  addresses  these  problems.  Included  are 
assessments  of  the  degradation  in  user  location  accuracy  that 
may  be  experienced  in  a battlefield  environment  as  well  as 
actions  which  may  be  taken  to  alleviate  or  counteract  these 
problems  when  encountered.  The  specific  goals  of  this  study 
are  to : 

• Assess  the  impact  on  the  Army  manpack 
GPS  user  navigation  performance  of  the 
use  of  directional  receiver  antennas. 

• Assess  the  vulnerability  of  the  Army 
GPS  user  to  electronic  warfare  and 
consider  methods  of  countering  this 
threat . 

• Assess  the  impact  on  Army  airborne 
GPS  user  navigation  performance  of 
the  receiver  clock  crystal  oscillator 
errors  induced  by  the  environment. 

This  study  is  an  extension  of  a previous  study  of  the  Army 
GPS  user  performance  described  in  Ref.  1.  This  reference 
is  frequently  cited  in  the  present  work  and,  wherever  pos- 
sible, the  same  notation  is  used. 

This  study  uses  analysis  tools  developed  under 
previous  studies.  In  particular,  the  previously  developed 
GPS  Direct  Simulation  (DSIM)  computer  program  (Ref.  2)  was 
used  here  for  the  assessment  of  user  navigation  performance 
in  both  the  directional  antenna  and  environmental  clock  mod- 
eling studies. 


1 . 1 BACKGROUND 

The  development  of  a high-accuracy,  world-wide, 
24-hour  navigation  satellite  system  has  been  in  progress  for 
over  a decade.  Extensive  studies  and  tests  have  proven  the 
feasibility  of  the  concept  and  that  it  could  be  configured 
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within  the  capabilities  of  current  technologies.  In  1973, 
TASC  initiated  a study  of  the  performance  of  Army  users  of 
the  GPS  (Ref.  1).  That  effort  consisted  of  the  evaluation 
of  the  relevance  of  previous  Air  Force  navigation  satellite 
test  programs  to  Army  users,  an  analysis  of  the  performance 
of  selected  Army  users  of  the  GPS,  and  a study  of  the  impact 
of  the  GPS  signal  structure  on  Army  users  including  consid- 
eration of  vulnerability  to  electronic  warfare. 

The  present  study  is  an  extension  of  the  previous 
Army  user  study  with  consideration  of  the  impact  on  navi- 
gation performance  of  directional  antennas  and  adverse  en- 
vironments. The  vulnerability  study  is  extended  with  a 
more  complete  model  of  the  propagation  environment  and  con- 
sideration of  possible  countering  techniques.  The  investi- 
gation of  the  environmental  effects  on  receiver  crystal 
oscillators  is  a new  topic. 


1.2  THE  NAVSTAR  GLOBAL  POSITIONING  SYSTEM 

The  NAVSTAR  Global  Positioning  System  is  a satellite- 
based  passive  radio  navigation  system  that , in  its  final  opera- 
tional phase,  will  be  capable  of  providing  highly  accurate 
navigation  information  to  properly  equipped  users  on  a con- 
tinuous world-wide  basis.  The  principal  elements  of  the 
system  are  a set  of  satellites  that  transmit  navigation 
signals,  a network  of  ground  monitoring  stations  and  the 
various  users.  The  GPS  is  to  evolve  in  three  phases  with 
the  final  phase  (Phase  III)  consisting  of  24  satellites 
in  a 3 x 8 (three  planes,  eight  satellites  per  plane) 
constellation  with  each  satellite  in  a 12-hour  (altitude 
of  3.16  earth  radii)  orbit  inclined  63  deg.  The  satellite- 
to-satellite  in-plane  spacing  is  uniform  (45  deg  spacing) 


with  a 15  deg  inter-plane  phasing.  The  plane-to-plane  angu- 
lar spacing  is  120  deg. 

The  ground  monitoring  system  tracks  the  satellites 
and  computes  the  ephemerides  of  each  satellite.  These 
satellite  ephemerides  and  related  data  are  then  "uploaded" 
to  each  satellite  once  a day  for  subsequent  transmission  to 
the  users. 

The  satellites  transmit  coded  navigation  signals 
which  enable  a user  to  make  range  measurements  between  him- 
self and  the  satellite.  These  measurements  are  made  by 
comparing  the  time  of  arrival  of  the  navigation  signal  with 
an  assumed  time  of  transmission  as  determined  by  the  user's 
reference  clock.  This  time  difference  is  converted  into  a 
range  measurement  by  means  of  a known  propagation  constant. 
The  user's  reference  clock  is  not  perfectly  synchronized 
with  the  clock  in  the  satellite  that  is  used  to  generate 
the  navigation  signal.  Thus,  this  range  measurement  dif- 
fers from  the  geometrical  range  between  user  and  satellite 
by  the  relative  offset  or  bias  between  the  user  and  satel- 
lite clock.  The  measured  range  is  called  "pseudo-range" 
to  denote  this  fact. 

The  user  obtains  a three-dimensional  position  fix 
by  making  at  least  four  independent  pseudo-range  measure- 
ments (to  four  separate  satellites).  These  four  measure- 
ments permit  the  simultaneous  determination  of  the  three 
unknown  position  coordinates  and  the  unknown  clock  offset. 

In  an  analogous  manner,  velocity  estimates  are  obtained 
from  pseudo-range-rate  measurements  by  extracting  the  dop- 
pler  shift  of  the  navigation  signal  carrier  frequency.  In 
order  for  the  user  to  convert  these  measurements  to  position 
and  velocity  estimates,  he  must  know  the  position  and  vel- 
ocity of  those  satellites  originating  the  navigation  signals 


used  in  the  measurements.  This  information  is  sent  by  the 
satellites  to  the  users  via  the  navigation  signals. 


1.3  REPORT  OVERVIEW 

The  major  results  of  this  study  are  contained  in 
Chapters  2,  3 and  4.  The  appendices  contain  additional  de- 
tails of  the  mathematical  models  and  analysis  techniques 
used  in  the  study. 

The  assessment  of  the  impact  of  directional  re- 
ceiver antennas  on  Army  user  navigation  performance  is  pre- 
sented in  Chapter  2.  The  number  of  satellites  visible  with 
directional  antennas  of  different  beamwidths  is  discussed 
first.  Then  performance  projections  are  presented  for 
users  employing  directional  antennas.  Next,  the  design 
of  directional  antennas  using  an  array  of  discrete  antenna 
elements  is  discussed.  The  optimum  method  of  eliminating 
directional  interference  is  presented  along  with  theoreti- 
cal performance  projections.  Applicability  to  the  GPS  man- 
pack  user  is  discussed.  Finally,  some  antenna  deployment 
procedures  are  considered. 

The  results  of  the  vulnerability  study  are  pre- 
sented in  Chapter  3.  The  threat  model  and  scenario  are 
outlined  with  a brief  discussion  of  the  propagation  loss 
model.  Threat  countering  techniques  are  then  considered 
and  their  effectiveness  evaluated.  Finally,  an  evaluation 
of  Army  GPS  user  vulnerability  to  electronic  warfare  is 
presented  for  the  threat  model  and  scenario  previously 
outlined.  An  assessment  of  the  relative  impact  of  the 
threat  countering  techniques  is  also  made. 

The  effects  of  the  environment  on  the  user  re- 
ceiver, in  particular,  the  receiver  crystal  oscillator  is 


discussed  in  Chapter  4.  The  environmental  clock  model  is 
presented  along  with  user  performance  projections  for  us- 
ers employing  this  clock. 

The  report  is  summarized  and  some  conclusions  pre- 
sented in  Chapter  5.  Several  suggestions  and  recommenda- 
tions, based  on  insights  gained  while  generating  the  data 
of  the  previous  chapters,  are  presented . 

Characteristics  of  directional  antennas  are  con- 
sidered in  Appendix  A with  emphasis  placed  on  parabolic-type, 
antennas.  The  GPS  Phase  III  configuration  and  performance 
data  used  in  the  user  performance  projections  are  discussed 
in  Appendix  B.  The  error  modeling  technique  used  in  devel- 
oping a model  for  intrinsic  clock  stability  is  discussed  in 
Appendix  C.  Details  of  the  propagation  loss  model  used  in 
the  vulnerability  evaluation  study  are  presented  in  Appen- 
dix  D. 
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2. 


DIRECTIONAL  RECEIVER  ANTENNAS 


GPS  users  may  use  directional  antennas  in  order 
to  overcome  some  of  the  problems  which  could  be  encountered 
in  a battlefield  environment.  In  particular,  they  can  be 
used  to  increase  the  received  signal  strength  when  operating 
under  heavy  foliage  and  to  reduce  interfering  signals  such 
as  intentional  signal  jamming.  However,  directional  antennas 
restrict  the  user's  visibility  to  satellites  and  thus  may  de- 
grade his  navigation  accuracy  and/or  increase  the  time  re- 
quired to  obtain  high  navigational  accuracy. 

This  chapter  considers  the  use  of  directional  re- 
ceiving antennas  by  GPS  users.  A discussion  of  the  number 
of  satellites  visible  with  directional  antennas  of  various 
beamwidths  and  orientations  is  presented  in  Section  2.1. 
Performance  projections  for  users  with  directional  antennas 
are  presented  in  Section  2.2.  The  results  presented  in  this 
section  illustrate  the  advantages  of  using  Kalman  filters  to 
perform  the  processing  for  computing  user  position.  Null 
and  beam  steering  techniques  are  discussed  in  Section  2.3. 

A general  theoretical  discussion  of  electronic  steering  is 
presented  as  well  as  an  illustration  of  the  performance 
that  can  be  expected  with  simple  antenna  arrays.  Finally, 
in  Section  2.4,  antenna  deployment  procedures  are  discussed 
for  typical  manpack  GPS  user  scenarios.  Discussion  of  the 
directional  antenna  characteristics  used  in  this  chapter 
and  the  equations  relating  the  various  antenna  parameters 
are  presented  in  Appendix  A. 
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2.1 


SATELLITE  VISIBILITY  WITH  DIRECTIONAL  ANTENNAS 


The  satellites  that  are  visible  to  a GPS  user  depend 
not  only  upon  his  location  and  antenna  orientation,  but  also 
upon  the  satellite  configuration.  This  study  uses  the  GPS 
Phase  III  satellite  configuration  that  was  discussed  briefly 
in  Section  1.2.  Further  details  of  this  satellite  configura- 
tion are  presented  in  Appendix  B. 

The  number  of  satellites  visible  to  users  with 
vertical  boresight  (or  main  lobe)  directional  antennas  for 
various  beamwidths  is  illustrated  in  Section  2.1.1.  These 
satellite  visibility  results  are  extended  to  off-vertical 
boresight  directional  antennas  in  Section  2.1.2. 

2.1.1  SATELLITE  VISIBILITY  WITH  VERTICAL 
BORESIGHT  ANTENNAS 

The  number  of  satellites  visible  to  a GPS  user  with 
a directional  receiving  antenna  are  determined  by  his  longi- 
tude and  latitude,  antenna  orientation,  and  the  time  of  day. 
Also,  as  the  user  makes  measurements,  visible  satellites 
may  drop  from  view  and/or  new  satellites  appear.  In  order 
to  reduce  the  number  of  variables  involved  but  still  pro- 
vide some  idea  of  satellite  visibility,  the  average  number 
of  satellites  visible  to  users  at  different  latitudes  with 
vertical  boresight  directional  antennas  was  computed.  This 
was  accomplished  by  averaging  over  a 24-hour  period  the 
number  of  satellites  that  a user  could  "see"  at  different 
latitudes  with  different  directional  antenna  beamwidths. 

The  results  of  this  computation  are  illustrated  in  Fig.  2.1-1. 
These  results  indicate  that  the  antenna  beamwidth  must  be  ap- 
proximately 140  deg  (or  greater)  in  order  to  insure  the  vis- 
ibility of  four  satellites  simultaneously.  This  confirms 
similar  results  reported  in  Ref.  1.  For  users  at  latitudes 
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Figure  2.1-1 
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Satellite  Visibility  of  Vertical 
Boresight  Directional  Antennas 


greater  than  25  deg,  a 120  deg  beamwidth  directional  antenna 
could  be  deployed  and  still  retain  visibility  to  four  satel- 
lites . 

As  can  be  expected,  the  narrower  the  beamwidth,  the 
fewer  the  satellites.  Thus,  a beamwidth  of  90  deg  usually 
results  in  fewer  than  three  visible  satellites.  However, 
as  observed  in  Ref.  1,  a user  can  still  obtain  a position 


fix  with  fewer  than  four  satellites  if  he  is  willing  to  wait 


for  a significant  change  in  satellite 


These  results  are  averages  over  a 24-hour  period; 
thus,  there  may  be  times  and  locations  where  four  satellites 
are  visible  with  narrower  beam  (e.g.,  less  than  120  deg) 
directional  antennas.  However , -since  these  regions  and 
times  are  isolated,  they  are  difficult  to  rely  on  when 
considering  a continuous,  24-hour  navigation  capability. 


9 


Thus,  these  average  results  provide  a realistic  picture  of 
satellite  visibility  with  directional  antennas. 

2.1.2  Visibility  With  Off-Vertical 
Boreslght  Antennas 

The  visibility  results  of  the  preceding  section  can 
be  extended  to  an  off-vertical  boreslght  directional  antenna 
by  referencing  Its  beamwidth  and  elevation  angle  to  a suit- 
ably located  vertical  boreslght  directional  antenna.  This 
enables  the  evaluation  of  any  potential  advantages  of  off- 
vertical  antenna  orientations  using  the  results  of  Section 
211.  The  relationship  between  vertical  boresight  (refer- 
ence) and  off-vertical  boresight  (user)  directional  antennas 
with  the  same  satellite  visibility  or  satellite  orbit  cover- 
age is  illustrated  in  Fig.  2.1-2.  This  figure  shows  that 
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Figure  2.1-2  Antenna  Beamwidth  Comparisons 


the  "sky  coverage"  of  an  off-vertical  (user)  antenna  can  be 
adjusted  so  that  it  is  the  same  as  that  of  a vertical  (ref- 
erence) antenna  of  different  beamwidth  located  directly 
under  the  "sky  coverage"  area. 

As  the  user  antenna  elevation  angle  decreases,  the 
latitude  difference  between  the  user  and  reference  antennas 
increases  and  the  user  antenna  beamwidth  decreases  slightly 
in  order  to  maintain  the  same  satellite  visibility  or  "sky 
coverage"  as  the  reference  antenna.  The  relation  between 
the  user  and  reference  antenna  beamwidths  is  shown  in  Fig. 
2.1-3.  This  figure  indicates  that  the  user  antenna  beam- 
width  does  not  change  appreciably  as  the  latitude  difference 
increases.  For  a given  reference  antenna  beamwidth,  a point 
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Figure  2.1-3 


Vertical  and  Off-Vertical  Boresight 
Directional  Antenna  Relationship 


will  be  reached  where  the  horizon  interferes  with  the  user 
antenna  main  lobe  or  beam  as  the  latitude  difference  in- 
creases; this  point  is  indicated  by  the  dashed  line  in  Fig. 
2.1-3. 


Because  the  equivalent  beamwidth  of  an  off-vertical 
boresight  (user)  directional  antenna  does  not  change  apprec- 
iably (within  the  applicable  elevation  angles)  from  that  of 
a vertical  reference  directional  antenna,  the  satellites 
visible  to  an  off-vertical  boresight  antenna  are  approxi- 
mately the  same  as  those  visible  to  a suitably-located 
vertical  boresight  antenna  of  the  same  beamwidth.  Thus, 
any  advantage  gained  by  using  off-vertical  antennas  would  be 
derived  by  pointing  to  an  advantageous  sector  of  the  sky. 
However,  Fig.  2.1-1  indicates  only  a small  advantage  in  point- 
ing the  user  antenna  to  more  advantageous  "sky  latitudes" 
(i.e.,  the  curves  are  relatively  flat).  It  thus  appears  that 
little  is  gained  by  the  deployment  of  off-vertical  boresight 
directional  antennas  for  the  simultaneous  reception  of  four 
GPS  navigation  signals. 


2.2  USER  PERFORMANCE  WITH  DIRECTIONAL  ANTENNAS 

As  illustrated  in  the  preceeding  section,  the  GPS 
users'  satellite  visibility  is  restricted  with  directional 
antennas.  Depending  upon  the  beamwidth  of  these  antennas,  the 
users'  visibility  may  be  restricted  to  four  satellites  that 
provide  a relatively  large  GDOP  (Geometric  dilution  of  Preci- 
sion, Ref.  1)  or  to  fewer  than  four  satellites.  Since  the 
satellite-user  geometry  (as  reflected  in  the  GDOP)  directly 
affects  the  accuracy  with  which  a user  may  determine  his  posi- 
tion, a large  GDOP  may  result  in  poor  estimates  of  user  posi- 
tion. If  fewer  than  four  satellites  are  visible  with  a parti- 
cular directional  antenna  (e.g.,  with  a very  narrow  beamwidth), 


12 


the  user  may  have  to  wait  for  a significant  change  in  the 
position  of 'the  satellites  to  obtain  a navigation  fix  (Ref.  1) 
or,  alternatively,  point  the  antenna  separately  to  four  satel- 
lites . 


This  section  considers  the  navigation  performance 
of  a user  employing  directional  antennas.  Two  directional 
antennas  are  considered,  one  which  restricts  visibility  to 
four  satellites  and  one  which  restricts  visibility  to  a 
single  satellite.  The  scenario  and  models  used  in  this  study 
are  discussed  in  Section  2.2-1.  The  satellite-user  geometry 
effects  or  GDOP  for  this  scenario  are  discussed  in  Section 
2.2.2  along  with  the  satellite  visibility  data.  User  per- 
formance projections  are  then  presented  in  Section  2.2.3 
for  three  user  antenna  beamwidths.  These  performance  pro- 
jections (as  well  as  those  of  Chapter  4)  were  obtained  us- 
ing the  DSIM  computer  program  (Ref.  2)  and  the  GPS  tracking 
system  performance  data  summarized  in  Appendix  B.  DSIM 
consists  of  a "real  world"  simulator  as  well  as  a simula- 
tion of  the  user  navigation  filter.  The  "real  world" 
simulator  generates  actual  satellite  ephemerides,  synthe- 
sizes noise-corrupted  pseudo-range  (and  pseudo-range-rate, 
if  required)  measurements  and  generates  the  user  trajectory, 
thus  making  possible  an  assessment  of  the  performance  of 
the  user  navigation  filter  implementation  by  providing  time 
histories  of  estimated  navigation  states  (e,g.,  position 
and  velocity)  for  comparison  with  the  true  values.  The 
real  world  models  used  in  the  DSIM  program  are  discussed 
in  Ref.  2,  except  for  the  satellite  clock  model,  which  is 
outlined  in  Appendix  B,  and  the  user  clock  model,  which 
is  discussed  in  Section  2.2-1. 
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2.2.1  User  Scenario  and  Models 
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Scenario  — The  scenario  assumed  for  this  investiga- 
tion consists  of  a stationary  manpack  user  with  one  of  three 
antennas:  conventional  (baseline),  wide  beam  directional  and 

narrow  beam  directional.  The  conventional  antenna  is  a fixed 
hemispherical  antenna  with  170  deg  beamwidth  (5  deg  masking 
angle)  directed  vertical.  The  wide  beam  directional  antenna 
has  a beamwidth  of  120  deg  directed  vertically  (main  beam). 
With  these  two  antennas,  measurements  are  made  sequentially  to 
the  four  satellites  (nos.  1,  2,  3,  5,  which  provide  the  best 
GDOP)  at  a 10  sec/satellite  rate.  The  narrow  beam  directional 
antenna  is  assumed  to  provide  visibility  to  only  one  satellite, 
and  is  slewed  sequentially  by  the  user  to  selected  satellites. 
This  antenna  dwells  on  each  satellite  for  50  seconds  during 
which  time  measurements  are  made  at  a rate  of  10  sec/measure- 
ment. The  direction  in  which  this  narrow  beam  antenna  is  to 
be  pointed  (which  depends  upon  the  selected  satellite)  is  de- 
termined by  the  users'  GPS  computer.  The  user  for  the  simula- 
tion runs  presented  here  is  assumed  to  be  located  at  50  deg 
north,  100  deg  west. 


Receiver  — For  purposes  of  this  study,  a manpack  re- 
ceiver was  modeled  as  a sequential  set  utilizing  pseudo- 
range measurements  (no  pseudo  range  rate).  These  measurements 
are  corrupted  by  several  error  sources  including  receiver 
mechanization  errors,  signal  propagation  path  delays  and  user 
clock  errors,  in  addition  to  the  satellite  epheineris  and  clock 
errors.  The  receiver  mechanization  errors  modeled  in  this 
study  include  thermal  noise  and  signal  quantization.  Thermal 
noise  errors  are  modeled  as  a normally  distributed,  zero-mean, 
random  sequence  while  the  signal  quantization  errors  are  uni- 
formly distributed.  The  parameter  values  associated  with  the 
receiver  mechanization  error  sources,  summarized  in  Table  2.2-1, 
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TABLE  2.2-1 


GPS  RECEIVER  MECHANIZATION  ERRORS 


ERROR  SOURCE 

STANDARD  DEVIATION 
OR  RANGE 

Pseudo-Range  Thermal  Noise 

1.4  m 

(Normally  Distributed) 

Pseudo-Range  Quantization 

±0.3  m 

(Uniformly  Distributed) 

are  based  on  previous  studies  (Ref.  2 and  31),  concurrent 
studies  (Ref.  32)  and  Air  Force  Specification  (Ref.  24). 

The  propagation  effects  modeled  in  this  study  include 
residual  ionospheric  and  tropospheric  delays.  The  residual 
ionospheric  delay  errors  are  those  remaining  after  two- 
frequency  ionospheric  compensation  (Ref.  2).  Residual  tropo- 
spheric delay  errors  are  those  remaining  after  analytical  com- 
pensation. The  residual  tropospheric  delay  errors  are  most 
significant  near  the  horizon  and  thus  are  a function  of  the 
user-to-satellite  elevation  angle.  The  numerical  values  used 
in  this  study  for  the  propagation  path  errors  have  been  ob- 
tained from  previous  studies  (Ref.  2 and  31)  and  are  summarized 
in  Table  2.2-2. 


TABLE  2.2-2 

PROPAGATION  PATH  DELAY  ERRORS 


STANDARD 

CORRELATION 

ERROR  SOURCE 

DEVIATION 

TIME 

Residual  Ionospheric  Delay 

0.76  m 

7200  sec 

Residual  Tropospheric  Delay* 

0.076  m 

00 

♦This  value  is  multiplied  by  the  cosecant  of  the  user-to- 
satellite  elevation  angle. 


User  Clock  — This  clock  is  a quartz  crystal  oscil- 
lator with  the  intrinsic  stability  model  developed  in  Appen- 
dix C.  The  fractional  frequency  stability  plot  for  this 
clock  model  is  shown  in  Fig.  2.2-1.  This  clock  error  model 
as  depicted  in  Fig.  2.2-2,  can  be  written  in  state  variable 
form  (Ref.  5)  as 

*c(t)  = Fc  ^c(t)  + *c(t)  (2.2-1) 

where 

xc(t)  = l <5<j> , 6f  , 6 f , xlf,x2flTis  the 
5x1  vector  of  error  states 
defined  in  Table  2,2-3. 

F is  the  clock  error  model  dynamics 

matrix 

wc(t)  = [ w^  0 0 w2  W3]  T is  the  5x1 

vector  of  white  noise  such 
that  E l wc ( t ) ] = 0,  Cov  [w  (t)]  = 

E[wc(t)  w T(X)]  . Qc( t ) 6(t-x) 
and  6(t)  is  the  Dirac  delta  func- 
tion . 


f- AVERAGING  TIME  (wc) 


Figure  2.2-1  Fractional  Frequency  Stability 

Characteristic  for  Quartz  Crystal 
Clock  Model 
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Figure  2.2-2  User  Clock  Model 
TABLE  2.2-3 

USER  CLOCK  ERROR  STATE  DEFINITIONS 
SYMBOL  ERROR  SOURCE  DEFINITION 


Phase  Error 
Frequency  Error 
Aging  Error 
Flicker  Noise  State  1 

Flicker  Noise  State  2 


m/sec' 


The  clock  error  model  dynamics  matrix,  F , is  given  by 


F 

c 


0 1 0 f1  0 
0 0 1 0 0 
ooooo 

OOOfg1 
0 0 0 f3  o 


(2.2-2) 


where  f^,  f2  and  f3  are  defined  in  Eq.  (C.3-7).  The  clock 
error  model  process  noise  matrix,  Q , is  given  by 


Q 


c 


1.026 
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CO 
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9.234 
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10-6 
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9.234 
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1.005 
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(2.2-3) 

Units  of  the  Qc~matrix  elements  are  consistent  with  those 
given  in  Table  2.2-3  for  the  corresponding  state  variables. 


User  Filter  — The  user  navigation  filter  is  a Kalman n 


mechanization  (Ref.  5)  which  assumes  the  user  is  stationary. 
The  user  filter  error  model  in  state  variable  form  can  be 
written  as 


x(t ) = F x( t ) + w(t ) 


(2.2-4) 


where 


x(t).  = 


xT 

_c 


is  the  8x1  column  vector 
of  user  error  states 


is  the  user  error  dynamics  8x8  matrix 
18 


F 


w(t  ) 


= [0,0,0,  w£]  is 

vector  of  white 


the  8x1  column 
process  noise 


such  that  E[w(t)l  = Q, 

Cov  [w(  t )]  = F[w(  t ) wT(  t )]  = 
Q(t)  S(t-t) 


The  user  position  error  states  are  given  by 

T 


6 R = 6R  , 6R  , 6h 
— u L n ’ e ’ J 


(2.2-5) 


where  <5R  , 6R  and  6h  are  the  north,  east  and  altitude  com- 
n e 

ponents.  The  user  error  dynamics  matrix,  F,  is  given  by 


“ 

“ 

°3x3 

°3x5 

°5x3 

F 

c 

(2.2-6) 


where  F is  given  by  Eq.  (2.2-2).  The  user  process  noise 
c 

matrix,  Q,  is  given  by 


| 

Q = 

°3x3  ' 

°3x5 

°5x3 

_ 

where  Q is  given  by  Eq.  (2.2-3). 

The  inputs  to  the  filter  are  the  pseudo-range  mea- 
surements from  the  receiver.  These  measurements  were  assumed, 
by  the  filter,  to  contain  1.6  m white  measurement  error.  This 
value  is  assumed  to  account  for  propagation  delays,  receiver 
errors,  and  satellite  clock  and  ephemeris  errors.  The  filter 
uses  these  measurements  along  with  the  initial  values  summar- 
ized in  Table  2.2-4  to  estimate  the  user's  position  (both  hori- 
zontal and  altitude)  and  the  five  states  of  the  clock  model. 


During  the  course  of  this  study,  it  was  noted  that 
the  originally  developed  manpack  navigation  filter  summarized 
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in  the  preceding  paragraphs,  could  produce  significant  errors 
due  to  the  unmodeled  errors  in  the  satellite  ephemeris  and 
clock  data.  Two  methods  of  augmenting  user  filters  to  ac- 
count for  unmodeled  errors  were  investigated  during  this 
study.  One  method  uses  process  noise  on  the  position  states, 
while  the  other  involved  the  addition  of  "pseudo-states"  to 
the  user  filter.  The  process  noise  method  of  filter  augmen- 
tation effectively  accounts  for  the  satellite  ephemeris  and 
clock  errors  as  uncertainty  in  user  position.  This  method 
changes  the  user  filter  process  noise  matrix,  Eq.  (2.2-7), 
as  follows: 


(2.2-8) 


where  q is  the  spectral  density  of  the  user  position  un- 
certainty. The  value  of  q is  best  determined  by  simula- 
tion experiments. 

The  pseudo-state  method  of  filter  augmentation  adds 
to  the  user  filter  model  a pseudo-state  for  each  satellite  to  ef 
fectively  represents  the  combined  ephemeris  and  clock  errors 
for  that  satellite.  The  user  position  error  states,  6Ru, 
are  augmented  with  pseudo-states,  £,  as  follows  (assuming 
m satellites  are  used  for  measurements): 


6R 

-u 

E 


(2.2-9) 


6R  is  the  3x1  column  vector  of  user 
u position  error  states 

T 

p = [ p i , p2,  pm]  is  the  mxl  column 

vector  of  pseudo-states 

Those  pseudo-states  are  modeled  as  biases  in  the  user  filter 
and  are  dynamically  uncoupled  from  the  other  filter  states. 
The  user  filter  error  dynamics  matrix,  Eq . (2.2-6),  is  thus 
modified  as  follows, 


°( 3+m  ) x (3+m) 

• °(3+m)  x 5 

°5  x (3+m) 

i p 

c 

The  user  filter  process  noise  matrix,  Eq.  (2.2-6), 
fied  as  follows, 


(2.2-10) 


is  modi- 


Q 


(3+m)  x (3+m) 


(3+m)  x 5 


5 x (3+m) 


Q, 


(2.2-11) 


Since  the  pseudo-states  are  used  to  represent  the  combined 
satellite  ephemeris  and  clock  errors  which  occur  in  the 
pseudo-range  measurements,  they  are  used  during  the  filter 
updates  with  these  measurements;  that  is,  if  the  true  pseudo- 
range, fiz^,  to  the  itn  satellite  is  given  by  (Ref.  1) 

6z.  = lTi  5R  +60+  n,  (2.2-12) 

i -p  -u  i 

u 

where 

T 

1 i is  the  unit  vector  between  user  and 
pu  i^h  satellite 

n.  is  the  pseudo-range  measurement  error, 

1 including  the  i^h  satellite  ephemeris 
and  clock  errors 


22 


lilittliiifa-Ai 


WM  M - 1 I I I ljum 


then,  the  estimated  pseudo-range  measurement  error,  Sz^, 
for  this  measurement  is  given  by 

6z.  = lTi  6R  + 54>  + p.  (2.2-13) 

i -P  -u 
u 

The  difference  between  these  two  quantities  (i.e.,  (Sz.^-62.^)  is 
used  to  update  the  filter  measurements.  Note  that  the  pseudo- 
states in  Eq.  (2.2-13)  are  basically  accounting  for  the  satel- 
lite ephemeris  and  clock  errors  in  the  measurement  error  of 
Eq.  (2.2-12)  (assumed  to  be  included  in  the  measurement  error, 
ni).  Thus,  this  is  a better  method  of  accounting  for  the 
unmodeled  satellite  ephemeris  and  clock  errors  than  to  assume 
uncertainty  in  user  position  through  use  of  process  noise. 
However,  this  method  increases  the  filter  complexity  and  com- 
putational load.  Also,  whenever  the  satellite  ephemeris  and 
clock  data  is  uploaded  by  the  ground  tracking  system,  there  is 
an  abrupt  change  in  the  resulting  errors  in  the  ephemeris  and 
clock  data  in  the  measurements  6zit  This  uploading  could  be 
modeled  in  the  user  filter  by  resetting  the  pseudo-states  to 
zero  and  the  pseudo-state  variances  and  covariances  to  their 
initial  (a  priori)  values  after  the  upload.  However,  this 
might  unacceptably  complicate  the  software  implementation,  al- 
though this  aspect  of  the  pseudo-state  argumentation  technique 
has  not  been  investigated. 

The  pseudo-states  were  initialized  with  the  values 
summarized  in  Table  2.2-4.  These  values  were  chosen  on 
the  basis  of  the  magnitudes  of  the  ephemeris  and  clock  er- 
rors illustrated  in  Appendix  B.  Comparison  of  user  per- 
formance projections  using  the  two  methods  of  filter  aug- 
mentation discussed  above  are  presented  in  Section  2.2.3. 


2.2.2  GDOP  and  Satellite  Visibilit 


uation.  GDOP  depends  upon  the  satellite-user  geometry  and 
thus  changes  as  the  satellites  and  user  move  relative  to  one 
another.  The  lower  the  GDOP,  the  lower  and  sensitivity  of 
the  pseudo-range  measurements  to  measurement  noise.  Refer- 
ence 1 contains  a discussion  of  GDOP  as  used  in  this  study. 


The  CEP*  GDOP  reflects  the  horizontal  position  fix 
sensitivity  to  the  pseudo-range  measurement  errors.  This  is 
illustrated  in  Fig.  2.2-3  for  a user  located  at  50  deg  N, 


MUM 


♦Circular  Error  Probable  (CEP)  resulting  from  uncorrelated 
unity-variance  pseudo-range  measurement  errors  from  each 
satellite  (using  all  visible  satellites). 
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100  deg  W using  the  five  satellite  subset  of  the  24-satellite 
Phase  III  configuration  visible  to  the  user,  as  discussed  in 
Appendix  B.  Both  the  conventional  and  wide  beam  directional 
antenna  cases  are  shown.  Note  that  GDOP  (as  used  in  this  study) 
is  only  defined  when  four  or  more  satellites  are  visible. 

A user  with  the  conventional  antenna  (170  deg)  can 
"see"  all  five  satellites,  so;  that  he  can  select  four  which 
give  the  best  GDOP  as  illustrated  (Satellite  Nos.  1,  2,  3 
and  5).  However,  the  wide  beam  directional  antenna  (120 
deg)  user  only  "sees"  four  satellites  which,  in  this  in- 
stance, provide  relatively  poor  GDOP.  However,  as  illustra- 
ted in  the  next  section,  a user  employing  a Kalman  filter 
may  not  be  seriously  penalized  by  this  poor  GDOP,  depending 
on  the  time  available  for  obtaining  the  GPS  position  fix. 

2.2.3  Performance  Projections 

Typical  navigation  performance  projections  of  a 
GPS  manpack  user  with  a directional  antenna  are  summarized 
in  this  section.  These  results  are  presented  in  \he  form 
of  time  histories  of  the  navigation  filter  estimation  error 
as  well  as  the  navigation  filter  standard  deviation  bound- 
ary for  an  error  state.  The  standard  deviation  boundaries 
are  computed  by  taking  the  square  root  of  the  diagonal  ele- 
ment of  the  user  navigation  filter  covariance  matrix  cor- 
responding to  the  particular  error  of  interest  and  plotting 
its  plus  and  minus  values.  According  to  estimation  theory 
(Ref.  5),  if  the  real  world  physical  processes  are  modeled 
properly  in  the  filter,  the  estimation  errors  should  not  go 
outside  the  standard  deviation  boundaries  more  than  31.7%  of 
the  time.  If  this  condition  is  not  satisf ied , then  the  filter 
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does  not  properly  model  the  physcial  process,  and  a filter 
redesign  is  in  order. 

The  time  required  to  acquire  the  satellite  signals 
and  to  obtain  the  satellite  ephemerides  and  other  data  from 
these  signals  has  not  been  accounted  for  in  these  results. 
It  is  assumed  that  the  signals  have  been  acquired  and  data 
removed  before  the  navigation  filter  solution  is  initiated 
at  t = 0.  Acquisition  performance  was  analyzed  previously 
(Ref.  1). 


This  discussion  of  user  performance  is  initiated 
with  a comparison  of  the  two  methods  of  filter  augmentation 
described  in  Section  2.2.1.  User  performance  projections  for 
the  scenario  and  models  of  Section  2.2.1  are  presented  next 
for  three  types  of  antennas:  conventional  (the  baseline), 

wide  beam  directional  and  narrow  beam  directional. 

Filter  Augmentation  — To  investigate  the  two  methods 
of  augmenting  the  navigation  filter  to  account  for  satellite 
ephemeris  and  clock  errors  as  discussed  in  Section  2.2.1, 
the  conventional  antenna  stationary  user  (baseline)  scen- 
ario was  used.  A comparison  of  the  root  sum  squared  (rss) 
horizontal  position  error  for  these  two  methods  is  shown  in 
Fig.  2.2-4.  In  order  to  clearly  show  the  steady-state  behavior 
of  the  filter,  the  large  initial  transient  error  which  oc- 
curs before  four  measurements  have  been  taken  is  not  shown. 

The  addition  of  process  noise  in  the  filter  gain 
computation  causes  the  sawtooth  pattern  of  the  standard  de- 
viation boundary  (illustrated  in  the  top  graph)  and  unrealis- 
tically limits  the  reduction  of  the  uncertainty  of  the  position 
estimates  as  additional  measurements  are  taken. 
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Figure  2.2-4  Filter  Augmentation  Comparison 

The  performance  of  the  pseudo-state  augmented  filter 
(bottom  graph)  is  more  realistic  for  the  stationary  user,  since 
the  uncertainty  in  the  estimates  decreases  with  additional 
measurements . These  results  indicate  that  the  pseudo-state 
augmentation  is  the  better  method  of  accounting  for  satellite 
ephemeris  and  clock  data  error,  although  the  relatively  small 
improvement  shown  in  Fig.  2.2-4  may  not  justify  the  additional 
filter  complexity.  Nonetheless,  subsequent  studies  included 
in  this  section  incorporate  pseudo-states. 

Baseline  — In  order  to  provide  a standard  with  which 
to  compare  performance  of  users  with  directional  antennas, 
performance  projections  for  a stationary  user  with  a conven- 
tional (hemispherical)  antenna  was  generated,  using  the  scen- 
ario and  models  (e.g.,  receiver,  clock,  filter)  given  in  Sec- 
tion 2.2.1.  Measurements  were  made  sequentially  to  Satellite 
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Nos.  1,  2,  3 and  5,  at  a 10  sec/satellite  rate.  As  illustra- 
ted in  Fig.  2.2-3,  the  CEP  GDOP  for  this  user  varies  between 
1 and  3. 

Performance  for  this  baseline  case  is  illustrated  in 
Fig.  2.2-5.  Again,  the  initial  transient  errors  are  not  shown 
in  order  to  clarify  the  steady-state  performance  of  the  filter. 
Further  design  of  the  filter  and,  in  particular,  the  selection 
of  the  initial  pseudo-state  uncertainty,  may  improve  this  base- 
line performance.  However,  detailed  filter  design  was  not  the 
objective  of  this  study. 

Wide  Beam  User  — The  directional  antenna  for  this 
user  has  a beamwidth  of  120  deg  and  is  pointed  vertically. 

Four  satellites  are  visible  to  this  user  throughout  the  10- 
minute  time  duration  of  this  case.  However,  this  user  experi- 
ences poor  GDOP,  as  shown  in  Fig.  2.2-3,  with  a corresponding 


TIME  (AFTER  DATA  ACQUISITION)  <MC|  TIME  (AFTER  DATA  ACQUISITION)  (mc) 

Figure  2.2-5  Baseline  Performance 


degradation  in  user  performance  as  illustrated  in  Fig.  2.2-6. 
The  uncertainty  in  the  estimates,  as  indicated  by  the  standard 
deviation  boundary,  is  strongly  dependent  upon  GDOP.  This  is 
particularly  evident  in  the  horizontal  position  error  results. 
The  horizontal  position  error  standard  deviation  after  four 
sequential  measurements  and  at  the  end  of  the  simulation  are 
summarized  in  Table  2.2-5.  After  four  measurements,  this 
standard  deviation  for  the  wide  beam  user  is  a factor  of  8.7 
larger  than  that  of  the  baseline  user,  because  of  the  poor 
GDOP.  However,  as  additional  measurements  are  taken  this 
error  standard  deviation  is  reduced  to  a value  comparable  to 
that  for  the  baseline  user. 


x - ESTIMATION  ERROR 
0 - STANDARD  DEVIATION  BOUNDARY 


Figure  2.2-6  Wide  Beam  User  Performance 
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TABLE  2.2-5 


PERFORMANCE  SUMMARY 


USER 

HORIZONTAL 

STANDARD 

POSITION  ERROR 
DEVIATION 

40  sec 

600  sec 

Baseline 

15  m 

9.4  m 

Wide  Beam 

131  m 

9.4m  i 

Narrow  Beam  User  — An  Army  manpack  GPS  user  may 
be  subject  to  an  environment  in  which  foliage  attenuates 
the  GPS  navigation  signal  and/or  interference  from  elec- 
tronic jamming  is  present,  preventing  acquisition  and  use 
of  the  GPS  navigation  signals.  One  method  of  overcoming 
these  problems  is  through  the  use  of  a narrow  beam  direc- 
tional antenna  which  can  be  physically  pointed  by  the  user 
in  the  direction  of  a GPS  satellite.  The  navigation  signal 
strength  would  be  enhanced  by  the  gain  of  the  antenna,  while 
reducing  the  effects  of  interference  and  jamming  by  virtue 
of  the  spatial  properties  of  the  antenna. 


In  order  to  obtain  a navigation  solution,  pseudo- 
range measurements  to  at  least  four  satellites  are  required 
(assuming  no  a priori  knowledge  of  position  or  time).  Thus, 
the  antenna  must  be  separately  pointed  physically  to  at 
least  four  satellites  in  sequence.  The  user  computer  would 
indicate  approximately  where  to  point  the  antenna,  and  the 
receiver  could  indicate  when  the  desired  signal  has  been 
acquired. 


In  this  study  of  a user  with  a narrow  beam  directional 
antenna,  the  pointing  of  the  antenna  to  each  satellite  is  simu- 
lated by  taking  five  measurements  (at  a 10  sec  rate)  to  one 
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satellite  and  switching  to  the  next  satellite  for  five  measure- 
ments, and  so  forth.  The  satellite  selection  sequence  was  1, 

2,  3,  4,  and  5.  Holding  on  individual  satellites  for  five 
measurements  is  reasonable  in  view  of  the  fact  that  the  antenna 
must  be  reoriented  by  the  user  to  acquire  a new  satellite.  With 
this  technique,  all  satellites  above  the  horizon  are  available 
to  the  user.  Time  required  to  reposition  the  antenna  to  a new 
antenna  has  not  been  taken  into  account  in  this  simulation. 
During  this  repositioning  period,  the  filter  could  be  in  a 
"hold"  condition  in  which  no  measurements  or  processing  is 
performed.  In  addition,  the  time  required  to  acquire  the 
satellite  signal  and  data  has  not  been  included.  The  time 
required  to  perform  both  antenna  repositioning  and  satellite 
signal  and  data  acquisition  would  simply  add  to  the  time  re- 
quired to  obtain  a navigation  fix  as  illustrated  in  the  results 
presented  below. 


Typical  performance  for  this  user  is  illustrated  in 
Fig.  2.2-7.  The  horizontal  position  error  has  been  plotted 
with  a vertical  logarithm  scale  in  order  to  illustrate  the 
wide  variation  in  dynamic  range  of  this  error  as  the  measure- 
ments are  processed.  This  curve  indicates  that  measurements 
to  the  fourth  satellite  reduce  the  horizontal  position  error 
to  less  than  10  m.  After  the  measurements  to  the  fourth  satel 
lite  the  horizontal  position  error  falls  outside  the  standard 
deviation  boundary.  This  is  indicative  of  a minor  filter  de- 
sign problem  which  was  not  pursued  in  this  study.  The  other 
user  errors  (altitude,  clock  phase  and  clock  frequency)  shown 
in  Fig.  2.2-7  are  also  strongly  dependent  on  acquisition  of 
the  fourth  satellite.  Additional  measurements  to  the  same  set 
of  five  satellites  does  not  significantly  reduce  the  error 
after  the  initial  set  of  measurements  has  been  processed. 


Performance  for  the  narrow  beam  antenna  user  compares 
favorably  with  that  of  the  baseline  user.  The  time  required 
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SATELLITE  No,  4 


TIME  (AFTER  DATA  ACQUISITION)  (mc) 


TIME  (AFTER  DATA  ACQUISITION)  fuel 


TIME  (AFTER  DATA  ACQUISITION)  (we) 


TIME  (AFTER  DATA  ACQUISITION)  (we) 


Figure  2.2-7  Narrow  Beam  User  Performance 


to  obtain  a satisfactory  navigation  solution  in  this  scenario 
is  about  3 min  plus  the  time  required  to  reposition  the  antenna 
three  times  and  acquire  satellite  signals  and  data. 


2.2.4  Performance  Summar 


The  performance  of  a GPS  user  whose  satellite  visi- 
bility is  restricted  with  a directional  antenna  is  not  adversely 
affected  by  this  restriction  as  long  as  the  measurements  to  a 
sufficient  number  of  satellites  are  processed  with  a Kalman 
navigation  filter.  Augmentation  of  this  filter  with  pseudo- 
states may  be  advantageous  for  a stationary  user.  Large  GDOPs 
can  result  in  a significant  degradation  in  performance  although 
movement  of  the  satellites  alleviates  this  problem  within  10  min, 
The  use  of  a very  narrow  beamwidth  antenna  that  is  physically 
pointed  to  the  satellites  can  provide  accurate  position  esti- 
mates, although  the  time  required  to  obtain  these  estimates  is 


■y.-rr- 


also  affected  by  the  time  required  to  position  the  antenna. 

In  general,  position  estimation  errors  of  less  than  10  m are 
possible  for  users  with  directional  antennas.  The  accuracy  of 
the  estimates  generated  here  could  be  improved  with  further 
filter  design  refinements. 


2.3  ANTENNA  ARRAY  SYSTEMS 

Directional  antennas  can  be  constructed  from  an  ar- 
ray of  discrete  isotropic  antennas  or  elements.  Depending 
upon  the  relative  position  of  the  elements  with  respect  to 
each  other  and  the  signal  frequency,  the  antenna  array  gain 
will  be  greater  in  some  directions  than  in  others.  If  the 
output  signal  from  each  element  is  delayed  with  respect  to 
a reference  element  by  a given  amount,  the  direction  of  this 
maximum  array  gain  can  be  changed.  This  results  in  a direc- 
tional antenna  with  a main  beam  that  can  be  pointed  in  any 
direction  simply  by  changing  the  delays.  Furthermore,  it 
is  possible  to  reduce  the  effects  of  directional  interfer- 
ence by  placing  a null  (direction  in  which  the  antenna 
gain  is  low)  of  the  antenna  beam  pattern  in  the  direction  of 
the-  interference. 

Antenna  arrays  can  provide  both  spatial  and  fre- 
quency filtering  to  enhance  the  desired  signal  while  reduc- 
ing undesired  directional  interference  or  jamming.  This 
section  presents  a discussion  of  antenna  arrays  and  their 
applicability  and  feasibility  for  use  by  a GPS  user,  par- 
ticularly the  Army  manpack  user.  A basic  discussion  of 
antenna  array  processing  is  presented  in  Section  2.3.1. 

The  reduction  of  directional  interference  is  considered 
in  Section  2.3.2  with  an  outline  of  the  optimum  array  pro- 
cessor. Finally,  in  Section  2.3.3  the  application  of  array 
antennas  by  GPS  users  is  discussed. 


2.3.1  Antenna  Array  Processing 


Antenna  arrays  consist  of  a set  of  antennas  or 
elements  whose  outputs  are  delayed  relative  to  some  refer- 
ence element  signal  phase  such  that  the  signals  arriving 
from  a given  direction  will  be  in  phase  after  the  delaying 
operation.  If  these  delayed  outputs  are  summed,  the  signals 
that  are  in  phase  (i.e.,  those  from  the  given  direction)  will 
be  enhanced  or  increased  in  strength  by  the  number  of  array 
elements.  Signals  arriving  from  other  than  the  given  direc- 
tion will  not  be  in  phase  after  this  delaying  operation, 
and  thus  will  be  less  enhanced  (or  degraded),,  depending  upon 
their  phase  (i.e.,  direction).  The  delays  can  be  changed  to 
"steer"  the  main  beam  of  the  array  in  any  desired  direction, 
thus  providing  spatial  discrimination  or  filtering. 

A simple  example  will  help  to  illustrate  the  con- 
cept. Consider  the  linear  array  with  four  elements  in  Fig. 
2.3-1.  Assume  a signal  is  arriving  with  a planar  wavefront 
as  shown.  If  the  outputs  of  the  elements  are  delayed  by 


Figure  2.3-1  Four  Element  Array 
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x . = — sin  <p  (2 . 3-1 ) 

1 c 

where 

d . is  the  spacing  from  reference  element  1 
1 to  element  i 

c is  the  velocity  ot  propagation 

the  signal  at  each  delayed  output  will  then  be  in  phase  with 
the  signal  at  element  1.  If  these  delayed  signals  are  then 
summed,  the  signal  strength  will  be  four  times  as  large  as 
for  a single  element.  Any  signal  arriving  at  a direction 
other  than  <J>  will  not  be  in  phase  at  the  delayed  outputs  and 
thus  will  not  be  enhanced.  This  array  can  be  steered  in  any 
direction  (<jj)  simply  by  adjusting  the  delays  as  given  by 
Eq.  (2.3-1). 

The  power  or  beam  pattern  (Appendix  A)  for  this  four 
element  linear  array  (assuming  each  element  is  an  ideal  isotro- 
pic antenna  with  element  spacing  d = A/2 , where  A = signal  wave- 
length) is  illustrated  in  Fig.  2.3-2  for  a steered  direction  of 
<|>  = 30  deg.  This  figure  is  a polar  plot  of  the  array  gain  as 
a function  of  angle  relative  to  the  center  of  the  array  in 
the  plane  of  the  array.  The  shape  of  the  beam  pattern  (and 


Figure  2.3-2  Typical  Beam  Pattern  of  Four 

Element  Linear  Array 
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the  width  of  the  main  beam)  will  vary  with  steering  direc- 
tion. Additional  elements  would  reduce  the  width  of  the 
main  beam  while  increasing  its  gain.  Because  of  array  sym- 
metry, the  beam  pattern  revolves  around  the  array  axis  . Several 
secondary  beams  or  sidelobes  are  also  evident  for  this  lin- 
ear array.  The  relative  amplitude  of  these  sidelobes  as 
well  as  the  width  of  the  main  beam  can  be  controlled  by 
"shading,"  that  is,  by  controlling  the  amplitude  of  the  sig- 
nal from  each  element  (Ref.  15). 

Array  antenna  beam  patterns  also  exhibit  nulls,  as  evi- 
dent in  Fig.  2.3-2.  If  an  interference  signal  arrived  from 
the  direction  of  one  of  these  nulls,  this  interference  would 
be  reduced  or  eliminated  at  the  output  of  the  array.  With 
additional  phase  shifting  and  further  signal  processing,  a 
null  can  be  placed  in  any  desired  direction,  as  discussed 
in  the  next  section. 

2.3.2  Interference  Reduction  With  Arrays 

The  design  of  antenna  arrays  to  provide  interfer- 
ence reduction  or  elimination  has  been  investigated  by 
several  authors  (Refs.  16  to  21).  The  optimum  array/pro- 
cessor configuration  for  the  reduction  of  a directional 
interference  has  been  shown  to  utilize  two  beams,  one 
steered  in  the  desired  signal  direction  (main  beam)  and 
the  other  steered  in  the  direction  of  the  interference 
(interference  beam)  (Refs.  17  and  20).  The  output  of  the 
interference  beam  is  passed  through  a weighting  filter 
and  then  subtracted  from  the  signal  beam  output,  as  illus- 
trated in  Fig.  2.3-3.  The  beamformers  consist  of  delay 
elements  and  frequency  filters.  The  weighting  filter 
modifies  the  output  of  the  interference  beam  so  that  it 
just  equals  the  contribution  of  the  interference  in  the 
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Figure  2.3-3  Interference  Reduction  Processor 

main  beam  output.  Thus  the  output  of  the  weighting  filter 
will  just  cancel  the  interference  signal  in  the  main  beam. 
This  weighting  filter  is  a function  of  the  directior  to  which  the 
main  and  interference  beams  are  steered  as  well  as  the  inter- 
ference power  ( I )-to-isotropic  noise  power  (N)  ratio,  (i.e., 
I/N).  The  isotropic  noise  power  (N)  is  the  power  of  iso- 
tropic noise  at  each  element  of  the  array,  assuming  the  noise 
is  uncorrelated  between  elements. 


The  beam  pattern  of  a four  element  array  with  this 
interference  reduction  processor  is  illustrated  in  Fig.  2.3-4 
for  an  array  steered  to  a desired  signal  at  30  deg.  and  an 
interference  at  60  deg  with  I/N  = 100  dB.  This  pattern  is 
symmetrical  about  the  array  axis,  as  in  Fig.  2.3-2;  how- 
ever, the  antenna  pattern  is  significantly  different.  This 
beam  pattern  indicates  that  the  processor  basically  steers 
a beam  pattern  null  in  the  direction  of  the  interference,  thus 
reducing  or  eliminating  its  effect  from  the  antenna  output. 
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Figure  2.3-4  Beam  Pattern  With  Interference  Reduction 

If  the  interference  direction  is  close  to  the  direc- 
tion of  the  desired  signal,  the  array  gair  in  the  desired 
signal  direction  must  be  compromised.  However,  if  the  in- 
terference is  outside  the  main  beam  (that  is,  at  or  beyond 
the  main  beam  nulls),  the  antenna  array  gain  or  directiv- 
ity of  the  main  beam  (in  the  direction  of  the  desired  sig- 
nal) is  not  significantly  affected  (Ref.  20).  The  influ- 
ence of  the  interference  direction  on  the  array  gain  for  a 
four  element  array  is  illustrated  in  Fig.  2.3-5.  This  fig- 
ure shows  the  gain  of  the  array  in  the  direction  of  the  de- 
sired signal  at  30  deg,  as  the  interference  direction  is 
changed.  The  array  gain  is  not  significantly  affected  un- 
less the  interference  is  within  30  deg  of  the  desired  sig- 
nal direction.  The  shape  of  this  gain  function  is  dependent 
upon  the  desired  signal  direction,  I/N  ratio,  and  the  num- 
ber and  spacing  of  array  elements.  The  interference  direction 
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Figure  2.3-5  Array  Gain  With  Interference  Reduction 

could  be  closer  to  the  desired  signal  direction  without  af- 
fecting array  gain  if  the  array  contained  more  elements  or 
the  array  was  steered  closer  to  broadside  (<J>  = 0 deg.). 

The  primary  interest  in  employing  interference 
reduction  with  antenna  arrays  is  to  improve  the  tolerance 
of  the  receiving  system  to  directional  interference  or  jam- 
ming. The  theoretical  reduction  in  interference  provided  by 
the  optimum  interference  reduction  processor  is  illustrated 
in  Fig.  2.3-6.  This  shows  the  reduction  in  the  strength 
of  the  interference  at  the  antenna  array  processor  output 
as  a function  of  the  interference  direction  with  the  array 
steered  to  30  deg  for  several  values  of  interf erence-to- 
noise  (I/N)  ratio.  This  interierence  reduction  is  the 
ratio  of  the  gain  of  the  array  (with  the  processor)  in  the 
direction  of  interference  to  the  gain  in  the  steered 
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Figure  2.3-6  Theoretical  Interference  Reduction 


direction.  The  reduction  in  interference  is  dependent  upon 
the  strength  (power)  of  the  interference;  that  is,  the 
stronger  the  interference,  the  greater  the  reduction.  How- 
ever, if  the  interference  is  close  to  the  direction  of  the 
desired  signal,  then,  as  shown  in  Fig.  2.3-5,  the  array 
gain  is  decreased  and  the  interference  reduction  technique 
becomes  ineffective,  as  indicated  in  Fig.  2.3-6  in  the 
vicinity  of  the  30  deg  value. 

The  theoretical  reduction  in  interference  obtain- 
able with  the  interference  reduction  technique  adds  to  the 
interference  margin  (i.e.  , the  jamming-to-signal  power 
ratio  ( I/S ))  of  the  receiver.  Thus,  if  the  receiver  noise 
level  (N)  is  -201  dBW-Hz  (Ref.  1)  and  the  jammer  power  (I) 
at  the  antenna  is  -110  dBW  (corresponding  to  a 1 kW  ground 


:S 

| 

t 


i 


a 


1 


jammer  lo  km  from  user, 
to-noise  ratio  (I/N)  is 
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-110  - (201)  dB  (2.3-2) 

= 91  dB 

Prom  Fig.  2.3-6,  this  value  of  I/N  corresponds  to  a reduction 
interference  at  the  output  of  the  processor,  of  199  dB  for  an 
interference  at  60  deg  and  signal  at  30  deg.  The  GPS  re- 
ceiver has  an  additional  Interference  margin  of  35  to  49  dB 
(Ref.  28).  ihus,  the  total  theoretical  interference  margin 
of  the  interference  reduction  processor-receiver  system 
would  be  234  to  248  dB. 


Caution  should  be  exercised  in  using  these  perform- 
ance indices  as  they  are  based  on  an  optimum  interference  re- 
duction processor  which  may  not  be  possible  to  implement.  The 
finite  nature  of  any  implementation  scheme  (e.g.,  computer 
word  size),  non-ideal  antenna  beam  patterns  and  radio  fre- 
quency circuitry  implementation  problems  may  make  it  diffi- 
cult to  obtain  these  large  reductions  in  interference. 


he  optimum  processor  assumes  that  the  directions 
Of  the  desired  signal  and  interference  as  well  as  the 
interference-to-noise  ratio  are  known.  It  may  be  unrealis- 
tic to  estimate  either  the  interference  direction  or  strength, 
owever,  the  direction  of  the  desired  signal  is  usually  known 
Several  authors  (Refs.  16  and- 21)  have  proposed  adaptive  an-  ' 
tenna  array  systems  that  are  similar  to  the  optimum  interfer- 
ence reduction  processor  and  do  not  require  knowledge  of  in- 
terference direction  or  strength,  but  do  require  knowledge 
Of  signal  structure  or  direction.  Such  a system  has  been 
ilt  and  tested  using  a four  element  linear  array  (Ref,  20). 


The  experimental  performance  of  this  system  was  close  to 
theoretical  optimum.  For  I/N  = 21.5  dB,  the  reduction  in 
i nterf erence-to-signal  ratio  was  42  dB.  Thus,  significant 
reductions  in  interference  were  obtained  with  modest  I/N 
rat ios . 


Another  possible  approach  to  implementing  the  opti- 
mum processor  is  to  use  an  adaptive  tracking  technique  (e.g. , 
a phase  lock  loop)  to  keep  the  interference  beam  directed 
toward  the  interference,  while  the  GPS  navigation  computer 
directs  the  main  beam  in  the  direction  of  the  desired  signal. 
This  technique  may  be  easier  to  implement  than  the  adaptive 
antenna  array  system  of  the  last  paragraph,  although  it  has 
not  been  investigated  here. 

2.3.3  Applicability  to  GPS  Manpack  User 

The  GPS  manpack  user  operating  in  a battlefield 
environment  is  likely  to  encounter  directional  interference 
such  as  jamming  by  the  adversary.  The  interference  reduc- 
tion technique  using  an  array  of  antennas  as  discussed  in 
the  preceding  section  is  capable  of  providing  a significant 
reduction  in  a user's  vulnerability  to  jamming.  The  four 
element  array  of  the  last  section  is  less  than  30  cm  long 
at  the  GPS  navigation  signal  frequencies  and  thus  physically 
feasible.  However,  the  implementation  of  this  array  may 
require  a considerable  amount  of  signal  processing  equip- 
ment in  order  to  obtain  the  desired  performance.  This 
equipment  could  be  simplified  if  a two  element  array  was 
employed,  at  the  expense  of  a reduction  in  performance. 

The  beam  pattern  of  a two  element  array  employing  inter- 
ference reduction  is  illustrated  in  Fig.  2.3-7.  This  fig- 
ure indicates  that  the  two  element  array  has  fewer  nulls 
and  thus  may  be  able  to  receive  the  navigation  signals  from 
several  satellites  at  the  same  time.  However,  the  gain  of 
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Figure  2.3-7  Two  Element  Beam  Pattern  With 

Interference  Reduction 

the  array  varies  considerably  with  the  direction  of  inter- 
ference as  illustrated  in  Fig.  2.3-8.  Because  of  the  low 
margin  in  the  GPS  navigation  signal  power  budget,  this  var- 
iation in  antenna  gain  may  preclude  the  use  of  a two  ele- 
ment array.  The  reduction  in  interference  which  can  be 
provided  by  the  two  element  array  is  illustrated  in  Fig. 
2.3-9.  This  figure  indicates  that  significant  reduction 
can  be  achieved  as  long  as  the  difference  in  directions  of 
the  desired  signal  and  interference  is  greater  than  30  deg. 

Two  and  four  element  "steerable  null  antenna  pro- 
cessors" for  Army  field  communication  systems  are  presently 
being  investigated  by  U.S.  Army  Electronics  Command  (Refs. 
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Although  these  processors  are  designed  for  lower  frequencies 
than  that  of  the  GPS  navigation  signal,  their  design  could 
be  modified  to  operate  at  L-band.  Thus  modified,  these  pro- 
cessors appear  to  be  feasible  for  GPS  manpack  users,  al- 
though details  of  their  performance  remain  to  be  investi- 
gated. 


It  thus  appears  that  interference  reduction  using 
antenna  arrays  are  feasible  for  GPS  manpack  users.  Addi- 
tional studies  of  possible  methods  of  implementing  the 
optimum  interference  processor  are  required  in  order  to 
determine  the  feasibility  of  obtaining  performance  that 
is  suitable  for  the  GPS  user  as  well  as  to  estimate  the 
costs  of  such  a processor. 


2.4  ANTENNA  DEPLOYMENT  PROCEDURES 

The  directional  antennas  considered  in  this  re- 
port require  proper  orientation  and  deployment  in  order  that 
the  user  computer  can  select  those  satellites  that  are  within 
the  main  beam  of  the  antenna.  This  section  discusses  these 
deployment  procedures  for  three  types  of  directional  anten- 
nas considered  in  this  report.  These  are  the  wide  beam  and 
narrow  beam  directional  antennas  and  the  array  antenna. 

2.4.1  Wide  Beam  Antenna  Deployment 

The  wide  beam  antenna  used  in  the  preceding  per- 
formance studies  is  assumed  to  be  a vertically  pointed  para- 
bolic type  with  a beamwidth  of  120  deg.  This  antenna  should 
be  mounted  on  an  adjustable  tripod  with  a two-axis  leveling 
device,  such  as  plumb  bob  or  bubble  level,  which  is  capable 
of  positioning  the  antenna  to  point  vertically.  The  receiver 
is  then  activated. 
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2.4.2  Narrow  Beam  Antenna  Deployment 


The  narrow  beam  antenna  is  also  assumed  to  be  a 
parabolic  type  but  with  a very  narrow  beamwidth.  This  an- 
tenna is  manually  positioned  by  the  user  to  sequentially 
point  to  individual  satellites.  This  will  again  require 
mounting  the  antenna  on  a tripod  with  a two-axis  leveling 
device.  In  addition,  a measure  of  both  elevation  and  azi- 
muth angles  must  be  provided. 

The  deployment  procedure  for  the  narrow  beam  an- 
tenna again  consists  of  unpacking  and  assembling  the  antenna 
and  tripod  and  connecting  to  the  receiver  with  the  cable. 

The  tripod  and  antenna  are  then  leveled  and  the  antenna 
azimuth  aligned  with  a reference  direction  (e.g.,  magnetic 
north)  using  a compass.  The  user  GPS  navigation  computer 
is  then  turned  on  and  queried  as  to  the  vertical  and  azi- 
muth angles  at  which  the  antenna  should  be  positioned  to 
receive  the  desired  satellite  navigation  signal.  The 
antenna  is  then  set  to  these  angles  and  the  receiver  turned 
on.  If  the  navigation  signal  is  received,  the  user  can  be 
so  notified  and  the  measurement  taking  commenced.  If  the 
signal  is  not  received,  the  measurements  can  be  skipped, 
and  an  attempt  made  to  locate  another  satellite.  After  a 
sufficient  set  of  measurements  has  been  taken,  the  antenna 
can  be  repositioned  to  the  next  satellite  as  directed  by 
the  computer.  The  exact  procedure  for  operating  the  GPS 
navigation  computer/receiver  and  antenna  will  depend  upon 
their  implementation;  however,  the  procedure  outlined  above 
is  typical  of  what  must  be  done  to  use  a narrow  beam  antenna. 

2.4.3  Array  Antenna  Deployment 

The  deployment  procedures  for  the  array  antenna 
will  depend  upon  the  number  of  elements  and  the  physical 
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construction  of  the  array.  Assuming  the  array  elements  are 
arranged  in  either  a line  or  planar  geometry,  the  deploy- 
ment procedure  will  be  similar  to  that  for  the  wide  beam 
antenna.  However,  the  axis  of  the  array  must  be  aligned 
with  a reference  direction  (e.g.,  magnetic  north). 

The  deployment  procedure  for  the  array  antenna  con- 
sists of  unpacking  and  assembling  the  antenna  and  tripod 
and  connecting  the  antenna  to  the  receiver.  The  antenna  is 
then  leveled  and  the  axis  aligned  with  the  reference  direc- 
tion. The  receiver  is  then  activated  and  measurements  pro- 
cessed. Some  implementations  of  the  interference  reduction 
processor  may  require  operator  assistance  in  identifying 
probable  directions  of  jammer  signals,  in  which  case  the 
operator  would  input  the  direction  of  the  battlefront  or 
his  best  estimate  of  the  jammer  direction.  Sequencing  from 
satellite-to-satellite  could  be  accomplished  automatically, 
under  computer  control. 


3. 


USER  VULNERABILITY  AND  THREAT  COUNTERING 


The  ground-based  users  of  GPS,  such  as  the  Army  man 
pack  user,  are  likely  to  be  operating  in  a battlefield  envi 
ronment  where  electronic  jamming  of  the  GPS  navigation  sig- 
nal is  present.  The  degree  of  vulnerability  of  the  GPS 
user  will  depend  upon  his  location  relative  to  the  jammer, 
the  power  of  the  jammer,  the  jamming  signal  structure,  etc. 
The  user  may  be  able  to  employ  threat  countermeasures  to 
reduce  his  vulnerability  to  jamming  and  other  directional 
interferences . 

This  chapter  considers  the  vulnerability  of  ground 
based  users  of  GPS  and  some  possible  threat  countering  tech 
niques.  The  scenario  and  threat  models  used  to  evaluate 
user  vulnerability  are  presented  in  Section  3.1.  Included 
in  this  section  are  a discussion  of  the  effectiveness  of 
various  jamming  signals  and  the  jamming  signal  degradation 
due  to  the  propagation  path  between  jammer  and  user.  Sev- 
eral techniques  for  countering  the  jammer  threat  are  pre- 
sented in  Section  3.2.  Finally,  an  evaluation  of  the 
vulnerability  of  a GPS  manpack  user  to  electronic  jamming 
for  the  threat  models  and  scenarios  of  Section  3.1  is  pre- 
sented in  Section  3.3.  The  effectiveness  of  the  threat 
countering  techniques  of  Section  3.2  is  also  discussed. 


3.1  SCENARIOS  AND  MODELS 

The  scenarios  and  models  used  in  the  evaluation 
of  the  GPS  manpack  user  vulnerability  to  electronic  jamming 
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are  discussed  in  this  section.  These  scenarios  and  models 
were  designed  to  provide  a realistic  assessment  of  user 
vulnerability  to  jamming. 

3.1.1  Scenarios 

This  study  is  concerned  with  an  Army  manpack  GPS 
user  in  a battlefield  environment  where  jamming  is  present. 
The  user  is  assumed  to  be  stationary  with  his  antenna  lo- 
cated 2 m above  ground  level.  This  antenna  is  initially 
assumed  to  be  the  conventional  antenna  defined  in  Chapter  2; 
however,  directional  antennas  will  also  be  considered.  The 
receiver  is  a Y-model  (Ref.  24)  capable  of  operating  with  a 
jammer  power-to-navigation  signal  power  ratio  (J/S)  of 
34  dB  to  49  dB,  depending  upon  operating  mode  (i.e.,  ini- 
tial signal  acquisition,  carrier  tracking,  etc.),  and  a 
receiver  input  signal  level  (P-code)  greater  than  -163  dBW. 

The  interference  rejection  capability  of  a receiver 
or  its  threshold  jammer-to-signal  power  ratio  is  heavily 
dependent  upon  the  receiver  design.  Since  this  study  was 
not  concerned  with  receiver  design,  the  values  of  this 
threshold  used  are  the  minimum  values  as  specified  by  the 
Air  Force  (Ref.  24).  Thus,  the  vulnerability  results  pre- 
sented here  are  for  worst  case  threshold  values.  Improved 
values  may  be  possible  with  particular  receiver  designs, 
thus  reducing  the  user  vulnerability  presented  in  this 
report . 


The  electronic  jammer  is  assumed  to  emit  a radio 
frequency  continuous  wave  (CW)  signal  at  the  frequency  of 
the  primary  GPS  navigation  signal  (i.e.,  1575.42  MHz), 
which  is  directed  toward  the  GPS  user.  Two  jammers  will 
be  considered  in  this  report.  One  is  a ground-based  jam- 
mer with  an  antenna  height  of  10  m and  the  other  a standoff 
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airborne  jammer  with  an  altitude  of  3 km.  The  jammer  ef- 
fective isotropic  radiated  power  (EIRP)  will  be  treated  as 
a parameter,  but  within  limits  of  1 W to  100  kW.  Further 
details  of  the  jammer  are  given  in  Section  3.1.2. 

The  signal  path  between  the  jammer  and  the  user  is 
initially  assumed  to  be  a clear  (i.e.,  unobstructed)  path. 
However,  consideration  is  given  to  the  effects  on  user 
vulnerability  of  obstruction  in  the  signal  path,  such  as 
an  RFI  (radio  frequency  interference)  screen.  Further  de- 
tails of  the  signal  path  are  given  in  Section  3.1.3. 

3.1.2  Threat  Models 

There  are  several  signal  structures  that  a jammer 
can  radiate  in  order  to  deny  user  access  to  GPS.  The  effec- 
tiveness of  these  signals  depends  upon  the  user  antenna  and 
receiver  configuration.  The  GPS  navigation  signal  is  a 
continuous  sine  wave  bi-phase  modulated  with  a ranging  and 
data  code  sequence.  The  GPS  receiver  mixes  or  correlates 
this  signal  with  a similar  reference  code  sequence  and,  if 
these  code  sequences  match,  correlation  is  achieved,  thus 
removing  the  phase  modulation  and  producing  the  original  CW 
signal  with  power  equal  to  that  of  the  received  signal  (S). 

The  jammer  could  attempt  to  emulate  the  GPS  naviga- 
tion signal  by  transmitting  a similar  bi-phase  modulated  sig- 
nal. However,  unless  the  jammer  code  sequence  is  identical 
to  the  navigation  signal  code  sequence  that  the  receiver  is 
searching  for  or  tracking,  no  correlation  is  achieved.  In 
fact,  when  this  or  similar  noise-like  jamming  or  interference 
signals  are  mixed  with  the  receiver  reference,  the  range  of 
their  frequency  spectrum  is  spread  or  increased  by  an  amount 
(Ref.  25)  equal  to  the  bandwidth  of  the  reference  code 
sequence  (which  is  equal  to  twice  the  ranging  code 
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sequencing  or  chipping  rate),  and  thus  its  power  is  distrib- 
uted over  a wider  frequency  spectrum.  Since  the  correlator 
is  followed  by  a very  narrow  band  filter,  the  effective  jam- 
mer power  will  be  reduced  by  the  amount  its  spectrum  is 
spread . 


Since  interference  signals  are  spread  across  the  GPS 
band,  an  effective  jamming  signal  is  one  whose  power  is  con- 
centrated in  a very  narrow  frequency  spectrum  within  this  band. 

A continuous  sine  wave  (CW)  signal  satisfies  this  requirement 
and  is  thus  an  effective  jamming  signal  (disregarding  potential 
countering  techniques  such  as  the  notch  filter  described  in 
Section  3.2).  The  jamming  threat  used  in  this  study  is  thus 
assumed  to  be  radiating  a CW  signal  in  all  directions  with  an 
EIRP  of  from  1 W to  100  kW . Other  jamming  signals  (e.g.,  pulse) 
may  be  equally  effective.  The  objective  here  was  to  evaluate  a 
number  of  jammer  deployment  methods  and  countering  techniques  for 
a single,  representative  jammer  signal  structure. 

3/.1.3  Propagation  Path  Models 

To  provide  a realistic  assessment  of  the  vulnera- 
bility of  GPS  users  to  jamming,  the  factors  which  influence 
the  propagation  of  radio  waves  between  the  jammer  and  user 
must  be  considered.  These  factors  limit  the  effective  range 
within  which  the  jammer  can  deny  a user  access  to  GPS  and  are 
particularly  important  to  the  ground-based  users. 

The  earth  is  the  most  important  factor  influencing 
radio  wave  propagation,  especially  for  ground-based  GPS 
users  and  jammers.  It  provides  a surface  to  reflect  the 
waves  and,  because  of  its  curvature,  limits  line-of-sight 
distances  for  direct  propagation  of  waves.  The  principle 
mechanism  for  propagation  of  radio  waves  at  GPS  frequencies 
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is  the  tropospheric  wave,  which  travels  in  the  lower  15  km 
of  the  earth's  atmosphere.  There  are  several  paths  by  which 
this  wave  may  travel  between  a transmitter  and  receiver; 
these  paths  can  be  divided  into  two  regions,  within  line-of- 
sight  and  beyond  line-of-sight . Within  line-of-sight , the 
wave  normally  has  two  components.  One  travels  directly  from 
transmitter  to  receiver,  while  the  other  reaches  the  receiver 
as  a result  of  reflection  from  the  surface  of  the  earth. 
Beyond  line-of-sight,  the  wave  is  diffracted  around  the 
curvature  of  the  earth.  Well  beyond  the  line-of-sight  the 
wave  can  be  reflected  from  scatterers  in  the  upper  tropos- 
phere. In  addition,  within  either  region,  the  wave  is  in- 
fluenced by  the  composition  of  the  atmosphere  and  by  foli- 
age, if  present.  A summary  of  these  propagation  paths  is 
presented  in  Table  3.1-1. 


TABLE  3.1-1 
PROPAGATION  PATHS 


REGION 

PATH 

Within 

• Direct 

Line-of-Sight 

• Reflection 

Beyond 

• Diffraction 

Line-of-Sight 

• Tropospheric  Scatter 

Either  of  above 

• Atmosphere 

• Foliage 

The  amount  of  energy  loss  by  a wave  traveling  through 
these  different  paths  is  a function  of  distance  traveled,  an- 
tenna heights,  wavelength,  composition  of  the  earth's  sur- 
face and  atmosphere,  etc.  This  transmission  or  propagation 
loss  can  be  computed  using  theoretical  models  and  empirical 
data.  Details  of  the  models  and  data  used  to  compute  propa- 
gation path  losses  is  given  in  Appendix  D. 
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The  propagation  loss  for  jammer  signals  to  the 
ground-based  user  is  illustrated  in  Fig.  3.1-1.  Two  curves 
are  presented;  one  for  the  ground  jammer  and  one  for  the 
airborne  jammer.  Also  illustrated  is  the  propagation  loss 
that  occurs  in  free  space.  Within  line-of-sight , the  total 
propagation  loss  can  be  less  than  the  free  space  loss  due  to 
the  reinforcement  of  the  direct  wave  by  the  reflected  wave. 
Beyond  line-of-sight,  propagation  loss  increases  rapidly, 
although  high-powered  jammers  may  still  be  effective.  Fur- 
ther details  on  the  origin  of  these  curves  are  given  in 
Appendix  D.  These  propagation  loss  characteristics  will  be 
used  in  the  evaluation  of  GPS  manpack  vulnerability  to  elec- 
tronic jamming,  in  Section  3.3.  These  curves  assume  a smooth 
earth's  surface  and  do  not  include  attenuation  due  to  foli- 
age . Detailed  data  on  foliage  attenuation  at  the  GPS  navi- 
gation signal  frequencies  is  presently  not  available. 


Figure  3.1-1  Propagation  Path  Loss 
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3.2  THREAT  COUNTERING  TECHNIQUES 


In  a severe  jamming  environment,  the  GPS  user  may  be 
able  to  emoloy  threat  countering  techniques  to  reduce  his 
vulnerability  to  jamming.  Several  of  these  techniques  are 
discussed  and  evaluated  in  this  section.  It  is  assumed  that 
the  GPS  satellite  power  and  navigation  signal  structure  are 
fixed.  The  techniques  considered  in  this  discussion  are  re- 
stricted to  those  a user  can  influence  through  design  of  the 
receiver  or  through  design  of  auxiliary  equipment.  The  for- 
mer techniques  are  classified  as  internal  to  the  receiver, 
while  the  latter  are  classified  as  external  to  the  receiver. 

3.2.1  Internal 


The  baseline  GPS  user  receiver  characteristics  as- 
sumed for  this  discussion  are  specified  in  Ref.  24.  There 
are  several  modifications  which  can  be  made  to  this  baseline 
receiver  to  reduce  its  vulnerability  to  electronic  jamming 
by  increasing  its  thresnold  jammer  power-to-signal  power 
ratio.  Two  techniques  that  appear  promising  are  notch  fil- 
ters and  reduced  tracking  filter  bandwidth. 

Notch  Filter  - This  technique  reduces  the  jamming 
or  interference  signal  by  placing  a "notch"  in  the  passband 
of  tne  receiver  s input  filter  (hence  the  term  notch  fil- 
ter) at  the  jamming  signal  frequency,  thus  attenuating  the 
jamming  signal.  If  the  notch  is  very  narrow,  only  a small 
amount  of  the  desired  signal  energy  would  be  reduced,  with 
a small  degradation  of  user  performance  (Ref.  26).  This 
technique  is  most  effective  against  a CW  jammer,  allowing 
a 30  to  40  dB  reduction  in  jammer  power.  Notch  filters 
require  a mechanism  for  acquiring  and  tracking  the  jamming 
signal.  The  most  promising  mechanism  for  doing  this  is  a 
phase-lock-loop  (Ref.  26).  However,  this  adds  complexity 
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to  the  receiver  and  may  increase  its  cost  substantially. 
Furthermore,  the  jammer  may  defeat  a notch  filter  by  modu- 
lating his  CW  signal,  although  this  increases  the  jammer  com- 
plexity and  may  decrease  its  effectiveness.  A detailed  in- 
vestigation of  this  technique  of  threat  countering  has  yet 
to  be  performed. 

Reduced  Bandwidth  — The  bandwidths  of  the  baseline 
GPS  receiver  in  the  various  operating  modes  are  chosen  to 
provide  reasonable  performance  for  dynamic  users.  However, 
the  manpack  user  will  normally  be  stationary,  so  that  his 
receiver  bandwidth  could  be  narrower  than  that  required  for 
dynamic  users.  The  bandwidth  of  the  carrier  tracking  loop 
could  easily  be  reduced.  A representative  proposed  design 
calls  for  a 21  Hz  carrier  tracking  loop  bandwidth  (Ref.  27). 
This  bandwidth  could  be  reduced  to  1 or  2 Hz  (and  perhaps 
less*)  for  the  manpack  user.  Reducing  the  carrier  tracking 
loop  bandwidth  to  1 Hz  from  21  Hz  would  result  (Ref.  1)  in 
a 13  dB  increase  in  the  receiver's  threshold  jammer-to- 
signal  power  ratio.  (See  Section  3.3.1  for  further  dis- 
cussion of  this  ratio.)  However,  some  receiver  designs  may 
not  allow  this  much  improvement,  while  others  may  provide 
for  an  even  greater  improvement.  Thus,  for  evaluation  pur- 
poses, the  range  of  10  to  15  dB  will  be  used.  Reducing 
the  receiver  bandwidth  during  initial  acquisition  would 
result  in  a longer  navigation  signal  acquisition  time  per- 
iod. For  example,  using  the  search  strategy  discussed  in 
Ref.  1,  reducing  initial  acquisition  bandwidth  from  a nom- 
inal 1 kHz  to  10  Hz  would  result  in  the  acquisition  time 
increasing  by  a factor  of  100.  Although  this  would  result 
in  a 20  dB  increase  in  interference  margin,  the  time  factor 
is  considered  to  be  unrealistic.  Thus,  reducing  the  receiver 
bandwidth  during  initial  acquisition  is  not  considered 
feasible . 

♦Usable  bandwidths  are  highly  dependent  on  oscillator  in- 
stabilities . 
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3.2.2  External 


Several  techniques  have  been  suggested  for  reducing 
user  vulnerability  with  auxiliary  equipment.  The  basic  idea 
is  to  reduce  the  sensitivity  of  the  GPS  user  antenna  to  the 
directional  interference  or  jamming  signal,  thus  reducing 
the  power  level  of  this  interference  at  the  receiver  input. 
The  techniques  to  accomplish  this  that  are  considered  in 
this  discussion  are  directional  antennas  and  RFI  screens. 

Directional  Antennas  — The  user's  vulnerability  to 
jamming  can  be  reduced  if  his  antenna  is  insensitive  to  sig- 
nals coming  from  the  direction  of  the  interference  or  jam- 
ming. This  can  be  accomplished  through  use  of  directional 
antennas  such  as  parabolic  reflectors  or  antenna  arrays,  as 
discussed  in  Chapter  2.  If  the  parabolic  antenna  can  be 
positioned  so  that  the  interference  is  outside  the  main 
lobe,  then  typically  the  interference  power  level  will  be 
reduced  20  to  25  dB  (Ref.  15),  relative  to  that  of  the  de- 
sired navigation  signal.  If  an  antenna  array  with  null 
steering  were  employed,  the  interference  power  level  could 
be  reduced  35  to  40  dB,  as  discussed  in  Section  2.3.  Thus, 
depending  upon  the  type  of  directional  antenna  deployed, 
reductions  in  the  receiver  jamming-to-signal  power  ratios 
of  25  to  40  dB  are  possible.  Multiple  jammers  may  present 
problems.  However,  an  antenna  array  with  null  steering 
can  be  effective  against  these  multiple  jammers  at  the 
expense  of  antenna  gain  and  increased  processor  complex- 
ity. 


RFI  Screens  — If  the  direction  of  the  interference 
were  known,  an  RFI  screen  could  be  placed  between  the  user 
antenna  and  the  interference.  This  screen  would  act  to 
reduce  the  power  level  of  the  interference,  the  amount  of 
reduction  depending  upon  the  size  of  the  screen  and  distance 
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between  this  screen  and  the  user  antenna.  According  to 
knife-edge  diffraction  loss  theory  (Ref.  15),  if  this  screen 
were  at  least  3 m high  and  10  m long  with  the  user  6 m from 
the  screen,  the  interference  power  level  would  be  reduced 
by  approximately  25  dB.  In  practice,  the  actual  loss  may 
be  less  than  this,  due  to  reflections  from  buildings,  hills, 
etc. , but  the  size  of  the  screen  could  be  increased  to  com- 
pensate. Empirical  data  (Ref.  15)  indicates  that  if  the 
user  antenna  is  behind  a hill  15  m high,  the  loss  in  inter- 
ference signal  strength  would  be  about  23  dB.  Thus,  a 20 
to  25  dB  loss  with  an  RFI  screen  may  be  feasible.  However, 
deployment  of  such  a screen  may  be  difficult  and  thus  limit 
its  usefulness.  This  was  not  investigated  in  this  study. 

3.2.3  Threat  Countering  Effectiveness 

The  effectiveness  of  the  threat  countering  techniques 
discussed  in  the  preceding  sections  is  summarized  in  Table 
3.2-1.  Threat  countering  effectiveness  is  measured  by  the 
amount  of  reduction  that  is  possible  in  the  user  receiver 
jamming-to-signal  power  ratio.  Several  of  these  techniques 
could  be  combined  to  provide  further  reduction  in  interference 

TABLE  3.2-1 

THREAT  COUNTERING  EFFECTIVENESS 


TECHNIQUE 

THRESHOLD  JAMMING-TO-SIGNAL 
POWER  RATIO  IMPROVEMENT  (dB) 

Notch  Filter 

30-40 

Reduce  Bandwidth 

10-15 

Directional  Antennas 

20-25 

Null  Steering  Antennas 

35-40 

RFI  Screens 

20-25 
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levels.  Further  evaluation  of  the  effectiveness  of  these 
threat  countering  techniques  will  be  made  in  the  following 
section,  where  the  reduction  in  the  range  at  which  the  re- 
ceiver is  jammed  will  be  illustrated. 


3.3  USER  VULNERABILITY  EVALUATION 

The  vulnerability  of  a GPS  manpack  user  to  elec- 
tronic jamming  is  evaluated  in  this  section  for  the  scenario 
and  models  of  Section  3.1.  The  criteria  used  for  this  eval- 
uation is  the  critical  Jammer  range.  This  is  defined  as  the 
range  (distance)  from  the  jammer  within  which  a user  would 
be  denied  access  to  the  GPS  navigation  signals,  since  the 
jamming  signal  is  effectively  "overpowering"  the  naviga- 
tion signals  in  the  user  receiver. 

The  jammer  or  interference  level  at  the  user  re- 
ceiver is  computed  in  Section  3.3.1  using  the  propagation 
path  losses  given  in  Section  3.1.3  for  the  ground  and  air- 
borne jammers.  The  critical  jammer  range  is  then  computed 
from  this  data  and  presented  in  Section  3.3.2.  The  effec- 
tiveness of  the  threat  countering  techniques  discussed  in 
Section  3.2  in  reducing  the  critical  jammer  range  is  dis- 
cussed in  Section  3.3.3.  Finally,  the  results  of  the  eval- 
uation are  summarized  in  Section  3.3.4. 

3.3.1  Interference  Level 


The  GPS  navigation  signal  structure  provides  in- 
herent interference  rejection  because  of  the  spread  spectrum 
modulation  technique  used.  However,  the  amount  of  inter- 
ference rejection  available  also  depends  upon  the  receiver 
and  the  particular  operating  mode  employed.  There  are 
basically  three  operating  modes  for  a GPS  receiver: 
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• Initial  signal  acquisition 

• Data  acquisition 

• Carrier  tracking  (normal  navigation  mode) 


Each  of  these  operating  modes  will  have  a particular  interfer- 
ence rejection  capability  depending  upon  the  receiver  design. 
The  first  of  these  modes  usually  has  the  least  amount  of  in- 
terference rejection/,  the  second  an  intermediate  capability, 
while  the  third  has  the  greatest.  Thus,  subsequent  discus- 
sion only  considers  the  interference  rejection  capability 
of  these  two  extremes,  that  is,  initial  signal  acquisition 
and  carrier  tracking.  The  data  acquisition  interference  re- 
jection capability  occurs  within  these  extremes. 


A measure  of  receiver  interference  rejection  capa- 
bility is  the  threshold  jammer-to-signal  power  ratio.  As 
long  as  the  power  of  the  jamming  signal  is  such  that  this 
power  ratio  (or  threshold)  is  not  exceeded,  the  receiver 
will  operate  in  the  mode  that  corresponds  to  that  threshold. 
However,  if  the  jamming  signal,  in  relation  to  the  naviga- 
tion signal  level,  is  such  that  a particular  threshold  is 
exceeded,  the  receiver  will  not  operate  in  the  correspond- 
ing mode.  Thus,  the  threshold  jammer-to-signal  power  ratio 
is  an  effective  measure  of  user  vulnerability  to  jamming. 


The  jammer-to-signal  power  ratio  (J/S)  at  the  re- 
ceiver input  is  defined  as  follows  (all  terms  are  in  dB) 
(Ref.  1): 


(§)  - PJ  + Lp  + Guj  - S dB  (3-3-1> 

where 

P.  is  the  effective  isotropic  radiated 
J power  (EIRP)  of  the  jammer  in  the 
direction  of  the  user  (includes 
jammer  antenna  gain) 
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L is  the  propagation  path  loss  (given 

p in  Fig.  3.1-1) 

G . is  the  gain  of  the  user  antenna  in 

J the  direction  of  jammer 

S is  the  navigation  signal  power  at 
the  receiver 

The  GPS  system  specifications  (Ref.  28)  require  the  P-code 
navigation  signal  level  (S)  at  the  user  antenna  equal  or 
exceed  -163  dBW.  Assuming  the  user  antenna  is  isotropic 
so  that  G . = 0 dB,  then  Eq.  (3.3-1)  becomes 

“ J 


P.  + L + 163 
j P 


dB 


(3.3-2) 


This  jammer-to-signal  power  ratio  for  a 1W  (EIRP)  jammer 
(i.e.,  P = 1W  or  0 dBW)  is  illustrated  in  Fig.  3.3-1.  For 
higher  jammer  powers,  the  vertical  scale  moves  down  by 
10  dB  for  every  order  of  magnitude  increase  in  the  jammer 
power  over  1W.  Note  that  this  figure  is  the  same  as 
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Fig.  3.1-1,  except  the  vertical  scale  has  been  increased  by 
a factor  of  +163  dB  to  reflect  the  value  of  S in  Eq.  (3.3-2). 

The  effectiveness  of  the  jammers  tends  to  diminish 
rapidly  beyond  line-of-sight , especially  for  the  airborne 
jammer.  The  effectiveness  of  the  ground  jammer  decreases 
more  rapidly  than  that  of  the  airborne  jammer  due  to  the 
higher  interaction  between  direct  and  reflected  waves  in 
the  former  case. 

Because  the  army  manpack  user  may  be  jammed  most 
of  the  time  in  a battlefield  environment,  his  receiver's 
most  critical  operating  mode  is  the  initial  acquisition  of 
the  navigation  signal.  In  this  mode  the  receiver  band- 
width is  quite  wide  because  of  the  initial  large  uncer- 
tainty in  the  navigation  signal  frequency.  Thus,  the  re- 
ceiver interference  rejection  capability  or  threshold 
jammer-to-signal  power  ratio  is  lowest  in  this  operating 
mode.  A minimum  value  for  this  power  ratio  is  J/S  = +34  dB 
(Ref.  24).  In  Fig.  3.3-1,  it  can  be  seen  that  this  value 
of  J/S  is  exceeded  whenever  the  user  is  within  7 km  of  the 
1W  ground  jammer  or  90  km  of  the  1W  airborne  jammer. 

Larger  values  of  jammer  power  will  extend  these  ranges  as 
will  be  illustrated  in  the  next  section. 

If  the  jammer  signal  appears  while  the  receiver  is 
in  the  tracking  mode,  the  critical  jammer  range  is  less 
than  that  of  the  acquisition  mode.  The  minimum  threshold 
jammer-to-signal  power  ratio  for  a receiver  in  the  track- 
ing mode  is  J/S  = 49  dB  (Ref.  24),  since  the  receiver 
carrier  bandwidth  is  very  narrow.  This  value  of  J/S  is 
exceeded  whenever  the  receiver  is  within  3 km  of  the  1W 
ground  jammer  and  16  km  of  the  1W  airborne  jammer.  Thus, 
the  critical  jammer  range  is  decreased  significantly  when 
the  receiver  is  in  the  tracking  mode. 
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In  the  next  section,  the  critical  jammer  range  is 
illustrated  as  a function  of  jammer  power  for  the  two  re- 
ceiver operating  modes  discussed  above. 

3.3.2  Critical  Jammer  Range 

The  discussion  in  the  preceding  section  illustrated 
the  dependence  of  the  user's  vulnerability  on  receiver  operat- 
ing mode.  The  dependence  of  user  vulnerability  on  the  jam- 
mer power  is  illustrated  in  this  section.  The  parameter  of 
interest  is  the  critical  jammer  range  which  is  the  range 
from  the  jammer  at  which  a given  threshold  jammer-to-signal 
power  ratio  occurs.  Two  threshold  jammer-to-signal  power 
ratios  are  of  interest,  that  for  initial  acquisition,  which 
is  assumed  to  be  J/S  = 34  dB,  and  that  for  carrier  tracking, 
which  is  assumed  to  be  J/S  = 49  dB  (Ref.  24). 

The  critical  jammer  range  as  a function  of  jammer 
power  for  these  two  threshold  jammer-to-signal  power  ratios 
and  the  ground  and  airborne  jammers  is  illustrated  in  Fig. 
3.3-2.  The  ground  jammer  has  the  lowest  critical  jammer 
range  as  noted  in  the  previous  section.  An  especially 
interesting  result  from  this  curve  is  that  jammer  powers 
greater  than  1000  W does  not  result  in  a significant  increase 
in  the  critical  jammer  range  for  the  airborne  jammer. 

3.3.3  Vulnerability  Reduction 

Although  the  GPS  receiver  has  an  inherent  inter 
ference  rejection  capability,  a user  with  the  baseline 
receiver  (especially  the  Army  manpack  user)  is  still  vul- 
nerable to  jamming  with  relatively  simple  jammers,  as  the 
previous  sections  have  shown.  Therefore,  threat  counter- 
ing techniques  are  extremely  important.  Since  these 
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JAMMER  POWER  EIRP  (watts) 


Figure  3.3-2  Critical  Jammer  Range 

techniques  effectively  increase  the  user  receiver  threshold 
jammer-to-signal  power  ratio,  they  can  decrease  the  criti- 
cal jammer  range.  Assuming  1 kW  EIRP  jammers,  the  critical 
jammer  ranges  for  users  with  and  without  the  threat  counter- 
ing techniques  discussed  in  Section  3.2  are  summarized  in 
Table  3.3-1.  Critical  jammer  ranges  are  given  for  the  re- 
ceiver in  the  initial  acquisition  and  carrier  tracking  modes 
for  both  the  ground  and  airborne  jammers.  The  critical  jam- 
mer range  for  the  combined  utilization  of  reduced  receiver 
bandwidth  and  RFI  screen  is  also  given. 

These  results  indicate  that  against  a ground-based 
jammer,  an  RFI  screen  and/or  reduced  bandwidth  are  about  as 
effective  as  the  other  techniques.  Against  an  airborne 
jammer,  the  notch  filter  and  null  steering  antenna  are  the 
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TABLE  3.3-1 


THREAT  COUNTERING  EFFECTIVENESS 

/ 

J 


CRITICAL  JAMMER  RANGE  (km) 

THREAT 

GROUND  JAMMING  (1  kW  EIRP) 

AIRBORNE  JAMMER  (lkW  EIRP) 

COUNTERING 

TECHNIQUE 

Initial 
Acquisition 
(J/S  - 34  dB) 

Carrier 
Tracking 
(J/S  - 49  dB) 

Initial 
Acquisition 
(J/S  » 34  dB) 

Carrier 
Tracking 
(J/S  - 49  dB) 

Baseline 

24 

14 

210 

170 

Notch  Filter 

4-7 

2-3 

30-90 

5-16 

Reduced  Bandwidth 

N/A 

7-9 

N/A 

90-120 

Directional  Antenna 

9-11 

4-5 

120-150 

30-50 

Null  Steering  Antenna 

4-5 

2 

30-50 

5-8 

RFI  Screen 

9-11 

4-5 

120-150 

30-50 

Reduced  Bandwidth 
and  RFI  Screen 

9-11 

2-3 

120-150 

5-16 

most  effective.  Reducing  the  bandwidth  of  the  carrier  loop  in 
the  receiver  is  effective,  and  is  the  simplest  threat  count- 
ering technique  for  a stationary  user,  although  it  can  only 
be  used  when  the  receiver  is  in  the  tracking  mode. 

3.3.4  Vulnerability  Evaluation  Summary 

The  GPS  user  is  quite  vulnerable  to  electronic  jam- 
ming from  both  a ground-based  and  an  airborne  jammer.  This 
is  especially  true  for  the  Army  manpack  user  since  he  will 
most  likely  be  attempting  to  acquire  the  GPS  navigation  sig- 
nals when  jamming  is  present.  The  results  of  this  study 
indicate  that  in  the  initial  acquisition  mode  with  a ground 
jammer  radiating  1 kW  EIRP,  the  user  will  not  be  able  to 
operate  within  24  km  of  the  jammer  without  special  threat 
countering  equipment.  If  the  user  were  behind  an  RFI  screen 
of  3 m height,  this  critical  jammer  range  could  be  reduced 
to  approximately  10  km. 
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If  the  jammer  is  airborne,  the  user  is  even  more 
vulnerable.  In  the  initial  acquisition  mode,  the  critical 
jamming  range  is  approximately  210  km.  The  most  effective 
threat  countering  technique  in  this  case  is  either  the  notch 
filter  or  the  null  steering  antenna. 


4. 


RECEIVER  CLOCK  ERROR  ANALYSIS 


This  chapter  develops  a model  describing  environ- 
mentally-induced errors  in  a quartz  crystal  oscillator.  The 
error  model  is  utilized  to  assess  the  degradation  in  GPS 
navigation  performance,  for  an  airborne  (helicopter)  user, 
arising  from  the  use  of  the  receiver  oscillator  (clock)  in  a 
realistic  (non-ideal)  environment. 


4.1  RECEIVER  CLOCK  ENVIRONMENTAL  ERROR  MODEL 
4.1.1  Introduction 

To  properly  assess  the  impact  of  environmentally- 
induced  quartz  crystal  oscillator  errors  on  GPS  navigation 
performance,  a comprehensive  model  of  the  oscillator  error 
behavior  must  be  developed.  This  section  describes  the  de- 
velopment of  an  error  model  which  considers  the  effects  of 
the  environment  on  a quartz  oscillator.  The  complete  error 
model  will  include  a description  of  the  intrinsic  stability 
of  the  oscillator  as  well  as  the  environmentally-induced 
errors.  The  model  for  intrinsic  stability  follows  directly 
from  a consideration  of  the  frequency  stability  of  a specific 
oscillator.  (This  model  is  developed  in  Appendix  C,  which 
also  includes  background  discussions  of  frequency  stability 
and  the  modeling  of  flicker  noise.)  The  Hewlett  Packard 
(HP)  10544A  quartz  crystal  oscillator  will  be  utilized 
throughout  this  study  as  a representative  device.  The 
model  describing  environmental  effects  is  based  upon  en- 
vironmental performance  data  supplied  by  Hewlett  Packard 
and  the  U.S.  Army  Electronics  Command. 
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The  prototype  for  the  quartz  crystal  oscillator  er- 
ror model  to  be  considered  is  shown  in  Fig.  4.1-1  with 


6<f>(t) 


«4>(t0) 

6f(t0) 


6f(t0) 


total  clock  phase  error  in  seconds 
at  time  t 

initial  phase  error  (seconds) 

initial  fractional  frequency  off- 
set (dimensionless) 

initial  fractional  frequency  drift 
or  "aging"  (seconds-*) 

spectral  level  of  white  noise 
driving  the  flicker  noise  model 
( dimensionless ) 


Figure  4.1-1  Prototype  Quartz  Crystal  Oscillator 

Error  Model 

The  initial  phase  and  fractional  frequency  errors, 

6 4> ( 1 0 ) and  6f(tQ),  respectively,  are  dependent  upon  the  spe- 
cific application  in  which  the  oscillator  is  to  be  used.  For 
this  reason,  specific  values  for  these  random  variables  will 
be  deferred  until  the  discussion  of  the  performance  projec- 
tions in  Section  4.2.  The  initial  condition  characterizing 
the  "aging"  of  the  oscillator  can  be  specified  and  is  com- 
puted in  Section  C.3  along  with  the  spectral  level  (q^)  of 


the  white  noise  in  the  flicker  noise  model.  The  environmental 
effects  model  is  developed  in  the  next  section. 

4.1.2  Modeling  of  Environmental  Effects 

Any  complete  time  standard  error  model  must  contain 
a description  of  the  effects  the  environment  has  on  the  per- 
formance of  the  oscillator.  A completely  general  error  model 
for  environmental  effects  would  be  prohibitively  complex. 
However,  if  certain  assumptions  $:re  made  about  the  form  and 
number  of  environmental  influences,  then  an  environmental 
effects  model  can  be  developed  for  a specific  device.  In 
the  following  paragraphs,  the  effects  of  temperature,  warm- 
up, acceleration,  vibration  and  shock  on  a quartz  crystal 
oscillator  (HP  10544A)  will  be  discussed. 

Temperature  - A quartz  oscillator  responds  to  a 
change  in  ambient  temperature  with  a shift  in  frequency. 

The  magnitude  of  this  nonlinear  effect  is  shown  in  Fig.  4.1-2 
for  an  HP  10544A  oscillator.  (The  fractional  frequency  change 
is  shown  as  an  excursion  from  the  oscillator  frequency  at 
30°C.)  When  the  oscillator  is  utilized  in  a GPS  receiver,  a 
linearized  coefficient  relating  the  frequency  shift  to  a 
temperature  change  can  be  obtained.  For  temperature  excur- 
sions about  the  nominal  operating  temperature  of  the  GPS  re- 
ceiver (assumed  to  be  =60°C,  from  Ref.  7),  this  temperature 
drift  coefficient  (Ktv)  is  linear  with  a maximum  value  (see 
oscillator  Serial  No.  1410),  K^v  = 1 x 10-^^/°C.  Figure  4.1-2 
illustrates  the  wide  variation  in  temperature  drift  coeffi- 
cients  for  differing  devices  of  the  same  type. 

i 

The  frequency  shift  which  results  from  a tempera- 
ture change  does  not  occur  instantaneously  due  to  a "thermal 
lag."  The  temperature-regulating  oven  surrounding  the  quartz 
crystal  is  the  main  contributor  to  this  lag.  The  thermal 
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Figure  4.1-2  Frequency-Temperature  Characteristic  of 

Selected  HP  10544A  Oscillators  (Ref.  6) 


response  shown  in  Fig.  4.1-3  can  be  adequately  described  by 
the  simple  first-order  system  shown  in  Fig.  4.1-4.  The  time 


TIME  (min) 

Figure  4.1-3  Response  of  Crystal  Oscillator  Frequency  to 

a Step  Ambient  Temperature  Change  (Ref.  7) 


R*  13971* 


FREQUENCY 

ERROR 


Figure  4.1-4 


Temperature  Variation  Error  Model 


constant,  , depends  upon  the  number  of  ovens  surrounding 

the  crystal.  For  a single  oven  device  (such  as  the  HP  10544A), 

this  time  constant  (t.  ) is  on  the  order  of  80  minutes  (see 

v tv'  v 

Fig.  4.1-3). 

The  error  model  shown  in  Fig.  4.1-4  represents  the 
response  of  the  crystal  oscillator  to  an  arbitrary  tempera- 
ture variation.  To  utilize  this  model,  the  input  temperature 
variation  must  be  described  quantitatively.  It  is  obvious 
that  differing  users  (e.g.,  manpack,  helicopter,  winged  air- 
craft) will  encounter  vastly  different  environmental  condi- 
tions; e.g.,  a manpack  user  moves  from  a heated  room  to  a 
subfreezing  field  location,  an  aircraft  climbs  from  sea 
level  to  cruising  altitude,  etc.  These  dissimilar  environ- 
mental scenarios  will  result  in  mission-dependent  tempera- 
ture fluctuations  outside  the  GPS  receiver.  These  variations 
are  not  transmitted  instantaneously  to  the  crystal  oscilla- 
tor due  to  the  moderating  influence  of  the  receiver  thermal 
mass.  Assuming  random  temperature  variations,  the  ambient 
temperature  surrounding  the  crystal  oscillator  is  reasonably 
modeled  as  a first-order  Markov  process.  The  parameters  of 
the  process,  rms  value  and  time  constant,  must  be  chosen  to 
reflect  the  temperature  variations  outside  the  GPS  receiver 
and  the  moderating  effect  of  the  receiver.  A composite 
thermal  error  model  is  given  in  Fig.  4.1-5.  For  a GPS 
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Figure  4.1-5  Quartz  Crystal  Therrral  Error  Model 

receiver,  the  rms  value  of  the  temperature  fluctuations  near 
the  crystal  oscillator  is  expected  to  be  on  the  order  of  3°C 
with  a time  constant  (t  ) of  100  minutes  (based  on  Ref.  7). 

Si 

The  corresponding  spectral  level  of  the  white  noise  (q&)  is 
3 x 10-3  °C2/sec . 

Warm-up  - Figure  4.1-6  shows  the  fractional  fre- 
quency error  versus  warm-up  time  characteristics  of  three 
HP  10544A  quartz  crystal  oscillators.  The  three  oscillators 
exhibit  vastly  different  behavior  which  is  not  amenable  to 
direct  analysis.  The  dashed  line  in  the  figure  describes 


Figure  4.1-6  Warm-up  Characteristic  of  HP  10544A 

Quartz  Crystal  Oscillator  (Ref.  6) 


an  exponentially  decaying  envelope  (the  plot  is  semi-log) 

which  can  be  viewed  as  a statistical  bound  on  the  frequency 

error.  The  warm-up  fractional  frequency  error  variance 

q 

(°6fw)  *”or  ^is  exPonential  model  is 


2 

°6fw 


-2t/t 


2 w 

o e 
w 


(4.1-1) 


2 -10 

with  o = 1 x 10  and  t ^ 100  seconds.  The  warm-up  time 
w w 

constant,  xw>  is  much  less  than  x since,  during  warm-up, 
the  response  is  primarily  that  of  the  crystal  without  the 
moderating  effect  of  the  oven. 


Acceleration  — A constant  acceleration  introduces 

a frequency  offset  in  a quartz  oscillator  proportional  to 

the  magnitude  of  the  acceleration,  i.e.,  6f  = K |a| . The 

a S 

coefficient  of  the  proportionality,  K , is  dependent  upon 

o 

the  orientation  of  the  acceleration  vector  with  respect  to 

the  crystallographic  axes  of  the  quartz  crystal.  A worst 

_g 

case  value  is  K = 1 x 10  /g  (Ref.  7). 

Vibration  - Vibration  introduces  sidebands  (at  the 
vibration  frequency)  into  the  spectrum  of  the  signal  from 
a quartz  crystal  oscillator.  These  sidebands  produce  a 
frequency  error  which  is  dependent  upon  the  oscillator  sen- 
sit  vity  at  the  vibration  frequency.  Measurements  have  been 
made  of  the  fractional  frequency  error  (averaged  over  a 
10-second  interval)  arising  from  the  lateral  vibration  of 
a quartz  crystal  oscillator  (Ref.  7).  Figure  4.1-7  shows 
the  measured  response  to  a sinusoidvl  input  similar  to  the 
vibration  specification  (Curve  M)  in  MIL-STD-810B.  The 
worst  case  response  is  adequately  modeled  by  the  resonant 
second-order  curve  which  is  superimposed  upon  the  data. 

In  lieu  of  a well-defined  spectrum  for  the  vibration  in  an 
operational  environment,  a baseline  model  will  be  developed 
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Figure  4.1-7  Frequency  Error  as  a Function  of 

Vibration  Frequency 

assuming  a wide  band  noise  input.  The  wide  band  noise  can 

be  represented  as  a first-order  Markov  process  with  a cor- 

-4 

relation  time,  x , of  3.2  x 10  seconds.  The  power  spec- 

m 

trum  of  this  process  is  essentially  flat  out  to  a frequency 

of  500  Hz.  This  wide  band  noise  is  undistinguishable  from 

white  noise  up  to  a vibration  frequency  of  300  Hz  which  is 

beyond  the  region  of  crystal  response  shown  in  Fig.  4.1-7. 

A worst  case  value  (cf,  MIL-STD-810B , Curve  AE)  for  the  rms 

level  of  the  input  noise  is  = 5.4g.  The  wide  band  noise 

m 

can  be  simplified  to  an  equivalent  white  noise  (as  shown  in 
Ref.  8)  with  a spectral  level 

q = 2t  = 0.02  g^  sec  (4.1-2) 

v mm 

The  model  representing  the  frequency  error  due  to 
vibration,  5fv,  is  shown  in  Fig.  4.1-8.  The  input  is  the 
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Figure  4.1-8  Crystal  Oscillator  Vibration  Error  Model 

white  noise  equivalent  to  the  assumed  wide  band  noise  and 
the  second-order  system  represents  the  response  of  the  oscil- 
lator as  shown  in  Fig.  4.1-7.  The  resonant  curve  in  the 
figure  corresponds  to  a natural  frequency,  u>n , of  754  rad/sec 

and  a damping  ratio,  £ = 0.1.  The  constant  K is  the  lateral 

v -10 

vibration  drift  coefficient  with  a value  of  6 x 10  /g  from 
Fig.  4.1-7.  As  expected,  the  vibration  coefficient  Kv  is 
very  close  to  the  g-sensitive  drift  coefficient  K . 

Shock  — When  a quartz  crystal  oscillator  is  sub- 
jected to  a shock,  a permanent  frequency  offset  results. 

No  specific  values  are  available  (at  this  time)  for  the 

magnitude  of  this  offset  during  the  shock,  but  step  frac- 

-9 

tional  frequency  shifts  on  the  order  of  1 x 10  have  been 
observed  (Ref.  7)  following  shocks  of  15g,  45g,  and  60g. 

A reasonable  model  for  the  effect  of  a shock  of  short  dura- 
tion (~  1 msec),  occurring  at  time  t , is  a step  fractional 

i -9  s 

frequency  shift  ( 6 f | ) of  1 x 10 

s t 

s 

4.1.3  Composite  Quartz  Oscillator  Error  Model 

In  the  previous  section,  models  were  developed  to 
describe  the  frequency  errors  which  arise  from  perturbations 
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in  the  oscillator's  environment.  The  total  fractional  fre- 


quency error  due  to  environmental  effects,  6f  , is  the  sum 
of  the  individual  errors  resulting  from  temperature  changes, 
oscillator  warm-up,  acceleration,  shock  and  vibration  ef- 
fects, viz, 


6f  = 6f.  + 
env  t 


5 f + 5f  + 
w a 


6f 


+ 6f 


(4.1-3) 


The  environmental  effects  model  can  be  combined  with 
the  intrinsic  stability  model  of  the  quartz  oscillator  (see 
Appendix  C)  as  shown  in  Fig.  4.1-1.  The  state  vector  for 
the  composite  quartz  crystal  oscillator  error  model  is  given 
in  Table  4.1-1.  Two  error  sources  are  not  included  explicitly 
in  the  state  vector.  The  error  arising  from  the  acceleration 
of  the  oscillator  (i.e. , 6f  = K |a|)  is  more  appropriately 
modeled  as  an  additive  quantity.  The  second  error  source 


TABLE  4.1-1 

QUARTZ  CRYSTAL  OSCILLATOR 
ERROR  MODEL  STATE  VECTOR 


STATE 

SYMBOL 

DESCRIPTION 

1 

6$ 

Phase  Error 

2 

6f 

Frequency  Error 

3 

6f 

Frequency  Rate  Error 

4 

xlf 

Flicker  Noise 

5 

X2f 

Intermediate  Flicker  Noise  State 

6 

6ft 

Frequency  Error  (Thermal  Model  ) 

7 

xt 

Temperature  Variation 

8 

6fw 

Frequency  Error  (Warm-up  Model) 

9 

fifv 

Frequency  Error  (Vibration  Model) 

10 

xv 

Intermediate  Vibration  Error  State 
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which  will  remain  unmodeled  is  the  effect  of  an  impulsive 
shock  to  the  oscillator.  In  any  practical  implementation  of 
a Kalman  filter,  a priori  knowledge  of  the  time  of  occurrence 
of  a shock  is  unlikely  and  such  behavior  cannot  be  described 
by  a Gauss-Markov  random  process  model.  However,  the  effect 
of  such  an  unanticipated  shock  can,  and  will,  be  examined  in 
this  study. 


The  error  model  can  be  written  in  state  variable  form 
as 

x(t)  = F x(t)  + G w(t)  (4.1-4) 

where 

x(t)  = 10  x 1 column  vector  of  the  error 

states  in  Table  4.1-1 


F 

ss 

10 

x 10 

error  dynamics  matrix 

w(t) 

= 

10 

X 1 

column  vector  of 

process 

noise 

G 

10  x 10 
matrix 

process 

poise  distribution 

The  dynamics  matrix, 

0 

1 

0 

0 

l 

0 

• l 

l 

0 

0 
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0 

0 

0 

0 
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0 
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0 
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0 

0 

0 

0 

0 

0 

0 

0 

0 

f2 

1 

0 

0 

0 

0 

0 

0 

0 

0 

f3 

0 

0 

0 
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(4.1-5) 
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is  based  upon  the  prototype  model  shown  in  Fig.  4.1-1,  the 
environmental  models  developed  above  and  the  model  of  intrin- 
sic stability  developed  in  Section  C.3.  Values  of  f , f 

X £t 

and  fg  are  defined  in  Eq.  (C.3-7).  The  noise  distribution 
matrix,  also  based  on  Section  C.3,  is 
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(4.1-6) 


where  glf  gg  and  gg  are  defined  in  Eq.  (C.3-8).  The  contin- 
uous process  noise  matrix,  Q , is  defined  by 

E^w(t)  wT(o]  = Qc  6(t-C)  (4.1-7) 

where  6(t-£)  is  the  Dirac  delta  function.  Table  4.1-2  lists 

the  non-zero  elements  of  Q . 

c 

TABLE  4.1-2 

CONTINUOUS  PROCESS  NOISE  MATRIX 
(NON-ZERO  ELEMENTS) 


ELEMENT 

SOURCE 

SYMBOL 

VALUE 

(1,1) 

Flicker  Noise 

qf 

1.14  x 10-19 

(7,7) 

Temperature  Variation 

qa 

3.0  x 10'3  °C2/sec 

(10,10) 

Vibration 

qv 

2.0  x 10~2  g2  sec 

L 
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4.2 


USER  PERFORMANCE  PROJECTIONS 
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4.2.1  Evaluation  Scenario 


The  analysis  of  the  airborne  user  performance  em- 
ployed the  GPS  Direct  Simulation  Program  (DSIM,  see  Ref.  2) 
in  order  to  investigate  the  effects  of  environmentally-induced 
errors.  The  airborne  user  trajectory  is  illustrated  in  Fig. 
4.2-1.  The  user  initially  processes  one  minute  of  GPS  mea- 
surements while  stationary  to  calibrate  the  receiver  clock. 
After  takeoff,  the  helicopter  makes  turns  at  a maximum  of 
15  deg/sec,  accelerates  at  a maximum  of  2.5  mps/sec  and 
changes  altitude  at  a maximum  of  3.0  mps. 
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Figure  4.2-1  Airborne  User  Trajectory 


Preliminary  investigation  indicated  that  a GPS/inertial 
configuration  would  be  more  sensitive  to  receiver  clock  errors 
than  other  mechanizations  such  as  GPS/air  data  or  GPS  alone. 

Thus,  in  order  to  evaluate  a worst  case  situation,  an  inertial 
navigation  system  (INS)  has  been  adopted  for  this  study.  Al- 
though a GPS/ INS  system  may  not  be  of  primary  Army  interest  at 
this  time,  the  findings  relative  to  clock  modeling  are  generally 
applicable  to  all  GPS  mechanizations  . A barometric  altimeter  is  also 
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included  to  bound  the  vertical  channel  of  the  INS.  The  measured 
altitude  is  used,  in  a closed-loop  fashion,  to  damp  the  alti- 
tude errors.  A detailed  error  model  for  the  INS  is  presented 
in  Section  4.2.2. 

A detailed  discussion  of  the  GPS  Phase  III  config- 
guration  used  in  this  study  is  presented  in  Appendix  B.  For 
the  airborne  user  performance  analyses  considered  in  this 
chapter,  pseudo-range  and  pseudo-range-rate  measurements  are 
obtained  simultaneously  from  four  of  the  satellites  in 
view  (numbers  1,  2,  3 and  5)  every  10  seconds.  Initial  state 
errors,  initial  filter  standard  deviation  values  and  process 
noise  variances  for  this  study  are  given  in  Section  4.2.2. 

4.2.2  Error  Models 

Two  error  models  are  required  for  the  generation  of 
user  performance  projections.  The  first  is  a model  describing 
the  "real  world"  errors.  This  model  would  ideally  be  a com- 
plete description  of  all  significant  error  sources  which 
impact  upon  user  navigation  performance.  The  second  model, 
which  is  a deliberately  simplified  description  of  the  "real 
world"  errors,  is  implemented  in  the  user  Kalman  filter. 

"Real  World"  Error  Model  — The  "real  world"  error 
model  includes  a description  of  GPS  system  errors,  the  user 
INS  errors  and  a user  clock  error  model.  The  GPS  system 
errors,  as  implemented  in  the  DSIM  program,  include  satellite 
clock  and  ephemeris  errors,  signal  propagation  delays  (iono- 
spheric, tropospheric  and  multipath),  and  receiver  noise 
and  quantization.  These  error  sources  and  their  models  are 
well  documented  in  Ref.  2 and  will  not  be  repeated  here.  The 
equations  describing  inertial  system  error  propagation  are 
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explicitly  stated  in  the  discussion  of  the  user  filter  later 
in  this  section.  The  barometric  altimeter  is  used  to  sense 
the  altitude  of  the  user  above  the  earth.  The  error  in  mea- 
suring the  altitude  is  modeled  as  a scale  factor  error  (first- 
order  Markov  process)  plus  a white  measurement  noise  (Ref.  2). 
Figure  4.2-2  is  a block  diagram  of  the  barometric  altimeter 
error  model.  The  altitude  error  is  computed  as 


6h  = h 6s  + v (4.2-1) 

m u a a 


where 

6h  = altitude  error 
m 

h = true  user  altitude 
u 

6s  = barometric  altimeter  scale  factor  error 
a 

v = measurement  noise  associated  with  the 
a altitude  measurement;  E[va(tn)]  = 0 

E [v  (t  ) v (t  )]  = a2  6 ; 6 is  the 

1 av  n7  av  m7  a mn’  mn 

the  Kronecker  delta  function,  i.e., 

6mn  - I1  m * n 

mn  )o  m f n 
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Figure  4.2-2  Barometric  Altimeter  Error  Model 
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The  user  clock  error  model  was  developed  in  Sec- 
tion 4.1.  The  10-state  clock  error  model  includes  a descrip- 
tion of  the  effects  of  temperature,  warm-up  and  vibration. 
Errors  which  arise  from  acceleration  and  shock  are  not  ex- 
plicitly included  in  the  error  state  vector  (see  Table 
4.1-1).  However,  these  two  error  sources  are  included  in  the 
description  of  the  "real  world"  and  their  effects  on  user 
navigation  performance  will  be  evaluated. 

User  Filter  — The  user  Kalman  filter  estimates  the 
20  user  navigation  and  clock  errors  summarized  in  Table  4.2-1. 
The  user  position  and  velocity  errors  are  referenced  to  a 
locally  level  [north  (n),  east  (e)  and  down  ( d ) ] coordinate 
system  located  at  the  instantaneous  user  position.  As  in- 
dicated in  Table  4,2-1,  INS  errors  are  characterized  by  gyro 
drift  (with  bias  and  random  walk  components)  and  accelerometer 
bias  error  terms.  This  model  will  provide  a reasonable  de- 
scription of  INS  errors  during  the  short  (=8  minute)  mission 
under  consideration.  Since  the  purpose  of  this  study  is  to 
evaluate  the  effect  of  user  clock  errors,  a more  detailed  INS 
error  model  would  be  undesirable.  The  user  filter  does  not 
include  a model  of  any  environmentally-induced  errors  in  the 
user  clock.  The  state  modeling  the  effects  of  aging  has  also 
been  deleted  in  the  filter  model,  since  the  total  mission  dura- 
tion is  shorter  than  the  time  required  for  the  effects  of  ag- 
ing to  become  significant. 

The  user  filter  error  model  can  be  written  in  state 
variable  form 


x (t)  = F (t)  x (t)  + G w (t)  (4.2-2) 

— u U — u u — u 

where 

x(t)  = 20  x 1 column  vector  of  user 

error  states 
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STATE 

NUMBER 


SYMBOL 


DESCRIPTION 


1 

6Rn 

North  Position  Error 

2 

6R 

e 

East  Position  Error 

3 

6h 

Altitude  Error 

4 

5Vn 

North  Velocity  Error 

5 

6Ve 

East  Velocity  error 

6 

6Vd 

Down  Velocity  Error 

7 

Platform  Misalignment  About  North 

8 

*e 

Platform  Misalignment  About  East 

9 

♦d 

Platform  Misalignment  About  Down 

10 

6sa 

Barometric  Altimeter  Error 

11 

en 

North  Gyro  Drift 

12 

ee 

East  Gyro  Drift 

13 

ed 

Down  Gyro  Drift 

14 

wn 

North  Accelerometer  Bias  Error 

15 

East  Accelerometer  Bias  Error 

16 

yd 

Down  Accelerometer  Bias  Error 

17 

64> 

Phase  Error 

18 

6f 

Fractional  Frequency  Error 

19 

xlf 

Flicker  Noise 

20 

X2f 

Intermediate  Flicker  Noise  State 

F^(t)  = 20  x 20  user  error  dynamics  matrix 

wu(t)  = 20  x 1 column  vector  of  white  pro- 
cess noise  such  that  E[wu(t)l  = 0, 
Cov[wu( t )]  = E[wu(t)  wuU)TJ 
= Qu(t)  <$(t-c);  6(t-c)  is  the  Dirac 
delta  function. 

G = 20  x 20  process  noise  distribution 

matrix 


The  dynamics  matrix,  F , may  be  partitioned  as 


F (t)  = 
uv  ' 


Fs(t)  ■ 

i 

0 

i 

0 . 

F 

c 

(4.2-3) 


where 

F (t)  = 16  x 16  inertial  navigation  sys- 

tem error  dynamics  matrix.  (This 
description  of  INS  errors  is  also 
included  in  the  "real  world"  error 
model . ) 

Fc  = 4x4  clock  error  dynamics  matrix 


The  F (t)  matrix  (as  developed  in  Refs.  2 and  8)  is  shown 
s 

in  Fig.  4.2-3  with 


V ,V 
n ’ e 


R 


r 

L 


the  user's  north  and  east  vel- 
ocities in  the  local  level  navi- 
gation frame 

position  vector  of  the  user 
measured  from  the  center  of  the 
earth 

111 

user  latitude 


ft  = earth's  angular  rate 

D1,D2  = vert::l-cal  channel  damping  coeffi- 

cients (0.0011  sec“2  and  0.0666 
sec_l  respectively) 
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= altimeter  inverse  time  constant 


f ,f  , f. 
n e ’ d 


nominal  acceleration  due  to 
gravity 

specific  force;  equal  to  a , a , 
ad  " go'  respectively 


The  clock  error  dynamics  matrix,  from  Eq. (4.1-5),  is 


0 1 


1.0  x 10‘ 


0 0 


0 0 -1.01  x 10 


(4.2-4) 


0 0 -1.0  x 10‘ 


The  noise  distribution  matrix  may  also  be  partitioned  such 
that 


G = 
u 


I16xl6  ' 


g 0 

fec  i 


(4.2-5) 


where  I^gx^g  denotes  a 16x16  identity  matrix  and,  from 
Eq.  (4.1-6), 


1.0  x 10‘ 


9.09  x 10 


(4.2-6) 


9.90  x 10' 


The  user  filter  initial  conditions  are  summarized 
in  Table  4.2-z,  The  data  contained  in  the  table  includes 
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AIRBORNE  USER  STATE  AND  FILTER  ERROR  VALUES 
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the  initial  filter  state  estimation  errors,  the  initial  filter 
estimation  error  standard  deviations,  and  values  assumed  for 
the  filter  noise  spectral  densities.  The  user  filter's  ini- 
tial covariance  matrix,  Pu(0),  is  assumed  to  be  diagonal  with 
the  exception  of  the  correlation  between  platform  level  align- 
ment (<p  , i|>  ) and  accelerometer  bias  errors  (y  , y ) and  be- 
tween  azimuth  alignment  (\p  ) and  east  gyro  drift  (e  ).  The 

Z 6 

level  alignment/accelerometer  bias  correlation  was  selected 
as  0.95  and  the  azimuth  alignment  and  east  gyro  drift  were 
assumed  to  be  perfectly  correlated.  This  represents,  in  a 
simplified  manner,  the  result  of  pre-flight  gyrocompass  align- 
ment. The  diagonal  elements  of  the  filter  covariance  matrix, 

P , are  set  equal  to  the  squares  of  the  initial  filter  stan- 
dard deviations.  The  continuous  spectral  density  matrix,  Qu, 
is  assumed  to  be  a diagonal  matrix  with  the  diagonal  equal 
to  the  spectral  densities  given  in  the  table.  The  noise  dis- 
tribution matrix  is  given  by  Eqs.  (4.2-5)  and  (4.2-6). 

The  clock  error  states  (phase,  frequency,  etc.)  in 
the  table  are  expressed  with  the  same  dimensions  as  position 
and  velocity  errors  (i.e.,  m and  mps).  The  conversion  from 
the  units  of  Section  4.1  to  those  considered  here  is  accom- 
plished by  multiplying  the  appropriate  clock  states  in  the 
Pu(0)  and  Qu  matrices  by  the  square  of  the  velocity  of  pro- 
pagation of  the  radio  signal  from  the  satellite  to  the  user. 
This  study  will  consider  the  signal  propagation  velocity  to 

be  constant  with  a value  equal  to  the  speed  of  light, 
ft 

c = 3 x 10  mps.  The  Fu  and  Gu  matrices  remain  unchanged. 

The  user  filter  formulation  assumes  that  the  user 
is  equipped  with  a multichannel  X-model  receiver  to  allow 
the  simultaneous  processing  of  pseudo-range  and  pseudo- 
range-rate information  from  four  of  the  satellites  in  view. 

The  filter-assumed  receiver  measurement  noise  levels  are 
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7.5  m (rms)  and  0.15  mps  (rms)  for  pseudo-range  and  pseudo- 
range-rate, respectively.  Errors  due  to  the  ionospheric 
and  tropospheric  effects  are  assumed  (in  the  filter)  to 
be  compensated  to  within  the  measurement  noise  and  are  not 
correlated  with  the  other  error  sources. 

4.2.3  Results 

This  section  will  present  user  performance  projec- 
tions which  predict  the  impact  of  environmentally-induced 
GPS  receiver  clock  errors  on  user  navigation  performance. 
The  results  are  given  in  the  form  of  user  estimation  errors 
and  the  rms  estimation  error  envelope  (standard  deviation 
boundary)  predicted  by  the  user  Kalman  filter.  Each  error 
source  modeled  in  Section  4.1.2  is  introduced  into  the 
"real  world"  model  separately  and  its  effect  on  user  navi- 
gation performance  is  assessed  by  comparing  the  resultant 
navigation  errors  with  those  of  a baseline  run  in  which 
environmentally-induced  errors  were  deleted.  The  follow- 
ing paragraphs  discuss  the  impact  of  the  various  error 
sources. 


Baseline  — The  baseline  performance  of  the  airborne 
user  is  indicative  of  the  navigation  errors  which  arise 
from  GPS  system  errors  (satellite  ephemeris  and  clock  errors, 
signal  propagation  errors,  etc.)  and  user  dependent  errors 
(INS  error  sources  and  the  non-environmentally-induced  clock 
errors).  This  performance  is  illustrated  in  Fig.  4.2-4 
for  the  trajectory  of  Fig.  4.2-1.  The  horizontal  position 
error  is  under  8 m and  the  altitude  error  is  below  1.5  m. 

The  horizontal  velocity  error  is  less  than  0.1  mps  except 
during  the  first  turn  when  it  peaks  to  0.26  mps.  The  peak 
is  due  to  a residual  heading  error  which  couples  the  accelera- 
tion experienced  during  the  turn  into  the  horizontal  velocity 


TIME  (tec) 


TIME  (tec) 


VERTICAL  MANEUVERS 


2S0 

TIME  (tec) 


Airborne  User  Errors  — Baseline 


X 

- ESTIMATION  ERROR 

a 

- STANOARD  DEVIATION  BOUNDARY 

Q — 

R-1B294 

error.  This  makes  the  heading  error  highly  observable,  such 
that  GPS  pseudo-range-rate  measurements  during  the  first  turn 
calibrate  the  heading  error  (see  Fig.  4.2-5).  The  vertical 
velocity  error  remains  below  0.05  mps.  The  user  clock 
phase  and  frequency  errors  settle  out  to  less  than  1.25  m 
and  0.01  mps,  respectively,  after  processing  250  seconds 
of  pseudo-range  and  pseudo-range-rate  measurements. 

Degraded  Quality  — The  "degraded  quality"  run  con- 
sidered the  case  in  which  the  intrinsic  stability  of  the 
real  world  oscillator  was  worse  than  that  chosen  as  the 
filter  design  model.  The  degraded  quality  oscillator  was 
obtained  by  increasing  the  process  noise  spectral  level  q^ 

[Eq.  (C.3-9)]  by  a factor  of  ten.  The  performance  projections 
for  the  degraded  quality  clock  are  shown  in  Fig.  4.2-6.  The 
horizontal  position,  altitude,  horizontal  velocity  and  verti- 
cal velocity  errors  remain  essentially  unchanged  from  their 
baseline  values.  The  user  clock  phase  and  frequency  errors 
do  change,  however.  Specifically,  the  clock  phase  error  is 
outside  the  standard  deviation  boundary  approximately  60%  of 
the  time.  In  addition,  filter  divergence  (Ref.  5)  occurs  at 
290  seconds  for  the  clock  frequency  error  state.  This  is  in- 
dicative of  the  fact  that  the  filter  model  is  optimistic  in 


Figure  4.2-5  INS  Heading  Error  During  Flight 


90 


i ir-'^f-  ~'"ir  i -'n  rn  n 


r J 


i 


its  description  of  the  user  clock  error  behavior.  This  sit- 
uation can  be  remedied  by  increasing  qf  in  the  filter  model, 
but  this  will  result  in  non-optimum  performance  if  the  real 
world  clock  performs  in  the  manner  described  by  the  baseline 
model.  This  type  of  effect  will  be  illustrated  below  when  a 
worst  case  filter  model  is  developed  and  evaluated. 

Temperature  — The  effects  of  ambient  temperature 
variations  on  user  navigation  performance  are  shown  in 
Fig.  4.2-7.  A comparison  of  these  results  with  the  baseline 
performance  indicates  that,  for  the  assumed  temperature 
fluctuations  of  the  GPS  receiver,  the  temperature  sensitivity 
of  the  quartz  crystal  oscillator  has  minimum  impact  on  user 
navigation.  Thus,  this  error  source  need  not  be  inrluded  in 
user  Kalman  filter  design  considerations. 
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Vibration  - The  effects  of  vibration  of  the  re- 
ceiver crystal  oscillator  on  user  navigation  performance 
are  shown  in  Fig.  4.2-8.  A comparison  of  these  results 
with  the  baseline  performance  indicates  that,  for  the  as- 
sumed level  of  random  vibration  of  the  crystal,  the  vibra- 
tion sensitivity  of  the  quartz  crystal  oscillator  has  mini- 
mum impact  on  user  navigation.  It  should  be  noted,  however, 
that  a strong  sinusoidal  vibration  near  120  Hz  (see  Fig. 
4.1-7)  would  result  in  a large  frequency  error  and  could 
lead  to  a degradation  in  user  navigation  performance. 

Warm-up  — The  impact  of  quartz  crystal  oscillator 
warm-up  on  user  navigation  performance  is  considered  next. 

It  was  assumed  that  the  oscillator  had  been  turned-on  12 
minutes  prior  to  the  start  of  the  flight  such  that  the 
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initial  frequency  offset  (at  t=0)  was  3.0  mps  (3.0  mps  = 
c x Af(12  min)/f0,  where  c is  the  speed  of  light  and 
Af/f0  is  obtained  for  oscillator  "c"  from  Fig.  4.1-6). 

The  12  minute  period  from  oscillator  turn-on  until  the 
start  of  the  flight  was  chosen  to  be  compatible  with  a 
five  minute  coarse  gyrocompass  and  420  second  fine-align  of 
the  INS. 


Figure  4.2-9  illustrates  the  effects  of  oscillator 
wartn-up  when  it  is  not  modeled  in  the  user  filter.  The 
frequency  offset  present  during  warm-up  appears  directly 
as  the  predominant  contributor  to  frequency  estimation 
error.  This  large  frequency  offset  results  in  an  increase 
in  velocity  error  since  the  filter  cannot  properly  allo- 
cate the  large  pseudo-range-rate  measurement  residual. 

The  clock  phase  error  also  increases  due  to  the  warm-up 
induced  frequency  offset  and  consequently,  user  position 
errors  are  increased. 

To  compensate  for  the  impact  of  oscillator  warm-up, 
a model  for  the  warm-up  was  included  in  the  user  filter. 

The  model  is  identical  to  the  warm-up  model  presented  in 
Section  4.1.3.  Figure  4.2-10  illustrates  the  resultant  im- 
provement in  user  navigation  performance.  User  position, 
horizontal  velocity  and  clock  phase  errors  are  comparable 
to  their  baseline  values  Vertical  velocity  and  clock 
frequency  errors  exhibit  a distinct  improvement  from  their 
values  in  the  unmodeled  case,  but  they  are  somewhat  larger 
than  the  baseline  performance  for  the  first  400  seconds  of 
flight.  Beyond  this  point,  these  errors  compare  favorably 
to  their  baseline  values. 

The  use  of  warm-up  initial  conditions  at  12  minutes 
after  oscillator  turn-on  is  arbitrary;  however,  it  is  con- 
sistent with  the  time  frame  imposed  by  the  calibration 
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of  the  INS.  If  the  oscillator  were  turned-on  at  a point 
in  time  closer  to  the  start  of  the  flight,  the  frequency 
offset  would  be  larger  and  consequently  user  navigation 
performance  would  be  degraded.  Conversely,  if  a longer 
warm-up  period  were  allowed,  user  navigation  errors  would 
tend  to  be  closer  to  their  baseline  values. 

Acceleration  — The  effects  of  acceleration  of  the 
user  clock  on  navigation  performance  are  shown  in  Fig.  4.2-11. 
The  frequency  error  arising  from  the  g-sensit ivity  of  the 
quartz  oscillator  is  readily  apparent  as  the  airborne  user 
traverses  the  trajectory  of  Fig.  4.2-1.  The  impact  of  this 
frequency  error  is  manifest  in  the  horizontal  and  vertical 
velocity  errors.  When  the  user  resumes  straight  and  level 
flight  following  the  third  turn,  the  frequency  error  is 
reduced  through  processing  of  pseudo-range-rate  measurements. 

Since  the  filter  model  does  not  contain  any  information  about 
the  acceleration  sensitivity,  the  filter  gains  are  small  (at 
this  point)  and  do  not  weight  the  pseudo-range-rate  measure- 
ments heavily  enough  to  remove  the  residual  frequency  error 
entirely.  One  consequence  of  this  residual  frequency  error 
is  the  large  residual  velocity  error  which  remains  after 
straight  and  level  flight  is  resumed. 

The  clock  phase  error  grows  (to  a maximum  value  of 
200  m)  as  the  user  proceeds  through  the  turns  and  begins 
to  recover  during  the  straight  and  level  portion  of  the 
flight.  However,  since  the  filter  gains  are  small  (the 
filter  has  no  information  about  the  large  clock  phase  er- 
ror), the  phase  error  is  not  reduced  rapidly.  The  large 
horizontal  velocity  error  (sil  mps)  drives  the  horizontal 
position  error  to  60  to  70  m.  The  altitude  error  (44  m at 
t = 500  sec)  results  from  the  misallocation  (by  the  filter) 
of  the  large  pseudo-range  measurement  residual  which  is 
due  to  the  clock  phase  error  (these  errors  are  highly  cor- 
related) . 
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Figure  4.2-11 


Airborne  User  Errors  — Acceleration 
(Not  Modeled  in  User  Filter) 


This  type  of  performance  is  clearly  not  desirable. 

One  method  which  will  afford  a reduction  in  user  navigation 

errors  is  to  add  white  noise  to  the  user  filter  to  increase 

the  clock  frequency  error  variance.  Adding  a white  noise 

3 2 

with  a spectral  level  of  1.83  x 10  (mps)  /sec  [discrete 
q = 1.83  x 10  (mps)  ],  driving  the  frequency  error  state 
in  the  user  filter, results  in  the  performance  shown  in 
Fig.  4.2-12.  The  horizontal  position  and  velocity  errors 
are  reduced  to  approximately  their  baseline  values.  The 
altitude  and  vertical  velocity  errors  exhibit  dramatic 
improvement,  but  remain  larger  than  in  the  baseline 
case.  (The  altitude  error  remains  well  outside  the  stan- 
dard deviation  boundary  from  280  seconds  until  the  end  of 
the  flight . ) 

A significant  reduction  in  user  clock  phase  error 
is  observed,  with  a peak  excursion  of  11  m evident  during 
the  third  turn.  The  clock  phase  estimation  error  exceeds 
the  filter  standard  deviation  boundary,  from  the  beginning 
of  the  second  turn  until  the  third  turn  is  completed,  as  a 
result  of  the  large  spikes  in  the  clock  frequency  error  at 
250,  260  and  300  seconds  (viz. , the  beginning  of  the  se- 
cond and  third  turns).  A somewhat  weaker  spike  in  the 
frequency  error  is  evident  at  the  beginning  of  the  first 
turn  (i.e.,  210  seconds).  These  spikes  are  due  to  the  un- 
modeled acceleration-dependent  sensitivity  of  the  quartz 
oscillator;  however,  since  the  frequency  error  state  is 
being  driven  by  white  noise,  the  filter  gains  are  large 
enough  to  effect  a reduction  in  this  frequency  error  at  the 
next  pseudo-range-rate  measurement.  Thus,  the  frequency 
error  does  not  build-up  as  it  did  in  the  unmodeled  case 
(see  Fig.  4.2-11).  This  approach  does  have  disadvantages 
in  that  the  clock  phase  and  frequency  estimation  errors  are 
much  "noisier"  (i.e.,  exhibit  larger  rms  errors)  in  the 
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nonaccelerating  portion  of  the  flight  than  in  either  the 
unmodeled  or  baseline  cases. 

The  value  selected  for  the  spectral  level  of  the 
white  noise  was  chosen  to  give  acceptable  performance  and 
should  not  be  considered  as  optimal.  Indeed,  better  user 
performance  may  be  obtained  with  careful  "tuning"  of  the 
filter  which  was  not  attempted  in  this  study. 

Shock  — The  effects  on  user  navigation  performance 

of  a shock  to  the  receiver  crystal  oscillator  are  shown  in 

Fig.  4.2-13.  The  sharp  jump  in  frequency  error  (0.3  mps  = 
-9-9 

c x 10  , where  1 x 10  is  the  fractional  frequency  offset 

due  to  the  shock  at  t * 350  sec)  is  quite  evident  as  is  the 

rapid  growth  in  clock  phase  error  due  to  this  frequency  shift. 
These  errors  are  subsequently  reduced  through  processing  of 
pseudo-range  and  pseudo-range-rate  measurements;  however,  the 
recovery  is  not  rapid  since  the  filter  gains  are  small  during 
this  latter  portion  of  the  flight,  and  the  position  and  velo- 
city errors  appear  after  to  diverge  after  t = 350  seconds. 

Since  this  type  of  error  cannot  be  described  by  a 

Gauss-Markov  random  process  model,  the  only  method  of  im- 

proving the  filter's  response  to  an  unanticipated  shock  is 
to  increase  the  variance  of  the  frequency  error  state  (as 
was  done  to  handle  the  acceleration  sensitivity  of  the 
oscillator).  Including  a white  noise  [with  spectral  level 
q = 5.625  x 10-^  (mps )^/sec]  driving  the  frequency  error 
state  in  the  user  filter  results  in  the  performance  shown 
in  Fig.  4.2-14.  The  horizontal  position  and  velocity  er- 
rors are  reduced  to  approximately  their  baseline  values. 

The  altitude  and  vertical  velocity  errors  exhibit  some 
improvement  over  those  in  Fig.  4.2-13,  but  remain  somewhat 
larger  than  in  the  baseline  case.  (The  altitude  error 


102 


x - ESTIMATION  ERROR 


0 - STANDARD  DEVIATION  BOUNDARY 


OCCURRENCE 
OF  SHOCK 


TIME  <mc) 


TIME  (svc) 


Figure  4.2-13  Airborne  User  Errors  — Shock 

(Not  Modeled  in  User  Filter) 


remains  well  outside  the  standard  deviation  boundary  from 
350  seconds  until  the  end  of  the  flight.) 

A significant  reduction  in  user  clock  phase  error 
is  observed,  with  the  estimation  error  well  bounded  by  the 
filter  standard  deviation  boundary.  Similarly,  the  fre- 
quency estimation  error  is  within  the  boundary  at  all  but 
a few  points.  The  shock  has  minimal  effect  on  the  fre- 
quency error  since  the  larger  filter  gains  (due  to  the 
increased  frequency  error  variance)  reduce  the  error  sig- 
nificantly at  the  next  pseudo-range-rate  measurement. 

Here  again,  the  disadvantages  of  this  approach  are  evi- 
dent as  the  clock  phase  and  frequency  estimation  errors 
are  much  "noisier"  (i.e.,  exhibit  larger  rms  values) 
throughout  the  flight  than  in  either  the  unmodeled  or  base- 
line cases. 

The  value  selected  for  the  spectral  level  of  the  white 
noise  was  chosen  to  give  acceptable  performance  and  should  not 
be  considered  as  optimal.  As  with  the  previous  cases,  better 
user  performance  may  be  obtained  with  careful  "tuning"  of  the 
filter  which  was  not  attempted  in  this  study. 

Modified  User  Filter  — In  light  of  the  results 
obtained  in  the  previous  paragraphs,  it  is  possible  to  re- 
formulate the  user  Kalman  filter  such  that  it  "acknowledges" 
the  simultaneous  presence  of  all  the  environmentally-induced 
error  sources.  In  fact,  this  filter  model  has  been  developed 
in  parallel  with  the  evaluation  of  the  environmental  error 
sources.  The  error  sources  which  lead  to  a serious  degrada- 
tion in  user  navigation  performance  (if  they  are  not  modeled 
in  the  user  filter)  are  warm-up,  acceleration  and  shock.  In 
the  discussion  of  warm-up,  it  was  shown  that  inclusion  of  a 
model  describing  oscillator  warm-up  (decaying  exponential) 
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in  the  user  filter  resulted  in  a substantial  improvement  in 
user  navigation  performance.  Similarly,  in  the  discussions 
of  acceleration  and  shock  effects,  it  was  determined  that 
the  addition  of  white  noise  driving  the  clock  frequency  er- 
ror state  in  the  user  filter  provided  acceptable  user  per- 
formance in  the  presence  of  acceleration  and  shock  of  the 
receiver  quartz  oscillator. 

Modifying  the  user  filter  to  include  the  model  of 

oscillator  warm-up  [decaying  exponential  with  initial  vari- 

2 

ance  of  (3.05  mps)  ] and  the  additional  white  noise  driving 

-3  2 

the  frequency  error  [spectral  level,  q = 1.83x10  (mps)  / 
sec)],  should  provide  adequate  user  navigation  performance 
under  all  environmental  conditions.  This  performance  is 
indeed  realized  as  shown  in  Fig.  4.2-15.  For  this  run,  all 
of  the  environmental  error  sources  discussed  individually 
in  the  preceding  paragraphs  were  included  in  the  "real 
world"  clock  model.  The  horizontal  position  and  velocity 
errors  are  indistinguishable  from  their  values  in  the  base- 
line run  (in  which  it  was  assumed  that  no  environmentally 
induced  errors  existed).  The  increase  in  altitude  error 
during  the  second  turn  is  a result  of  the  unmodeled  phase 
error  growth  in  the  turn.  Due  to  the  high  correlation  be- 
tween altitude  and  clock  phase  errors,  the  filter  cannot 
distinguish  between  the  large  phase  error  and  the  altitude 
error,  thus  causing  an  increase  in  the  altitude  error  at 
the  pseudo- range  update.  Following  the  turn,  the  altitude 
error  settles  out  at  4 to  5 m.  A similar  situation  arises 
between  vertical  velocity  and  clock  frequency  errors.  The 
vertical  velocity  error  increases  at  the  pseudo-range-rate 
measurement  as  a result  of  the  jumps  in  frequency  error  but 
is  reduced  to  about  0.05  mps  during  straight  and  level  flight. 
The  altitude  and  vertical  velocity  errors  are  much  smaller 
with  the  modified  filter  than  they  would  be  with  the  original 
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filter  (compare  with  the  errors  resulting  from  acceleration 
of  the  oscillator  prior  to  the  filter  modification,  shown  in 
Fig.  4.2-11).  Although  filter  design  is  not  the  objective 
of  this  study,  it  should  be  possible  to  further  reduce  these 
errors  by  including  process  noise  driving  the  altitude  and 
vertical  velocity  states  in  the  filter.  This  alternative 
was  not  pursued  here. 

The  user  clock  phase  and  frequency  errors  are  well 
behaved  and  remain  below  11  m and  0.6  mps,  respectively, 
throughout  the  flight.  These  errors  exhibit  an  rms  level 
larger  than  the  baseline  case  (Fig.  4.2-4),  but  the  error 
growth  which  would  be  present  due  to  acceleration  and  shock 
(Figs.  4.2-11  and  4.2-13)  does  not  appear. 

User  navigation  performance  in  the  presence  of  a 
worst  case  environment  is  shown  in  Fig.  4.2-15  with  the  user 
filter  modified  to  "acknowledge"  the  significant  error  sources. 
These  changes  to  the  user  filter  will  have  a minor  adverse 
effect  on  user  navigation  performance  when  the  environmental 
errors  are  not  present  (as  represented  by  the  baseline  case). 
This  is  illustrated  in  Fig.  4.2-16  in  which  the  baseline  run 
(see  Fig.  4.2-4)  was  repeated  with  the  user  filter  modified 
to  take  into  account  the  possibility  of  environmentally  in- 
duced errors.  The  apparent  degradation  in  vertical  channel 
performance  is  associated  with  the  particular  random  sequence 
generated  in  this  Monte  Carlo  run.  Throughout  the  flight, 
the  clock  phase  and  frequency  estimation  errors  exhibit 
larger  excursions  than  in  the  previous  baseline  run  due  to 
the  increased  uncertainty  attributed  to  these  states  in  the 
user  filter. 
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Figure  4.2-16  Airborne  User  Errors  — Modified  User  Filter 

(No  Environmental  Error  Sources) 


4.2.4  Airborne  User  Performance  Summary 


The  results  of  the  airborne  user  performance  analysis 
may  be  summarized  as  follows: 

• Acceleration,  shock  and  warm-up  of  the 
GPS  receiver  quartz  crystal  oscillator 
are  the  dominant  environmental  error 
sources.  The  temperature  and  vibration 
sensitivities  of  the  oscillator  and  the 
effect  of  a degraded  quality  oscillator 
have  minimal  impact  on  user  navigation 
errors. 

• The  environmental  error  sources  have  direct 
impact  upon  the  oscillator  phase  and  fre- 
quency errors.  Due  to  the  high  correlation 
between  clock  errors  and  vertical  channel 
errors,  large  unmodeled  phase  and  fre- 
quency errors  tend  to  increase  altitude  and 
vertical  velocity  errors  at  filter  updates; 
i.e.,  the  filter  cannot  properly  attribute 
the  large  pseudo-range  and  pseudo-range- 
rate  measurement  residuals  to  the  clock 
states  alone  and  thus  "divides"  the  large 
residual  between  clock  and  vertical  chan- 
nel errors.  The  horizontal  position  and 
velocity  errors  are  insensitive  to  all  but 
extreme  environmentally  induced  errors 
(e.g.,  acceleration). 

• A simple  modification  of  the  user  fil- 
ter to  model  oscillator  warm-up  and 
the  effects  of  acceleration  and  shock 
substantially  improves  user  navigation 
performance , if  these  error  sources  are 
present.  However,  in  the  absence  of 
these  errors,  the  modified  filter  in- 
troduces some  degradation  in  user 
navigation  performance.  Careful  "tun- 
ing" of  the  user  filter  or  an  adaptive 
filtering  mechanization  could  reduce 
this  sensitivity. 


✓ 


110 


5. 


SUMMARY  AND  CONCLUSIONS 


Performance  studies  of  Army  users  of  the  NAVSTAR 
Global  Positioning  System  (GPS)  in  a battlefield  environ- 
ment are  presented  in  this  report.  These  studies  represent 
a continuation  of  an  earlier  Army  user  performance  analysis 
conducted  by  TASC  and  reported  in  Ref.  1.  The  problems 
analyzed  in  this  report  include  an  assessment  of  the  impact 
on  the  Army  manpack  GPS  user  navigation  performance  of  the 
use  of  directional  receiver  antennas.  These  antennas  can 
be  employed  to  increase  the  GPS  navigation  signal  strength 
while  reducing  the  users'  vulnerability  to  jamming.  Methods 
of  countering  jamming  threats  are  considered  and  evaluated. 

An  assessment  is  also  made  of  the  environmentally-induced 
error  effects  of  the  user  receiver  crystal  oscillator  (clock) 
on  navigation  performance. 


5.1  STUDY  HIGHLIGHTS 

User  Performance  with  Directional  Antennas  - The 
extent  to  which  directional  antennas  restrict  user  satel- 
lite visibility  and  the  degradation  in  navigation  perform- 
ance which  results  with  their  use  were  investigated  in 
this  report.  The  navigation  performance  of  a stationary 
user  with  a 120  deg  beamwidth  antenna  is  illustrated  in 
Fig.  5.1-la.  After  the  initial  transient  period  (during 
which  the  initial  set  of  measurements  to  the  four  visible 
satellites  are  being  made),  the  uncertainty  of  the  posi- 
tion estimates  of  this  restricted  visibility  user  is 
greater  than  if  satellite  visibility  was  unrestricted, 
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>)  Widebeam  Antenna  b)  Conventional  Antenna 

Figure  5.1-1  Horizontal  Position  Error  with 

Wide  Beam  and  Omni-Directional 
Antennas 

as  illustrated  in  Fig.  5.1-lb.  This  is  attributed  to  a 
poor  GDOP  condition  for  the  four  satellites  available  with 
the  120°  beamwidth  antenna. 

A stationary  manpack  user  could  significantly  in- 
crease the  strength  of  navigation  signals  while  reducing 
his  vulnerability  to  jamming  by  employing  a very  narrow 
beam  directional  antenna  which  restricts  visibility  to  one 
satellite.  This  antenna  must  be  physically  pointed  by  the 
user  (with  directions  provided  by  his  GPS  navigation  com- 
puter) to  each  satellite  for  measurements.  This  antenna/ 
deployment  concept  was  evaluated  in  this  study,  indicating 
that  performance  comparable  to  that  shown  in  Fig.  5.1-1  is 
achievable  after  the  user  has  sequenced  through  five  satel- 
lites . 


As  a byproduct  of  the  directional  antenna  study, 
it  was  discovered  that  significant  improvement  in  manpack 
user  performance  could  be  achieved  if  one  additional  state 
variable  ("pseudo-state")  per  satellite  were  included  in 
the  user  filter  software  to  account  for  satellite  clock 
and  ephemeris  errors.  Comparisons  of  this  filter  design 
with  alternative  designs  are  included  in  this  study. 

112 


Antenna  Array  Systems  — Directional  antennas  can  be 
constructed  from  an  array  of  discrete  antennas  or  elements. 
These  antenna  arrays  can  be  designed  to  provide  both  spatial 
and  frequency  filtering  to  enhance  the  desired  signal  while 
reducing  directional  interference  or  jamming.  This  study 
indicates  significant  potential  reductions  in  interference 
are  possible  with  array  antennas  for  Army  GPS  manpack  users. 
Futher  investigations  of  methods  of  implementing  the  opti- 
mum array  are  required  to  determine  practically  realizable 
performance  and  cost. 

User  Vulnerability  — The  range  from  a continuous 
wave  jammer  within  which  a ground-based  user  could  be  denied 
access  to  the  GPS  signal  is  illustrated  in  Fig.  5.1-2  as  a 
function  of  jammer  type  and  power,  and  user  receiver  operat- 
ing mode.  The  user  receiver  is  most  vulnerable  to  jamming 
during  initial  acquisition  of  the  satellite  navigation 


Critical  Jammer  Range 
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Figure  5.1-2 


signals,  because  the  initial  large  frequency  uncertainty  of 
these  signals  require  a large  receiver  bandwidth.  If  the 
jamming  commences  after  the  receiver  has  acquired  the  navi- 
gation signals,  it  can  operate  in  the  tracking  mode  with  a 
significant  reduction  in  critical  jammer  range,  as  illustrated. 

Five  techniques  of  countering  a jammer  threat  that 
are  considered  in  this  report  include  notch  filters,  reduced 
receiver  bandwidth  (navigation  mode  only),  directional  an- 
tennas, null  steering  antennas  and  RFI  screens.  It  is  shown 
in  this  report  that  the  null  steering  antenna  and  notch  fil- 
ter are  most  effective  against  an  airborne  jammer  while  a 
combination  of  the  reduced  receiver  bandwidth  and  RFI  screen 
are  effective  against  a ground  jammer. 

Receiver  Clock  Error  Analysis  — An  examination  of 
the  potential  impact  of  environmentally-induced  errors  in 
the  GPS  receiver  quartz  crystal  oscillator  on  user  naviga- 
tion performance  was  conducted.  Errors  considered  in  this 
study  include  those  due  to  vibration,  shock,  warm-up,  tem- 
perature variations  and  acceleration  of  the  clock  platform. 

A comprehensive  model  of  the  oscillator  error  behavior, 
under  assumed  operating  conditions,  was  developed  and  em- 
ployed to  obtain  the  performance  projections  for  an  air- 
borne (helicopter)  user  employing  a four  channel  GPS  re- 
ceiver and  an  inertial  navigation  system. 

The  environmental  error  sources  which  proved  to 
be  dominant  are  acceleration,  shock  and  oscillator  warm- 
up . The  oscillator  temperature  and  vibration  sensitivi- 
ties and  the  effects  of  an  unknown  increase  in  the  oscil- 
lator flicker  noise  level  have  only  a minimal  effect,  on 
user  navigation  performance.  The  navigation  performance 
of  the  airborne  user,  illustrating  the  impact  of  oscil- 
lator acceleration  sensitivity  (not  modeled  in  the  user 
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Kalman  filter)  is  shown  in  Fig.  5.1-3a.  However,  the  navi- 
gation errors  arising  from  the  dominant  error  sources  can 
be  reduced  via  simple  modifications  to  the  user  Kalman  fil- 
ter. This  is  illustrated  in  Fig.  5.1-3b  for  the  acceleration- 
sensitive  case. 


TIME  (lec) 


a)  Not  Modeled  in  User 
Kalman  Filter 


b)  "Acknowledged"  in  User 
Kalman  Filter 


Figure  5.1-3  Impact  of  Quartz  Oscillator  Acceleration 

Sensitivity  on  Airborne  User  Navigation 
Performance  (see  Section  4.2.3) 


5 . 2 RECOMMENDATIONS 

This  study  has  indicated  the  need  for  additional 
analytical  effort  in  the  following  areas: 

User  Filter  Design  — Additional  analy- 
sis of  user  filters  should  be  performed 
with  particular  emphasis  on  elimination 
or  compensation  of  modeling  deficiencies. 
This  report  has  indicated  that  the  addi- 
tion of  pseudo-states  is  an  effective 
method  of  compensating  for  unmodeled 
satellite  ephemeris  and  clock  errors. 

This  technique  requires  further  analysis 
to  determine  its  full  potential  and  ad- 
vantages. A method  for  reducing  the 
impact  of  unmodeled  user  clock  errors 
is  the  inclusion  of  white  process  noise 
in  the  filter,  driving  the  user  clock 
frequency  errors.  Optimum  determination 
of  the  white  noise  spectral  level  (i.e., 
"tuning"  of  the  filter)  would  require 
further  study. 
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• Interference  Reduction  Array  Implemen- 
tation — This  report  has  shown  the  po- 
tential advantage  of  using  interference 
reduction  (or  null  steering)  arrays  to 
reduce  GPS  user  vulnerability  to  jamming. 
Methods  of  implementing  this  technique 
should  be  investigated  to  determine  per- 
formance, cost  and  computational  require- 
ments. A promising  alternative  to  this 
technique,  not  investigated  in  this  study, 
is  the  use  of  a tracking  filter  (e.g., 
with  phase  lock  loops). 

• User  Vulnerability/Receiver  Design  — 

This  report  has  illustrated  that  user 
vulnerability  to  jamming  is  highly 
dependent  upon  the  threshold  jammer- 
to-signal  power  ratio  of  the  user  re- 
ceiver. Further  analysis  of  receiver 
designs  are  required  to  determine  poten- 
tial improvements  in  this  threshold  for 
proposed  receiver  designs.  Simulation 
of  the  user  receiver  with  jamming  in  a 
typical  manpack  scenario  would  provide 
confirmation  of  user  vulnerability.  In 
addition,  this  analysis  could  determine 
the  feasibility  of  reducing  the  manpack 
user  receiver  carrier  loop  bandwidth  to 
improve  his  jamming  resistance. 


116 


APPENDIX  A 


DIRECTIONAL  ANTENNA  UTILIZATION  CONSIDERATIONS 


The  effectiveness  of  directional  antennas  depends 
upon  the  amount  of  gain  they  can  provide,  as  well  as  their 
beamwidth.  Increased  gain  and  decreased  beamwidth  provides 
greater  immunity  to  signal  interference  such  as  jamming,  mul- 
tipath, rotor  blade  modulation  and  foliage  attenuation.  To 
evaluate  the  results  of  Chapter  2,  the  achievable  gain  as  a 
function  of  beamwidth  must  be  determined.  That  is  the  sub- 
ject of  this  appendix.  In  particular,  equations  relating 
antenna  gain,  beamwidth  and  size  are  developed  in  Section 
A . 1 . These  equations  are  applied  to  a parabolic-type  an- 
tenna in  Section  A. 2.  This  type  of  antenna  is  often  used 
in  high  frequency  communication  systems  where  medium  to  high 
gain  is  required.  Thus  the  results  of  Chapter  2 can  be  given 
some  physical  interpretation  through  the  equations  developed 
in  this  section. 

The  definitions  in  this  appendix  are  applicable  to 
any  antenna,  although  some  of  the  quantities  (in  particular, 
effective  cross-sectional  area)  are  difficult  to  calculate 
for  most  antenna  types.  References  12  to  15  contain  addi- 
tional details  on  antennas.  In  particular,  Ref.  14  con- 
tains a good  discussion  of  parabolic  antennas. 


A . 1 ANTENNA  CHARACTERISTICS 

A measure  of  an  antenna's  ability  to  transmit  or  re- 
ceive radio  frequency  power  is  described  by  its  antenna  power 
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pattern  w(0,4>),  which  is  defined  as  the  power  transmitted  in 
(or  received  from)  a direction  (0,4>)  per  unit  solid  angle 
(Refs.  12,  13  and  15).  Figure  A. 1-1  defines  the  spatial 
geometry  of  the  antenna  used  in  this  report.  The  conceptu- 
ally simplest  antenna  is  the  ideal  isotropic  antenna  whose 
antenna  pattern  is  constant  (unity)  in  all  directions.  The 
isotropic  antenna  is  often  used  as  a reference  with  which 
to  compare  other  antennas. 


* -I 5579 


ANTENNA 

. direction 
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ANTENNA  CENTER 


Figure  A. 1-1  Antenna  Spatial  Geometry 


The  total  power,  PT>  (transmitted  or  received)  assoc- 
iated with  an  antenna  is  the  surface  integral  of  the  antenna 
pattern  over  the  surrounding  sphere  of  unit  radius,  i.e., 


= J"  w(  0 , <)) ) dfi 
sphere 


(A. 1-1) 
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The  power  associated  with  an  ideal  isotropic  antenna  is  thus 
just  the  surface  area  of  a unit  sphere  or  4tt  . 


The  antenna  pattern  is  often  normalized  to  form  the 
antenna  gain  function.  The  normalization  term  is  the  power 
per  unit  solid  angle  of  an  isotropic  antenna  with  the  same 
total  power  PT.  Thus,  the  antenna  gain  function  is 


gC  PT/4n 


(A. 1-2) 


The  gain  or  directivity  of  an  antenna  is  taken  as 
the  maximum  value  of  its  gain  function,  i.e., 


G = max  g(8,4>) 

9,4) 


(A. 1-3) 


or  in  terms  of  decibel  rating 


G(db ) = 10  log10  G 


(A. 1-4) 


The  actual  gain  of  an  antenna  is  lower  than  this  theo- 
retical gain  due  to  coupling  and  feed  losses.  The  relation 


between  the  theoretical  antenna  gain  G and  actual  gain  G is 

El 


the  antenna  gain  efficiency  factor  n , i.e. , 


= ng  G 


(A. 1-5) 


The  antenna  gain  efficiency  factor,  n , is  typically  50%,  thus 

& 


the  actual  antenna  gain  is  3 dB  less  than  the  theoretical  gain. 
However,  with  good  design,  this  efficiency  factor  may  be  only 
1 to  2 dB. 


The  antenna  solid  beamwidth  is  defined  as  the  solid 


angle,  fi,  through  which  all  of  the  radiated  power  would  pass 
if  the  gain  in  this  angle  were  constant  and  equal  to  its 
maximum  value.  Thus,  the  solid  beamwidth  is  given  by 
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The  effective  cross-sectional  area,  a(9,4>),  of  an 
antenna  is  the  total  power  collected  at  the  antenna  termin- 
als divided  by  the  incident  power  density  in  the  direction 
(9, 4> ) • This  is  a measure  of  the  amount  of  power  that  a re- 
ceiving antenna  can  absorb  from  a passing  wave.  The  physical 
antenna  area  A is  related  to  the  maximum  cress-sectional  an- 

cL 

tenna  area  A = max  a(9,<j>)  by  the  antenna  aperture  efficiency 

factor  n , i . e . , 
a 

A = na  Aa  (A. 1-7) 

Thus,  the  maximum  cross-sectional  area  of  an  antenna  is  al- 
ways less  than  the  physical  area  of  the  antenna  due  to  the 
aperture  efficiency  factor. 

The  maximum  cross-sectional  area  of  an  antenna  can 
be  expressed  in  terms  of  antenna  gain  and  signal  wavelength 
by  (Ref.  13) 


where  X is  the  radio  frequency  wavelength  of  the  signal. 

Equations  (A. 1-6)  and  (A. 1-8)  illustrate  the  relationship 

between  antenna  gain,  beamwidth  and  size.  Antenna  gain  is 

2 

directly  proportional  to  (frequency)  and  antenna  size.  An- 


tenna beamwidth  is  inversely  proportional  to  frequency  and 
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A.  2 PARABOLIC  ANTENNA  PARAMETERS 

The  gain  of  an  antenna  as  a function  of  beamwidth 
and  antenna  efficiency  factor  is  [from  Eqs.  (A. 1-5)  and 
(A. 1-6)] 


The  solid  angle  beamwidth  ft  is  the  surface  area  of  that 
segment  of  a unit  sphere  which  intersects  a solid  cone  with 
conical  angle  6 and  apex  at  the  center  of  the  antenna,  as 
i: lustrated  in  Fig.  A. 2-1.  From  spherical  calculus,  the 
relation  between  ft  and  conical  beamwidth  3 is 

ft  = 2tt [l  - cos  (3/2)]  (A. 2-2) 
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Figure  A. 2-1  Definition  of  Beamwidth  of 

a Parabolic  Antenna 

Combining  Eqs.  (A. 2-1)  and  (A. 2-2),  the  antenna  gain  G as 

£1 

a function  of  (conical)  beamwidth  3 is 


; - - - rir  Tin 


fl 


» 
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'1  - cos  (6/2)' 


(A. 2-3) 


This  relation  is  presented  in  Fig.  A. 2-2  for  two  gain  effi- 
ciency factors  expressed  in  decibels.  (The  -1  dB  and  -3  dB 
efficiency  factors  correspond  to  80%  and  50%  efficiencies, 
respectively . ) 


<J0  • ANTENNA  POLARIZATION 
AND  COUPLING  LOSSES 
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*o  40  to  to  m no  mo  ito  ito 
IEAMWI0TH  (DEG) 


Figure  A. 2-2 


Parabolic  Antenna  Gain  vs  Beamwidth 


The  gain  of  an  ideal  parabolic  antenna  as  a function 
of  size  is  obtained  from  Eqs.  (A. 1-5),  (A. 1-7)  and  (A. 1-8)  as 


II  A 

= n n — o A 
a g a ^2  a 


(\.2-4) 


The  physical  area  a parabolic  antenna  presents  to  a radio 
frequency  wave  is  the  area  of  the  front  of  the  antenna  or 


(A. 2-5) 
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where  d is  the  diameter  of  the  antenna  "dish."  Thus,  the 
gain  of  a parabolic  antenna,  in  terms  of  its  diameter,  is 


1 0.  - ng  na  (^)2  (A., -6 

where 

f is  signa)  frequency  in  hertz 
c is  propagation  constant  = 3 x 10®  meters/sec 
d is  diameter  of  parabolic  antenna  in  meters 

This  relation  is  presented  in  Fig.  A. 2-3  for  efficiency  fac- 
tors of  -1  dB  and  the  two  GPS  navigation  signal  frequencies. 


Figure  A. 2-3  Parabolic  Antenna  Gain  vs  Diameter 

The  beamwidth  of  a parabolic  antenna  as  a function 
of  its  size  or  diameter  is  obtained  from  Eqs.  (A. 1-6),  (A. 1-7), 
(A. 1-8),  (A. 2-2)  and  (A. 2-5)  as 

8 = 2 cos"1  fl -1  (A. 2-7) 


where 

n is  antenna  aperture  efficiency  factor 
a 

d is  antenna  size 
X is  radio  frequency  wavelength 

This  relationship  is  presented  in  Fig.  A. 2-4  for  an  aperture 
efficiency  factor  of  -1  dB  and  the  two  GPS  frequencies. 
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Figure  A.  2-4 


DIAMETER  (METERS) 

Parabolic  Antenna  Beamwidth 
vs  Diameter 


This  last  figure  illustrates  the  different  beamwidths 
that  occur  at  the  two  frequencies  for  a constant  antenna  size 
For  example,  if  a beamwidth  of  140  deg  at  1.6  GHz  is  desired, 
the  antenna  diameter  should  be  0.12  meters  (5  inches).  At 
1.2  GHz,  the  beamwidth  is  180  deg  for  this  size  antenna  and 
the  gain  is  2 dB  less  than  at  1.6  GHz.  Thus,  some  compromise 
in  performance  might  be  required  if  a single  antenna  is  em- 
ployed for  both  frequencies. 
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A. 3 PARABOLIC  ANTENNAS  FOR  GPS  USERS 


A parabolic  antenna  could  be  employed  by  a GPS  user 
in  order  to  increase  the  signal  strength  at  the  receiver  and 
to  reduce  interference  such  as  jamming.  The  data  presented 
in  Section  A. 2 can  be  used  to  determine  the  amount  of  signal 
gain  and  directivity  a user  can  obtain  using  parabolic  anten- 
nas under  various  operating  conditions. 

It  has  been  observed  in  Section  2.1  that  the  opera- 
tional Phase  III  GPS  will  provide  at  least  four  satellite 
visibility  to  a user  employing  a vertical  boresight  antenna 
with  a beamwidth  of  140  deg.  From  Fig.  A. 2-2,  an  antenna 
with  140  deg  beamwidth  has  a 1.8  to  3.8  dB  gain,  depending 
upon  the  efficiency  factor.  From  Fig.  A. 2-4,  this  beam- 
width  corresponds  to  an  antenna  size  of  0.12  meters  at 
1.6  GHz.  However,  at  1.2  GHz  the  beamwidth  is  180  deg  for 
this  antenna  size  and  the  gain  is  2 dB  less.  This  may  not 
be  a serious  problem  if  the  higher  frequency  is  used  for 
initial  acquisition  of  the  signals.  Acquisition  of  the  sig- 
nals at  1.2  GHz  is  then  relatively  easy  since  the  signals 
on  these  two  frequencies  are  synchronized.  For  120  deg 
beamwidth,  the  antenna  gain  is  3 to  5 dB. 

For  the  narrow  beam  antenna  of  Section  2.2,  if 
the  bandwidth  were  60  deg  then  a gain  of  8.7  to  10.7  dB 
would  be  provided  with  a diameter  of  0.26  m at  n = -1  dB. 

ct 

This  represents  a significant  gain  in  signal  strength  which 
may  be  desired  under  foliage  conditions.  However,  this 
antenna  must  be  physically  slewed  by  the  operator  to  each 
satellite  for  acquisition  and  pseudo-range  measurements, 
since  four  satellites  cannot  be  seen  simultaneously  by 
this  antenna. 

These  examples  illustrate  the  use  of  the  design  data 
developed  in  Section  A. 2 to  determine  user  performance  im- 
provements and  limitations  with  directional  antennas. 
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APPENDIX  B 


GPS  PHASE  III  CONFIGURATION  AND  PERFORMANCE 


Realistic  user  performance  projections  require  real- 
istic satellite  ephemerides  and  measurement  data.  This  data 
is  obtained  by  simulating  the  satellites  and  the  ground 
tracking  system  that  computes  the  satellite  ephemerides  which 
are  relayed  to  the  user.  The  user's  performance  is  thus  de- 
pendent upon  the  performance  of  the  tracking  system.  This 
appendix  presents  details  of  the  ground  tracking  system  sim- 
ulation that  was  performed  to  generate  the  data  required  in 
the  user  performance  projections  of  Sections  2.2  and  4.2. 
Further  details  of  the  simulation  program  are  given  in  Ref.  2. 

This  study  is  concerned  with  the  fully  operational 
Phase  III  satellite  configuration.  In  order  to  reduce  the 
computer  time  required  in  the  tracking  system  simulation,  a 
subset  of  the  Phase  III  satellite  configuration  was  used. 

The  satellites  selected  and  the  rationale  for  their  selec- 
tion are  discussed  in  Section  B.l.  The  error  models  that 
are  peculiar  to  this  study  (and  thus  not  in  Ref.  2)  are 
presented  in  Section  B.2.  Also,  the  initial  conditions 
used  in  the  tracking  filter  are  given.  Finally,  the  track- 
ing system  performance  results  are  summarized  in  Section  B.3. 
Included,  are  illustrations  of  typical  performance  for  each 
of  the  tracking  system  error  sources. 


B.l  SATELLITE  CONFIGURATION 

The  GPS  Phase  III  satellite  configuration  consists 
of  24  satellites  in  a 3 x 8 (three  planes,  eight  satellites 
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per  plane) constellation  with  each  satellite  in  a 12-hour  orbit 
inclined  63  deg.  The  orbital  elements  of  this  constellation  of 
the  initiation  of  the  tracking  process  are  summarized  in  Table 
B.l-1  with  the  satellite  geometry  shown  in  Fig.  B.l-1. 


TABLE  B.l-1 

GPS  PHASE  III  SATELLITE  POSITIONS 
AT  THE  START  OF  TRACKING 


SATELLITE 

NO. 

SEMI MAJOR 
AXIS 

ECCENTRICITY 

INCLINATION 

RICHT 
ASCENSION 
OF  THE 
ASCENDING 
NODE 

ARGUMENT 

OF 

PERIGEE 

MEAN 

ANOMALY 

a 

• 

1 

n 

k» 

U 

earth  radii 

dec 

dee 

deg 

deg 

1 

4.16449 

0 

63 

0 

0 

15 

2 

4.16449 

0 

63 

0 

0 

60 

3 

4.16449 

0 

63 

0 

0 

105 

4 

4.16449 

0 

63 

0 

0 

150 

S 

4.16449 

0 

63 

0 

0 

195 

6 

4.16449 

0 

63 

0 

0 

240 

7 

4.16449 

0 

63 

0 

0 

285 

8 

4.16449 

0 

63 

0 

0 

330 

9 

4.16449 

0 

63 

120 

0 

348 

10 

4.16449 

0 

63 

120 

0 

30 

11 

4.16449 

0 

63 

120 

0 

75 

12 

4.16449 

0 

63 

120 

0 

120 

13 

4.16449 

0 

63 

120 

0 

165 

14 

4.16449 

0 

63 

120 

0 

210 

15 

4.16449 

0 

63 

120 

0 

255 

16 

4.16449 

0 

63 

120 

0 

300 

17 

4.16449 

0 

63 

240 

0 

0 

18 

4.16449 

0 

63 

240 

0 

45 

19 

4.16449 

0 

63 

240 

0 

90 

20 

4.16449 

0 

63 

240 

0 

135 

21 

4.16449 

0 

63 

240 

0 

180 

22 

4 16449 

0 

63 

240 

0 

225 

23 

4. iS449 

0 

63 

240 

0 

270 

24 

4.16449 

0 

03 

240 

0 

315 

The  user  scenarios  in  this  report  are  relatively  short 
(less  than  one  hour)  and  thus,  only  a limited  set  of  satellites 
will  be  visible  to  the  user  during  this  period.  The  computer 
time  required  to  provide  a realistic  tracking  system  simulation 
is  dependent  upon  the  number  of  satellites  used.  This  time  can 
be  reduced  if  only  those  (four)  satellites  which  are  visible  and 
provide  the  best  GDOP  to  a user  (with  a conventional  antenna) 
during  the  scenario  time  period  are  included  in  the  tracking 
system  simulation . A wide  beam  (120  deg)  directional  antenna 
user  will  see  three  of  these  four  satellites,  plus  a fifth. 
Thus,  five  satellites  were  included  in  the  simulation. 
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The  satellites  chosen  to  include  in  the  simulation 
are  listed  in  Table  B.l-2  and  correspond  to  Satellite  Nos. 
10,  11,  12,  17,  and  18,  respectively,  in  Table  B.l-1.  These 


TABLE  B.l-2 

SATELLITE  ORBITAL  ELEMENTS  AT  THE 
START  OF  TRACKING 


SATELLITE 

NO. 

SEMI MAJOR 
AXIS 

ECCENTRICITY 

INCLINATION 

RIGHT 
ASCENSION 
OF  THE 
ASCENDING 
NODE 

ARGUMENT 

OF 

PERIGEE 

MEAN 

ANOMALY 

a 

e 

i 

B 

(l) 

M 

earth  radii 

deg 

deg 

deg 

deg 

1 

4.16449 

0 

63 

120 

0 

30 

2 

4.16449 

0 

63 

120 

0 

75 

3 

4.16449 

0 

63 

120 

0 

120 

4 

4.16449 

0 

63 

240 

0 

0 

5 

4.16449 

0 

63 

240 

0 

45 
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five  satellites  provide  a wide  range  of  satellite  visibility 
and  GDOP  conditions  to  users  located  at  latitude  +50  deg, 
longitude  -100  deg  as  illustrated  in  Fig.  B.l-2.  Ground 
tracks  of  these  five  satellite  subpoints  during  a three- 
hour  period  are  illustrated  in  Fig.  B.l-3.  The  satellite 
points  along  each  ground  track  at  one-hour  intervals  are 
depicted  by  hatch  marks. 
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Figure  B.l-2  User  GDOP  and  Visibility 

The  ground  tracking  system  was  simulated  for  49 
hours  of  simulation  time  at  which  point  the  user  scenarios 
were  initiated.  The  simulation  then  continued  for  three 
hours  or  to  52  hours  after  the  commencement  of  the  simula- 
tion. Thus,  the  user  scenarios  can  last  up  to  three  hours. 
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Figure  B.l-3  Satellite  Ground  Tracks 


B. 2 ERROR  MODELS 

The  five  satellites  are  tracked  by  an  "optimal" 
filter  in  that  it  is  assumed  that  all  unknown  quantities 
which  will  have  an  impact  on  the  orbital  estimation  process 
are  included  as  state  variables  in  the  filter.  Also  the 
statistics  associated  with  each  state  variable  are  assumed 
to  be  modeled  correctly. 

The  details  of  the  error  models  used  in  this  simu- 
lation are  given  in  Ref.  2 and  summarized  in  Ref.  1.  The 
only  significant  changes  to  these  error  models  were  a re- 
duction of  some  of  the  initial  estimation  errors  (especially 
a 2%  solar  radiation  force  uncertainty  and  a reduction 
in  gravitational  constant  uncertainty)  and  the  incorporation 
of  a new  model  for  the  satellite  and  monitor  station  clocks, 


y 
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The  new  clock  model  is  of  a high  performance  cesium 
beam  tube  time  standard  and  consists  of  three  states:  phase  error, 
frequency  error  and  a one-state  flicker  noise  modeled  as  a Markov 
process  with  a time  constant  of  approximately  one  day.  (Ref.  3 
and  Appendix  C).  This  clock  model,  illustrated  in  Fig.  B.2-1, 
is  based  on  data  contained  in  Ref.  4.  The  fractional  frequency 
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Figure  B.2-1  Cesium  Beam  Clock  Error  Model 

stability  plot  for  this  clock  model  is  shown  in  Fig.  B.2-2.  The 
error  states  for  the  model  are  defined  in  Table  B.2-1.  The 
spectral  density  (continuous)  matrix  associated  with  this  clock 
model  is  given  by  Eq.  (B.2-1).  The  parameters  k^,k2  and  g shown 
in  Fig.  B.2-1  have  values  0.1,0.03  and  lxlO-5,  respectively. 

The  satellites  are  tracked  by  four  ground  tracking 
stations  located  in  the  northern  hemisphere  in  United  States 
territory.  The  locations  of  the  tracking  stations  are  given 
in  Table  B.2-2  with  the  first  three  also  shown  in  Fig.  B.l-2. 

131 


R— 16596 


Figure  B.2-2  Fractional  Frequency  Stability  Data 

for  Cesium  Beam  Clock  Model 


TABLE  B.2-1 

CESIUM  BEAM  CLOCK  ERROR  STATE  DEFINITIONS 


SYMBOL 

ERROR  SOURCE 
DEFINITION 

UNITS 

5<j> 

Clock  Phase  Error 

m 

6f 

Clock  Frequency  Err^r 

mps 

xf 

Flicker  Noise 

mps 
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0 


5.95  x 10 


7.35  x 10 
0 

5.95  x 10 


-6 


-10 


1.11  x 10 


-15 


(B.2-1) 


■10 


5.35  x 10 
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TABLE  B.2-2 

GROUND  TRACKING  STATION  LOCATIONS 


STATION 

NO. 

NAME 

LATITUDE 

deg 

LONGITUDE 

deg 

ASSUMED 

ALTITUDE 

m 

1 

Vandenberg  AFB 

34.7 

-120.5 

0 

2 

Elmendorf  AFB 

61.3 

-149.5 

0 

3 

Hawaii 

20.0 

-158.0 

0 

4 

Guam 

13.0 

145.0 

0 

Uploading  of  the  satellites  with  current  ephemeris 
data  occurs  once  a day  when  the  individual  satellites  make 
their  closest  approach  to  the  master  tracking  station  at 
Vandenberg  AFB.  The  time  of  the  uploading  to  each  satellite 
is  given  in  Table  B.2-3. 


TABLE  B.2-3 

SATELLITE  DAILY  INFORMATION  UPLOADING  TIMES 


SATELLITE 

NO. 

TIME 

(min) 

1 

210 

2 

90 

3 

1395 

4 

105 

5 

1380 
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The  ionospheric  and  tropospheric  delay  errors  are 

assumed  to  be  accounted  for  perfectly  in  the  measurements, 

and  the  only  measurement  error  is  a normally-distributed 

white  measurement  noise  with  zero  mean  and  variance  of 
2 * 

(0.3m)  . Measurements  are  processed  at  15-minute  inter- 
vals for  all  satellites  in  view  of  the  ground  tracking  sta- 
tions, while  the  filter  states  are  propagated  between  these 
measurements  using  one-minute  integration  steps.  The  track- 
ing filter  was  initialized  with  the  filter  statistics  pre- 
sented in  Table  B.2-4. 


B . 3 PERFORMANCE  RESULTS 

Typical  results  of  the  satellite  tracking  system 
for  the  52  hours  of  simulation  time  are  summarized  in  this 
section.  These  results  are  presented  in  the  form  of  time 
histories  of  the  estimation  error  as  well  as  the  tracking 
filter  standard  deviation  boundary  for  an  error  state.  The 
standard  deviation  boundaries  are  computed  by  taking  the 
square  root  of  the  diagonal  element  of  the  covariance  matrix 
corresponding  to  the  particular  error  of  interest  and  plot- 
ting plus-and-minus  the  value.  Theory  says  (Ref.  5)  that  if 
everything  is  modeled  properly,  the  estimation  error  should 
not  go  outside  that  standard  deviation  boundary  more  than 
31.7%  of  the  time.  If  the  filter  estimate  is  seen  to  diverge 
outside  the  boundary,  then  the  filter  does  not  properly  model 
the  physical  process,  and  a filter  redesign  is  in  order. 

Typical  performance  of  the  tracking  system  in  esti- 
mating satellite  ephemerides  is  illustrated  in  Figs.  B. 3-1  and 
B.3-2.  These  show  the  error  in  estimating  the  Z-component^ 

♦This  value  assumes  prefiltering  of  the  tracking  data  over 
the  15  minute  update  intervals. 

+The  Z-component  lies  along  the  earth's  polar  axis. 
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of  the  position  and  velocity  of  Satellite  No.  4.  Satellite 
No.  4 is  tracked  for  the  first  five  hours  after  start  of 
tracking  and  then  again  from  the  14^^  and  16th  hour  after 
the  start  of  tracking.  These  tracking  periods  are  repeated 
every  24  hours.  After  the  second  period  of  tracking,  the 
estimation  errors  are  well  within  the  standard  deviation 
boundary,  indicating  that  the  filter  has  estimated  the 
satellite  position  and  velocity  accurately.  Similar  behavior 
was  noted  for  the  other  two  error  components.  The  behavior 
of  the  filter  in  estimating  the  position  and  velocity  of 
Satellite  No.  4 after  the  initial  transient  has  died  out  is 
illustrated  more  clearly  in  Figs.  B.3-3  and  B.3-4,  respec- 
tively. These  results  indicate  that  the  steady  state  per- 
formance of  the  tracking  filter  is  well  behaved  and  thus  the 
errors  are  properly  modeled. 

The  errors  in  estimating  the  states  of  the  clock 
in  Satellite  No.  4 are  illustrated  in  Figs.  B.3-5,  B.3-6, 
and  B.3-7  for  the  clock  phase,  frequency  and  flicker  noise, 
respectively.  The  clock  phase  and  frequency  error  appear 
to  become  stabilized  after  the  second  period  of  tracking  as 
was  evident  for  the  satellite  position  and  velocity  errors. 
The  flicker  noise  error  does  not  become  stabilized  since 
this  state  is  driven  by  random  disturbances  or  noise.  How- 
ever, the  flicker  noise  error  stays  mostly  within  the  fil- 
ter's standard  deviation  boundary  which  indicates  that  the 
error  is  bounded  and  thus  properly  modeled  in  the  filter. 

The  uncertainty  in  estimating  the  solar  radiation 
force  on  Satellite  No.  4 is  illustrated  in  Fig.  B.3-8. 

Solar  radiation  force  uncertainty  is  modeled  as  a random 
Markov  process  and  can  only  be  estimated  during  periods  of 
visibility  of  the  satellite  from  at  least  one  tracking  sta- 
tion. Since  the  initial  uncertainty  was  2%,  further  reduc- 
tions will  not  be  noticed  in  these  results. 
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Figure  B.3-5  Satellite  No 


Figure  B.3-6  Satellite  No. 
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Figure  B. 3-7  Satellite  No.  4 Clock 

Flicker  Noise  Error 


R— 16601 


Figure  B.3-8  Satellite  No.  4 Solar  Radiation 

Force  Uncertainty 
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The  uncertainty  in  the  ground  tracking  station  loca- 
tions is  modeled  in  the  tracking  filter  as  a bias.  Typical 
performance  of  the  tracking  filter  in  estimating  station  loca- 
tions is  illustrated  in  Fig.  B.3-9,  which  shows  the  error  in 
estimating  the  Z-component  of  Monitor  Station  No.  3 location. 
Similar  performance  results  occur  for  the  other  components 
of  the  location  error  and  for  the  other  tracking  stations. 


Figure  B.3-9  Monitor  Station  No.  3 Location  Error 

The  clocks  in  the  monitor  stations  are  the  same  as 
the  satellite  clocks  and  exhibit  similar  behavior  as  illus- 
trated in  Figs.  B.3-10,  B.3-11,  and  B.3-12  for  Monitor  Sta- 
tion No.  3.  The  estimation  errors  tend  to  be  within  the 
standard  deviation  boundaries  thus  indicating  that  the  fil- 
ter is  doing  a good  job  of  estimating  the  monitor  station 
clock  errors.  Similar  behavior  is  evident  with  the  other 
monitor  station  clocks. 
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Figure  B.3-10  Monitor  Static 

Clock  Phase  Ei 
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Figure  B.3-12  Monitor  Station  No.  3 

Clock  Flicker  Noise  Error 

The  gravitational  uncertainty  influences  all  satel- 
• lite  orbits  and  is  modeled  as  a bias.  The  error  in  estimat- 
ing this  uncertainty  is  illustrated  in  Fig.  B.3-13.  Because 
of  the  small  initial  uncertainty  in  this  parameter,  the  fil- 
ter does  not  significantly  reduce  the  uncertainty  during  the 
simulation  tracking  period. 

The  performance  of  the  tracking  filter  after  49 
hours  of  tracking  is  summarized  in  Table  B.3-1.  The  initial 
errors  and  their  standard  deviations  in  addition  to  the  er- 
rors and  filter  statistics  at  49  hours,  are  presented  for 
comparison  purposes.  These  results  confirm  those  presented 
previously  in  this  section  in  that  there  is  no  evidence  of 
filter  divergence  indicating  the  tracking  filter  is  operat- 
ing properly. 

The  satellite  position  and  velocity  (Z-component ) 
data  transmitted  by  Satellite  No.  4 to  the  users  is 
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Figure  B.3-13  Gravit?tional  Constant  Uncertainty 

illustrated  in  Figs.  B.3-14  and  B.3-15,  respectively.  Also 
shown  in  these  figures  is  the  corresponding  estimation  er- 
ror of  the  tracking  filter.  Since  this  satellite  is  up- 
loaded with  data  at  the  105  min  point  each  day,  the  satel- 
lite data  is  effectively  reset  to  correspond  to  that  of  the 
tracking  filter  at  that  time.  However,  because  there  are 
no  measurement  updates  associated  with  the  satellite  data, 
it  tends  to  diverge  slowly  from  the  tracking  filter  esti- 
mates after  the  uploading. 
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Figure  B.3-14  Ground  System  Estimation  and  Satellite 

Prediction  Errors  of  Satellite  No.  4 
Z Position 


Figure  B.3-15  Ground  System  Estimation  and  Satellite 

Prediction  Errors  of  Satellite  No.  4 
Z Velocity 
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APPENDIX  C 


CLOCK  ERROR  MODELING 


This  appendix  will  summarize  the  basic  approach  taken 
in  the  modeling  of  the  intrinsic  stability  of  a time  standard 
(clock);  a more  extensive  treatment  of  this  topic  may  be 
found  in  Ref.  3.  Intrinsic  stability  is  the  characteristic 
behavior  of  the  clock  when  it  is  operating  in  a benign  en- 
vironment. Thus,  the  model  for  intrinsic  stability  will  pro- 
vide a lower  bound  on  the  clock  errors.  Operation  in  a non- 
ideal environment  will  degrade  the  performance  of  the  clock. 

A detailed  discussion  of  environmentally- induced  errors  and 
their  modeling  is  presented  in  Section  4.1. 


C.l  FREQUENCY  STABILITY 

A clock  or  time  standard  is  basically  a frequency 
source  coupled  with  a counter  that  accumulates  the  total 
elapsed  phase  as  depicted  in  Fig.  C.l-1.  The  performance  of 
the  time  standard  is  directly  related  to  the  stability  of 
the  frequency  source.  Thus  it  is  common  to  deal  with  fre- 
quency stability  when  discussing  the  performance  of  time 
standards.  The  following  will  consider  a frequency  source 
whose  instantaneous  output  may  be  written  as 

V(t)  = Vo  sin  [2TrfQt  + 6<Kt)J  (C.l-1) 

where  it  is  assumed  that  the  frequency  fluctuations,  6<j>(t) 
are  very  small,  i.e., 
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Figure  C.l-1  Relationship  Between  a Frequency 

Standard  and  a Clock 
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(C. 1-2) 


o 


The  instantaneous  phase,  <j>(t),  is 


4>(t)  = 2ir  fQ  t + 54>(t) 


(C.l-3) 


w 


ith  the  instantaneous  frequency 


f(t)  - fo  + 2F  5*(t) 


(C.l-4) 


The  quantity  of  interest  is  the  fractional  frequency  fluc- 


tuation y(t)  = Af / f o defined  as 


y(  t) 


f(t)  - f 


o 6 <t>  ( t ) 

f 2 7T  f 

o o 


(C . 1-5 ) 


Two  measures  of  frequency  stability  are  now  pre- 
sented, one  in  the  frequency  domain  and  the  other  in  the 
time  domain.  The  frequency  domain  measure  of  frequency 
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stability  is  the  two-sided  power  spectral  density  of  the 
fractional  frequency  fluctuations,  Sy(u)).  The  power  spec- 
tral density  may  be  defined  as  the  Fourier  transform  of  the 
autocorrelation  function  of  y(t)  if  it  is  assumed  that  y(t) 
is  stationary.  This  assumption  has  been  verified  in  prac- 
tice, hence 

00 

/-jw£ 

RyU)  e dC  (C.  1-6) 


where 


T 

RV(C)  = Lim  f y(t)  y(t+c)  dt 

y t -*•  <»  J 

-T 


(C.  1-7) 


is  the  autocorrelation  function  of  y(t).  (It  has  been  im- 
plicitly assumed  that  both  stationarity  and  ergodicity  hold 
for  Eq.  (C.l-7)  to  be  valid.)  Equation  (C.l-7)  may  also  be 
written  symbolically  as 

Ry(0  = (y(t)  y(t+o)  (C.l-8) 

with  <•  > denoting  infinite  time  average. 


The  time  domain  measure  of  frequency  stability  most 
commonly  used  is  the  two-sample  Allan  variance  defined  as 


°y<T>  5 


with 


(C. 1-9) 


y 


k 


dt , 


tj^+i  = » k— 0 ,1,2,  . . . 

(C.l-10) 
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being  the  average  fractional  frequency  fluctuation  over  the 
averaging  interval 

It  has  been  shown  (Ref.  9)  that  a relationship  ex- 
o 

ists  between  S (w)  and  o (t)  for  a certain  class  of  fre- 

y y 

quently  encountered  noise  processes.  If  Sy(uj)  is  restricted 
to  a power  law  spectral  density,  i.e., 

2 

Sy(»)  = 2 hn  ^n  0 < 0)  ^ 0)h  (C.l-11) 

n=-l 

2 

then  conversions  between  Sy(oj)  and  ay(T)  are  possible  with 
the  aid  of  the  expressions  given  in  Table  C.l-1.  The  physi- 
cal characteristics  of  the  power  law  spectral  density  terms 
(with  -1  <.  n <.  2)  are  indicated  in  the  table. 

TABLE  C.l-1 

RELATIONSHIPS  BETWEEN  TIME  AND 
FREQUENCY  DOMAIN  STABILITY  MEASURES 


These  relationships  between  the  time  and  frequency 
domain  stability  measures  have  direct  application  in  time 


TYPE  OF  NOISE 
PERTURBING  THE 
OSCILLATOR 

SPECTRAL  DENSITY  OF 
FREQUENCY  FLUCTUATIONS 

sy(«) 

TWO-SAMPLE 
ALLAN  VARIANCE 

°y(*> 

White  Phase  Noise 

. 2 
h2  “ 

h2  6 wh 

Flicker  Phase  Noise 

hjM 

h. 

[9  ♦ 6 tn(u>hT) 

T - 2 in  2]  : 

White  Frequency 
Noise 

ho 

VT 

Flicker  Frequency 
Noise 

h.i/M 

h_x  4 in  2 

V >>  1 
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standard  error  modeling.  When  time  standard  error  models 


are  expressed  in  state-space  notation,  the  process  noise  in 
the  model  (see  Table  C.l-1)  is  characterized  by  its  spectral 
level;  i.e.,  the  appropriate  h-coef f icient . For  typical  time 
standard  applications,  the  low  frequency  terms  (with  n £ 0) 


dominate  and  time  domain  stability  measures  (i.e.,  the  Allan 

2 

variance)  are  used  exclusively.  Hence,  knowledge  of  a^(T) 
for  differing  values  of  t affords  an  estimate  of  the  h^. 


The  typical  method  of  presenting  time  domain  sta- 
bility data  is  via  a fractional  frequency  stability  or 
"sigma-tau"  plot.  Figure  C.l-2  is  typical  and  will  be 
utilized  to  obtain  error  model  parameters  for  the  example 


quartz  crystal  oscillator,  the  Hewlett  Packard  (HP)  10544A. 
From  the  figure,  it  can  be  seen  that  flicker  frequency 
noise(ay(x)  ~ constant)  is  the  dominant  noise  term  for  this 
quartz  oscillator.  The  portion  of  the  quartz  curve  with 
Oy(x)~  t is  indicative  of  "aging"  in  the  device. 

The  preceding  discussion  is  only  a brief  outline 
of  frequency  stability;  a more  complete  discussion  may  be 
found  in  Refs.  3 and  9. 

C . 2 FLICKER  NOISE  MODEL 

As  noted  in  Section  C.l,  flicker  noise  is  the  dom- 
inant error  source  for  the  quartz  crystal  oscillator  under 
consideration.  This  section  will  outline  the  development 
of  a model  for  flicker  noise;  Ref.  3 contains  a more  com- 
prehensive derivation.  The  flicker  noise  model  to  be  de- 
veloped will  be  linear,  for  ease  of  inclusion  in  a Kalman 
filter  estimation  algorithm  or  associated  error  covariance 
analysis.  It  will  also  be  an  approximation  since  the  ob- 
served power  spectral  density  (l/|w|)  cannot  be  realized 
exactly  by  a f jnHte  dimensional  linear  system.  (Only  power 
spectral  densitites  which  vary  as  even  powers  of  w can  be 
treated  exactly.) 

Since  its  most  natural  description  is  in  the  fre- 
quency domain,  the  methodology  outlined  below  will  be  to 
attempt  to  approximate  the  power  spectral  density  of  flicker 
noise  (Sy(oj)  = h ^/|u>|)  as  closely  as  possible  subject  to 
the  constraint  that  only  a finite  dimensional  linear  sys- 
tem may  be  used.  The  first  step  in  the  approximation  is 
to  pass  a known  random  process  through  a linear  frequency 
domain  "shaping  filter"  to  obtain  the  desired  random  pro- 
cess (i.e.,  flicker  noise).  This  requires  the  application 
of  a fundamental  result  of  linear  system  theory. 
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Consider  a time  invariant  linear  system  with  out- 
put y(t)  and  input  w(t),  which  is  a sample  function  of  a 
known  random  process.  It  can  be  shown  that  the  power  spec- 
tral densities  of  the  input  and  output  (Sw(to)  and  Sy(co), 
respectively)  are  related  by  Eq.  (C.2-1), 

sy  (03)  = I G ( j 03  ) ] 2 Sw(03)  (C.2-1) 

where  G(jw)  is  the  system  transfer  function. 


In  this  discussion  it  will  be  assumed  that  the  in- 
put random  process  is  white  noise  and  the  desired  output 
process  should  approximate  flicker  noise  Ti.e.,  S (w)  = 
q'  = constant  and  Sy(w)  ~ h ^/lul  over  several  decades  of 
angular  frequency,  w].  From  the  above  it  is  evident  that 
the  flicker  noise  model  transfer  function  Gp(jw)  should 
sat isf y : 


GF(  J03)  I 2 = h_1/(q'|a)|) 


(C.2-2) 


Equation  (C.2-2)  provides  the  characteristics  of 
the  shaping  filter  required  to  continue  with  the  develop- 
ment of  the  flicker  noise  model.  Since  Eq.  (C.2-2)  cannot 
be  realized  exactly  (it  is  not  an  even  power  of  co),  an  ap- 
proximation is  required.  The  approximation  to  Eq.  (C.2-2) 
proposed  by  Barnes  and  Jarvis  (Ref.  10)  is  now  given. 


Consider  approximately  Eq.  (C.2-2)  by  the  trans- 
fer function 


GF(ju))  = Go 


N 

TT 

i.=l 


ju>  + a. 


j ( u>/ 6 ) + ai 


(C.2-3) 


where  a . . 

l+l 


*±l&‘ 


and  6 is  a constant  less  than  unity.  A 
log-log  plot  of  Eq.  (C.2-3)  for  the  N-stage  flicker  noise 
model  transfer  function  is  given  in  Fig.  C.2-1. 
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Figure  C.2-1  General  N-Stage  Flicker  Noise 

Model  Transfer  Function 


The  transfer  function  of  Eq . (C.2-3)  has  four  dis- 
tinct parameters  which  may  be  varied  as  requirements  die- 
tage.  locates  the  lower  end  of  the  frequency  interval 

(for  a fixed  B)  over  which  the  approximation  is  valid.  This 
can  be  seen  by  reference  to  Fig.  C.2-1  in  which  the  approxi- 
mating interval  is 


to, 


03/2  „ , 02N-3/2 

= a.  B ' < (jo  < a.,  / B = u). 

1 l h 


(C.2-4) 


The  width  of  the  approximation  interval  varies  as  B 

2 

while  the  precision  of  the  approximation  increases  as  B 
tends  toward  unity.  Gq  adjusts  the  magnitude  of  the  trans- 
fer function  to  provide  the  correct  flicker  noise  spectral 
level.  In  any  application,  the  choice  of  values  for  a^,  B, 
N and  G^  will  characterize  the  flicker  noise  model. 
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For  error  modeling  purposes,  a time  domain  model 
for  flicker  noise  is  required.  The  transfer  function  of 
Eq.  (C.2-3)  can  be  associated  with  the  differential  system 

(aN  DN  + ...  + aQ) y = (bN  DN  + ...  + bQ) w 

(C. 2-5) 

where  D1  is  the  differential  operator,  Di(-)  = d*(.)/dt*. 
Equating  the  input-output  transfer  function  of  the  system 
defined  by  Eq.  (C.2-5)  with  the  flicker  noise  model  trans- 
fer function  (setting  s = ju>)  yields 


Gf(s)  = 


. N 
bN  S + 


+ b. 


N 

n 


s + a 


aN  sN  + ...  + aQ  !■!  (s'6>  + “i 


(C.2-6) 


The  coefficients  bi  and  a^  can  be  obtained  by  equating 


N 

bN  sN  + . . . + bQ  = IT  (s  + a ) (C.2-7) 

i=l 


and  setting 

ai  = bi/01 

The  constant  parameter  Gq  has 
term,  w,  for  this  discussion, 

EP]  E qf  = q " i gq  i 


(C.2-8) 

been  absorbed  into  the  input 
that  is, 

(C. 2-9) 


Zadeh  and  Desoer  (Ref.  11)  give  the  realization 
shown  in  Fig.  C.2-2  for  Eq.  (C.2-5).  Typically,  this  time 
domain  realization  of  the  flicker  noise  approximation  will 
be  used  as  a noise  input  at  the  frequency  level  in  a time 
standard  error  model.  In  this  case,  the  output  y is  the 


155 


Figure  C.2-2 


Time  Domain  Realization  of  Flicker 
Noise  Transfer  Function 


input  to  the  phase  error  (6<j>)  integrator  or 

6*<t>  = y 

The  state  equations  can  then  be  written  (with  reference  to 
Fig.  C.2-2)  as 


x = F x + £ w 


(C . 2-10) 


with 


(C.2-11) 


0 l/aN  0 0 

0 -aN-l^aN  1 ° 

0 "ftN-l/aN  ° 1 


(C.2-12) 


and 


N 


aN  bN-l  " aN-l  bN 


aN  bN-2  " aN-2  bN 


S = 


N 


aN  b0  " a0  bN 


(C.2-13) 


C . 3 MODELING  OF  INTRINSIC  STABILITY 

Intrinsic  stability  will  be  treated  as  the  perform- 
ance of  the  oscillator  in  an  environmentally  favorable  (be- 
nign) situation.  In  practice,  it  is  common  to  maintain  as 
favorable  an  environment  as  possible  when  measuring  the 
fractional  frequency  stability  of  an  oscillator  (see  Fig. 
C.l-2).  The  two  sources  of  error  which  arise  for  a quartz 
crystal  oscillator  in  a benign  environment  are  flicker  noise 
and  aging. 

Aging  is  a nearly  linear  change  in  resonance  fre- 
quency with  time.  The  a priori  magnitude  of  this  frequency 
drift  cannot  be  determined.  However,  by  observing  its 
effects  over  a long  period  of  time  (days),  the  aging  rate 
can  be  measured.  This  behavior  can  be  adequately  modeled 

as  a random  bias.  The  maximum  (3a)  aging  rate  for  the 

-15 

HP  10544A  is  stated  as  5.79  x 10  /sec. 

In  Section  C.2,  a model  for  flicker  noise  was  de- 
veloped. The  input  to  the  model  was  white  noise  with  a 
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spectral  level,  q^,  which  is  dependent  upon  the  parameters  of 
the  model.  If  Gq  in  Eq.  (C.2-3)  is  chosen  such  that 

I Gp(  ju> ) | 2 ~ 1/M  (C.3-1) 

then  q"  = h ^ [from  Eq.  (C.2-2)]  and  q^  can  be  evaluated  as 

qf  = h-liGo|2  (C. 3-2) 


In  Ref.  3 it  is  shown  that  a choice  of  N=2, 

-2  -1 

a1=10  and  3=10  results  in  a flicker  noise  model  which 
is  valid  over  the  approximation  interval 

2 sec  £ x < 2 x 104  sec  (C.3-3) 

This  interval  is  compatible  with  the  observed  flicker  noise 
region.  For  Eq.  (C.3-1)  to  be  satisfied,  with  this  choice 
of  N,  and  3,  requires  that  |Gq|  = 3.16  x 10  . Hence, 
the  spectral  level  of  the  white  noise  driving  the  flicker 
noise  model  is  calculated  as 

qf  = 3.16  x 103  h_1  (C.3-4) 


A value  for  h_1  may  be  obtained  from  Table  C.l-1  and  Fig. 
C.l-2  as 


o2(t ) 

. y 

4 In  2 


ii 


x 10 


-11 


2.77 


3.66  x 10"23 


(C . 3-5 ) 


The  state  equations  of  the  flicker  noise  model  for 
the  HP  10544A  crystal  oscillator  are  described  by  the  state 
vector , 


x = 


(C.3-6) 


and  the  system  matrices 
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F 


1 


(C.3-7) 


0 

0 


oj 


where 

f1  = 1.0  x 10-2 

f2  = -1.01  x 10-1  sec- 

f3  = -1.0  x 10”4  sec-2 

and 


£ = 


ei 

^2 

L K3 


where 

gl  = 1.0  x 10-2 

g2  = 9.09  x 10-1  sec-1 

g3  = 9.90  x 10“3  sec-2 

The  spectral  level  of  the  white  driving  noise  is 
q-  = 1.14  x 10"19 


(C . 3-8) 


(C . 3-9) 


APPENDIX  D 


PROPAGATION  PATH  MODELS 


The  propagation  of  radio  frequency  energy  from  a 
transmitter  to  a receiver  both  located  close  to  the  earth's 
surface  is  subject  to  many  disturbances,  most  significant  of 
which  are  the  earth's  surface  and  atmosphere.  There  are  a 
number  of  mechanisms  by  which  propagation  of  radio  frequency 
energy  can  take  place,  but  for  the  frequencies  of  interest 
in  this  report  (L-band,  1.6  GHz),  the  primary  mechanism  is 
tropospheric  waves. 

Energy  propagated  via  tropospheric  waves  travels 
between  the  transmitter  and  receiver  within  the  earth's 
troposphere,  which  is  that  portion  of  the  earth's  atmosphere 
within  15  km  of  the  surface.  There  are  several  paths  that 
this  wave  can  take  between  the  transmitter  and  receiver, 
depending  upon  their  separation.  These  paths  lie  in  two 
regions:  within  line-of-sight  and  beyond  line-of-sight . 

The  energy  loss  in  the  propagation  of  radio  waves 
is  the  subject  of  this  appendix.  In  particular,  models  of 
energy  loss  along  the  various  propagation  paths  (path  los- 
ses) that  were  used  in  the  development  of  the  results  in 
Section  3.1.3  are  discussed.  Models  of  these  losses  within 
line-of-sight  are  discussed  in  Section  D.l  while  those  for 
beyond  line-of-sight  are  discussed  in  Section  D.2.  In 
addition,  a model  of  path  loss  common  to  both  regions  is 
discussed  in  Section  D.3.  A typical  example  of  propagation 
path  energy  loss  using  these  models  is  presented  in  Section 
D.4. 


D.l  WITHIN  LINE-OF-SIGHT 


If  the  transmitter  and  receiver  are  within  line-of- 
sight  of  each  other,  then  the  radio  wave  can  travel  two  oaths. 
One  is  the  direct  path  between  the  transmitter  and  the  re- 
ceiver, while  the  other  is  a result  of  reflection  from  the 
earth's  surface.  Models  of  these  two  paths  are  presented 
below. 


D.1.1  Direct  Path 


The  loss  in  energy  through  the  direct  path  is  due 
to  spreading  of  the  energy  over  an  ever-increasing  (with 
distance)  spherical  surface  area.  This  is  also  called  the 
inverse  square  distance  loss  or  free  space  loss.  This  loss  is 
if  modeled  as  follows  (Ref.  15): 


(D.l-1) 


where 


Ld  is  the  direct  path  loss 

X is  the  wavelength  of  radio  wave 

d is  the  distance  between  the  transmitter 
and  receiver  antennas  in  same  units  as  X 

g.  is  the  transmitter  antenna  gain  in  direc- 
tion of  receiver 

g is  the  receiver  antenna  gain  in  direc- 
r tion  of  transmitter 


This  is  the  only  loss  that  would  be  present  in  free  space, 
where  there  are  no  objects  to  reflect  the  radio  wave  or 
otherwise  influence  it.  This  is  the  basic  loss  to  which  the 
other  losses  due  to  reflections,  etc.,  discussed  below  are 
referenced. 
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D.1.2  Reflected  Path 


When  either  the  transmitter  or  receiver  is  close  to 
the  earth's  surface,  the  wave  will  be  reflected  by  the  earth 
as  illustrated  in  Fig.  D.l-1.  The  loss,  Lr , in  this  reflected 
path  relative  to  the  free  space  loss  is  modeled  as  follows 
(Ref.  29): 


where 

Ik  is  the  effective  height  of  the  trans- 
t mitter  antenna  above  the  reflecting 
plane. 

h'  is  the  effective  height  of  the  receiver 
r antenna  above  the  reflecting  plane. 

This  model  assumes  the  earth  has  a perfectly  smooth  surface 
for  reflection.  The  reflecting  plane  is  tangent  to  the 
earth's  surface  at  the  point  of  reflection.  It  is  important 
to  note  the  distinction  between  effective  antenna  heights 
used  in  Eq.  (D.l-2)  and  actual  values.  Near  the  horizon, 
the  actual  values  may  differ  significantly  from  the  effec- 
tive values. 


2 sin 


2,  h;  h- 

i Xd 


(D.l-2) 


Figure  D.l-1  Line-of-Sight  Propagation 
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The  path  of  a radio  wave  is  bent  as  it  passes  through 
the  atmosphere.  In  order  to  be  able  to  represent  these  paths 
as  straight  lines,  an  "effective  earth's  radius"  is  used  in 
the  computations  of  effective  antenna  heights.  This  effec- 
tive earth's  radius  is  a function  of  antenna  height  and  the 
earth's  surface  refractivity  (zero  altitude)  as  follows  (Ref. 
30) : 

a = a [l  + J^SL  Ml hi]  1 
a ao  L1  + n(h)  3h  J 

n (h)  = 1 + N exp  [-  0.139  h] 

s 

where 

a is  the  effective  earth  radius  in  km 

a,_  is  the  actual  radius  of  the  earth  in  km 

° —6 
Ng  is  the  earth's  surface  refractivity  (300x10  ) 

h is  the  height  of  the  lowest  antenna  in  km 

The  reflected  wave  causes  reinforcement  of  the  direct 
wave,  accounting  for  the  factor  of  two  in  Eq.  (D.l-2).  Par- 
tial cancellation  may  also  occur,  especially  where  the  sep- 
aration between  transmitter  and  receiver  antennas  in  relation 
to  their  heights,  is  not  large.  These  effects  will  be  illus- 
trated in  Section  D.4. 


(D.l-3) 

(D.l-4) 


D.  2 BEYOND  LINE-OF-SIGHT 

If  the  transmitter  and  receiver  are  beyond  the  line- 
of-sight  relative  to  one  another,  the  direct  and  reflected 
waves  do  not  exist.  However,  some  energy  is  propagated  be- 
yond line-of-sight  by  the  phenomenon  of  diffraction  and  by 
scattering  of  the  radio  waves  in  the  upper  troposphere,  usu- 
ally identified  as  tropospheric  scatter  or  forward  scatter. 
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Models  of  the  losses  occurring  in  these  two  paths  are  pre- 
sented in  the  sections  below. 

D.2.1  Diffraction  Loss 


Radio  waves  are  bent  around  the  earth  by  the  phen- 
omenon of  diffraction.  This  phenomenon  occurs  when  the  di- 
rect path  of  the  radio  wave  approaches  tangency  to  the  earth's 
surface.  The  effect  of  diffraction  around  the  earth's  curva- 
ture is  to  make  possible  transmission  of  radio  waves  to  well 
beyond  the  horizon,  but  with  an  additional  loss  relative  to 
free-space  loss.  This  additional  loss  or  attenuation,  called 
diffraction  loss,  depends  upon  the  separation  and  heights  of 
transmitter  and  receiver  antennas,  electromagnetic  ground 
constants,  earth  radius,  the  radio  frequency  and  polariza- 
tion of  the  radio  wave.  No  simple  relation  exists  which 
expresses  diffraction  loss  in  these  various  factors.  Avail- 
able data  is  empirical  (Ref.  29).  The  diffraction  loss  com- 
putations for  this  study  assumed  a smooth  earth  surface, 
horizontal  wave  polarization  and  average  earth  surface 
properties  (i.e.,  e = 15,  a = 0.005  mhos/m).  Diffraction 
loss  was  assumed  to  start  at  the  line-of-sight  separation, 
which  is  standard  practice  in  propagation  loss  calculations, 
providing  a smooth  transition  across  the  line-of-sight 
boundary . 

D.2.2  Tropospheric  Scatter 

Well  beyond  the  horizon,  a propagation  path  between 
a transmitter  and  receiver  can  exist  due  to  scattering  of 
the  radio  wave  by  atmospheric  turbulence  or  incoherent  re- 
flections in  the  upper  troposphere.  This  propagation  mode 
is  usually  termed  forward  scatter  or  tropospheric  scatter. 

The  quantity  of  energy  propagated  in  this  mode  is  small 


relative  to  free-space  propagation,  and  subject  to  large 
variations  with  time.  Again,  no  simple  equation  relating 
this  forward  scatter  characteristic*  to  the  various  factors 
of  influence  is  available.  However,  a long-term  median  trans- 
mission loss  due  to  forward  scatter  can  be  computed  based  on 
empirical  data  (Ref.  29). 

For  propagation  paths  extending  only  very  slightly 
beyond  line-of-sight , diffraction  will  be  the  dominant  propa- 
gation mechanism  and  scattering  can  be  neglected.  Conversely, 
for  long  paths,  the  diffraction  field  may  be  much  weaker  than 
the  scattered  field,  and  thus  the  diffraction  mechanism  can 
be  neglected.  In  intermediate  cases,  both  mechanisms  have 
to  be  considered  and  the  results  combined  as  discussed  in 
Ref.  29. 


D. 3 ATMOSPHERIC  ATTENUATION 


At  the  GPS  navigation  signal  frequencies  (L-band, 
1.6  Grfz),  attenuation  of  the  radio  waves  due  to  absorption 
or  scattering  by  constituents  of  the  atmosphere  and  by  par- 
ticles in  the  atmosphere  may  be  significant.  At  these  fre- 
quencies, the  attenuation  is  due  mainly  to  water  vapor  and 
oxygen.  This  atmospheric  absorption  attenuation  or  loss 
relative  to  free  space  loss  is  given  by  (Ref.  29) 


L 

a 


1Q0 . 009  d 


(D.3-1) 


where 


La  is  the  atmospheric  loss  relative  to  free  space 

d is  the  distance  between  transmitter  and 
receiver  in  km. 


♦ 

Hereafter  referred  to  as  a loss  relative  to  free  space  loss. 
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This  loss  occurs  throughout  the  troposphere  and  thus  is  in- 
cluded in  all  propagation  path  loss  computations  for  this 
study.  Note  that  at  100  km,  this  atmospheric  attenuation  is 
less  than  1 dB,  but  at  1000  km  it  is  9 dB. 


D . 4 PROPAGATION  PATH  LOSS 

The  total  loss  in  transmitted  energy  of  a radio  wave 
is  the  product  of  the  losses  in  the  various  paths  that  the 
wave  can  follow.  It  is  standard  practice  to  express  this 
loss  in  decibels;  thus,  the  total  loss  is  expressed  as  the 
sum  of  the  individual  path  losses,  each  expressed  in  dB. 

The  free  space  propagation  path  loss  is  the  basic  loss  to 
which  all  other  losses  discussed  in  this  appendix  are  refer- 
enced. That  is,  these  other  losses  (i.e.,  reflected,  dif- 
fraction, scatter  and  atmospheric  attenuation)  occur  in 
addition  to  the  basic  direct  loss . 

Within  line-of-sight , the  total  propagation  path 
loss  in  dB  is 


L - L . + L + L dB  (D.4-1) 

p d r a 

where 

L is  the  total  propagation  path  loss 

p in  dB 

L,  is  direct  wave  or  space  loss  in  dB 
a [Eq.  (D . 1-1 )] 

L is  the  reflected  wave  loss  in  dB 
r [Eq.  (D . 1-2 )] 

L is  the  atmospheric  attenuation  in  dB 
a [Eq.  (D . 3-1 )] 


Beyond  line-of-sight,  the  total  propagation  path  loss  is 


Lp  Ld  + Ldif f + Lfs 


dB 


(D.4-2) 
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where 


Ldiff  *-s  the  diffraction  loss  (relative  to 
free  space  loss)  in  dB 

Lfs  is  the  forward  scatter  loss  (relative 
to  free  space  loss)  in  dB 

As  discussed  in  Section  D.2.2,  the  values  of  the  diffraction 
loss  and  forward  scatter  loss  are  complex  functions  of 
transmitter-receiver  separation,  expressed  empirically. 

A typical  propagation  path  loss  curve  is  illus- 
trated in  Fig.  D.4-1.  The  top  of  this  curve  is  +6  dB  above 
the  free  space  loss,  reflecting  the  factor  of  two  in  Eq. 
(D.l-2);  that  is,  the  reinforcement  of  the  direct  wave  by 
the  reflected  wave.  The  "spikes"  are  the  result  of  the  par- 
tial cancellations  discussed  earlier.  The  propagation  path 
loss  results  presented  in  Section  3.1.3  (i.e.,  Fig.  3.1-1) 
are  worst  case  in  that  the  minimum  loss  is  used,  which  is 
the  top  part  of  the  curve  in  Fig.  D.4-1. 
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